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Abstract 

^D ■ We review the theory of hadronic atoms in QCD + QED, based on 

a non-relativistic effective Lagrangian framework. We first provide an 
introduction to the theory, and then describe several applications: meson- 
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compare the quantum field theory framework used here with the traditional 
approach, which is based on quantum-mechanical potential scattering. 
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1 Introduction 

Consider a hydrogen atom, and replace the electron with a negatively charged 
pion. This compound is called pionic hydrogen and is the simplest example of 
a hadronic atom. In many aspects, the properties of hadronic atoms are similar 
to those of the hydrogen atom, because both bound states are predominantly 
formed by the static Coulomb force. The typical size of pionic hydrogen e.g. is 
characterized by its Bohr radius 

rs = (a/ic)~\ fic = mpM„{mp + A'Uy'^ , (1.1) 

where a ~ 1/137 denotes the fine-structure constant, and where fie is the reduced 
mass of the system. The mass of the proton and of the charged pion are denoted 
by rrip and M^^, respectively. The distance r^ ~ 220 fm is much smaller than the 
hydrogen radius, but still much larger than the range of strong interactions, which 
is typically of the order of a few fm. It is for this reason that strong interactions 
do not change the structure of the bound-state spectrum in a profound manner. 
At leading order in an expansion in a, the energy of 5"- wave states of pionic 
hydrogen is still given by the standard quantum-mechanical formula 

En = mp + M^-^, n = l,2,---. (1.2) 

However, in distinction to ordinary hydrogen, the ground state of pionic hydrogen 
is unstable. It decays mainly via the charge-exchange reaction 7r~p — >■ 7r°n, and 
via the electromagnetic channel 7T~p — > '-/n, with width Fi ~ 1 eV [1]. This 
corresponds to a lifetime ri ~ 10"^^ s, which is much smaller than the lifetime of 
the charged pion, r^r ~ 10"*^ s, so that the pion in the atom can be considered a 
practically stable particle. Despite of its short lifetime, pionic hydrogen can be 
considered a quasi-stable bound state, since the pion travels many times around 
the proton before decaying, as the ratio l/ica^/Fi ~ 10^ indicates. 

To understand the significance of experiments performed with hadronic atoms, 
we imagine a simplified picture, where the interaction between the constituents 
of a hadronic atom consists of only two terms: the long-range Coulomb part 
and the short-range strong interaction. Equation (1.2) describes the bound- 
state spectrum only approximately, at leading order in an expansion in a. The 
leading correction to Eq. (1.2), which emerges at order a^, is due to strong inter- 
actions only - there is no interference between Coulomb and strong interactions 
at this order. Consequently, measuring very accurately the difference between the 
true energy levels of pionic hydrogen and their pure Coulomb values Eq. (1.2), 
one can extract information about the strong interactions between the pion and 
the proton. Because the size of the atom is much larger than the characteris- 
tic radius of strong interactions, the bound-state observables can feel only the 
low-momentum behavior of the strong pion-nucleon S'-matrix - in other words. 



the energy shift must be expressed in terms of the threshold parameters of the 
pion-nucleon scattering amphtude: the scattering lengths, effective ranges, etc. 
Because the characteristic 3-momenta of the hadrons within the atom are of size 
|p^y| ~ r~^ = afXc oc a, the relative strength of the corrections which contain scat- 
tering lengths, effective ranges, ■ ■ ■ should be ordered according to 1 : a^ : ■ ■ ■ . As 
a result of this, the leading-order energy shift depends only on the pion-nucleon 
scattering lengths. The situation closely resembles a well-known example in 
classical electrodynamics, which tells us that an arbitrarily complicated charge 
distribution can be characterized by a few parameters of the multipole expansion, 
if the distance to the system is much larger than the size of the system itself. It 
is also clear, why the observables of the atom cannot depend on anything else 
but on the threshold parameters of the pion-nucleon amplitude: the distances 
characteristic for the atom are already asymptotic for the pion-nucleon system, 
where nothing but the S'-matrix elements can be observed. In other words, the 
measurement of the observables of the pionic hydrogen does not probe the inner 
region of the pion-nucleon interaction. 

The insensitivity of the pionic hydrogen observables to the short-range details 
of the pion-nucleon interaction is very fortunate: it provides us with the possi- 
bility to directly extract the values of the pion-nucleon scattering lengths from 
atomic experiments. A different method for determining the same quantities is 
to measure the scattering cross sections at different energies, and to extrapolate 
the result to threshold [2-7]. The former method, however, is free from the dif- 
ficulties which are related to this extrapolation procedure. This property is even 
more important in other hadronic systems, where the scattering amplitude near 
threshold is hardly accessible by other experimental technique. 

Deser, Goldberger, Baumann and Thirring (DGBT) were the first who derived 
the model-independent relation between the (complex) energy shift of the ground 
state in pionic hydrogen and the strong pion-nucleon scattering amplitude at 
leading order in a [8]. The result reads^ 

AEf --T, = |^io(0)r A(7r-p ^ n-p) + O(a^) , (1.3) 

where |\^io(0)p = 0^111/1: denotes the square of the Coulomb wave function 
of the atom at the origin and is, therefore, a measure of the probability that 
the charged pion and the proton in the atom come very close to each other. 
Further, A{'K~p -^ 7r~p) is the 7r~p elastic scattering amplitude at threshold, 
which describes strong interactions between these particles after they come close. 
The scattering amplitude is normalized so that 

Aiir'p^Tv'p) = a^^ + a^^ + --- . (1.4) 



"'^We use throughout the Landau symbols 0{x) [o{x)] for quantities that vanish Hkc x [faster 
than x] when x tends to zero. Furthermore, it is understood that this holds modulo logarithmic 
terms, i.e. we write also 0{x) for a; In a;. 



Here, the quantities on the r.h.s. denote the S'-wave vrA^ isospin even/odd scat- 
tering lengths, which are defined according to Ref. [2], and the eUipsis stands for 
isospin breaking terms. 

The imaginary part of the quantity y4(7r~j> -^ t^p) is given by unitarity. For 
simphcity, let us assume for a moment that the atom decays only into the strong 
channel 7r~p -^ n^n. Then, the imaginary part is^ 

lmA{7i~p -^ n'p) = 2p*\A{n~p -^ 7r°n)p , (1.5) 

where p* is the CM momentum of the outgoing tt^ in the scattering process 
7i~p — > n^n at threshold. The normalization is chosen such that, in the isospin 
limit, the threshold amplitude is A{7T~p — »■ ir^n) = — ao+. Finally, the decay width 
Fi of the ground state is calculated from Eqs. (1.3, 1.5). Therefore, measuring 
the spectrum of pionic atoms provides information on the pion-nucleon scattering 
lengths. 

Hadronic scattering lengths are important low-energy characteristics of quan- 
tum chromodynamics (QCD) for the following reason. At small momenta of order 
of the pion mass, QCD can be described by a low-energy effective field theory 
(EFT), which contains only hadronic degrees of freedom - the approach goes un- 
der the name chiral perturbation theory (ChPT) [9]. The crucial point is that, in 
certain cases, (e.g., for vrvr scattering [10-15]), it is possible to obtain a very ac- 
curate prediction for the hadronic scattering lengths in ChPT and, consequently, 
comparing the theoretical values with the results obtained in the experiment, 
one can perform a direct study of the properties of the low-energy QCD. For a 
method which does not make use of chiral symmetry to determine the scattering 
lengths, see Refs. [16]. Scattering lengths in the framework of pure QCD can 
now also be determined from numerical simulations of QCD on a lattice [17], see 
also the reviews by Leutwyler [18] and Colangelo [19]. Last but not least, the 
hadronic scattering lengths often serve as an input in the determination of other 
important low-energy characteristics of strong interactions (like the pion-nucleon 
cr-term [20] or the irNN coupling constant [21-23] in the case of the pion-nucleon 
scattering lengths). Below, we discuss the fundamental physics content of various 
hadronic scattering lengths. Here we merely note that the knowledge of the exact 
values of these quantities would greatly advance our knowledge of the fundamen- 
tal features of strong interactions, allowing for a variety of high-precision tests 
of QCD at low energy. 



^Note that, as we shall see below, adding the second decay channel n^p -^■ jn merely leads 
to the replacement hiiA(Tr~p —^ tt~p) ^ (1 + P~^)lTn A('k~ p -^ '^~p), where P — g(t:~p — > 
'K^n)/a{TT~p -^ -fn) is the so-called Panofsky ratio P ~ 1.546, which describes the relative 
probability for the "strong" and "electromagnetic" decay channels. Since the value of the 
Panofsky ratio is taken from other measurements and is thus considered to be an input in the 
analysis of the pionic hydrogen data, this analysis proceeds exactly similar to the 1-channel 
case. 



In order to carry out the comparison of the experimentally determined hadro- 
nic scattering lengths with the theoretical predictions, the accuracy of the lowest- 
order formulae, analogous to Eq. (1.3), is not sufficient in most cases. For this 
reason, there have been numerous attempts to improve the accuracy in these 
relations by including electromagnetic corrections and strong isospin breaking 
effects. These investigations had first been carried out within the framework 
of potential scattering theory: strong interactions are described by potentials, 
assumed to be isospin symmetric. Isospin violation then stems from Coulomb 
corrections, from the mass differences in the free Green functions (masses of all 
particles involved are put equal to their physical masses), from the coupling to 
the electromagnetic channels, etc. Unfortunately, the results of potential-model- 
based calculations turned out to be not unique - in some cases, they even differ 
to an extent that matters at the level of the experimental precision. 

The following two comments are in order at this point. 

i) The scattering lengths which one attempts to extract from the investigation 
of hadronic atoms are the ones in pure QCD, in the isospin symmetry limit 
m„ = rud- In particular, the electromagnetic effects are switched off. In 
potential scattering theory, the Coulomb potential can be easily turned on 
or off. It is however well known that splitting off electromagnetic effects 
in quantum field theory (QFT) is ambiguous and leads in general to scale 
dependent results. 

ii) Whereas the relation of hadronic atom observables to the threshold param- 
eters in the hadron-hadron scattering amplitude is universal, the relation 
of these threshold parameters to the purely hadronic scattering lengths is 
not universal and depends on a complicated interplay of strong and electro- 
magnetic interactions. Potential models do not include the full content of 
isospin breaking effects ab initio and therefore do not, in general, suffice to 
work out this relation. For example, the effects due to the short-range elec- 
tromagnetic interactions and the quark mass dependence of the hadronic 
scattering amplitudes - which naturally emerge in ChPT - are omitted in 
simple hadronic potentials. On the other hand, in many cases the bulk of 
the total correction is provided exactly by these omitted terms. 

These two points illustrate the major challenge to any theoretical investiga- 
tion: in order to be able to fully exploit the high-precision data from present 
and future experiments on hadronic atoms, one needs a systematic framework to 
relate hadronic atom characteristics to observables of the underlying theory. Chi- 
ral perturbation theory is an obvious candidate for such a framework. However, 
despite remarkable progress in describing low-energy scattering processes with 
elementary hadrons in ChPT, the study of hound states within the same theory 
remained terra incognita for a long time. In certain cases, the problem could be 
solved by brute force with the use of methods of QFT, namely, Bethe-Salpeter 



equation or quasipotential reductions thereof [24-30]. On the other hand, the 
matter turned out to be very involved in general - it became obvious that new 
ideas are needed to overcome the difficulties on the way towards a general theory 
of hadronic atoms in QCD+QED. 

A crucial observation was that the calculation of the spectrum of hadronic 
atoms can be treated separately from the chiral expansion [31-36]. The reason 
for this consists in the large difference of the pertinent momentum scales: whereas 
the characteristic momenta in ChPT are of order of the pion mass, the typical 
bound-state momentum is the inverse of the Bohr radius. This difference in 
scales prompts using a non-relativistic effective Lagrangian approach to bound 
states, which has ffist been introduced by Caswell and Lepage in the framework 
of QED [37]. 

Let us shortly describe the pertinent low-energy EFT, which is valid at typical 
momenta of size |p| <^ M^. The structure of this EFT is much simpler than the 
structure of ChPT, for the following reasons, i) As just mentioned, the system 
can be described by means of an effective non-relativistic Lagrangian. Relativis- 
tic corrections are included in a perturbative manner. In bound-states, these 
corrections will obviously contribute at higher orders in a, because p is a quan- 
tity of order a. ii) The bound states in the non-relativistic theory are described 
by a Schrodinger equation, instead of the Bethe-Salpeter equation. This leads 
to dramatic simplifications and allows one to design a systematic perturbative 
approach to the calculation of the bound-state observables on the basis of the 
Feshbach formalism [38,39]. Finally, the parameters of the effective Lagrangian 
are determined by matching this theory to ChPT in the scattering sector. At 
the end of the calculation, the reference to the (auxiliary) non-relativistic EFT 
completely disappears from the final results, which express the observables of 
hadronic atoms it terms of the parameters of ChPT. 

The application of the above simple idea for the description of hadronic atoms 
turned out to be very fruitful. Indeed, within a few years, most of the (theoret- 
ically) relevant hadronic bound states have been revisited using this approach 
(see, e.g., Refs. [36,40-49]. For calculations using field theoretical methods dif- 
ferent from the ones just described, see Refs. [50-62]). Data analysis now rests 
on the basis of low-energy effective field theories of QCD+QED, and is in general 
free of any model assumptions'^. 

The article is organized as follows. In section 2 we start with a detailed 
discussion of the fundamental physics background behind hadronic atom exper- 
iments. In sections 3, 4, 5 and 6, we construct the general theory of hadronic 
atoms in QCD+QED. We keep this part rather self-contained, avoiding extensive 
references to existing literature. For this reason, these sections also contain the 



■^Non-relativistic Lagrangians in the spirit of Ref. [37] are heavily used for the description 
of heavy quark bound states in QCD, see Ref. [63] for the introduction of the method. Two 
recent reviews on the subject are Refs. [64,65]. 



Atom 


Name 


Symbol 


TT + TT^ 


pionium 


A2n 


7r^J^± 


TiK atom 


A^K 


■K~p 


pionic hydrogen 


ttH 


K~p 


kaonic hydrogen 


KH 


TT^d 


pionic deuterium 


ttD 


K-d 


kaonic deuterium 


KD 



Table 1: Hadronic atoms investigated in this report. The corresponding names 
and symbols used in the text are also displayed. 

essentials of non-relativistic effective field theories, to an extent that is needed 
for understanding the material related to hadronic atoms. In section 7 we col- 
lect useful general information on the relation between spectra and scattering 
lengths. We then describe in sections 8-12 several hadronic atoms, which have 
been treated so far in the framework of QFT, see Table 1. 

In section 13 we briefly discuss the relation of our approach to conventional 
potential models. Finally, section 14 contains a summary and concluding remarks. 
Appendix A collects some notation, appendix B discusses the unitarity condition 
for non hermitian Hamiltonians, and in appendix C we discuss matching for 
non-relativistic coupling constants in vrif, using unitarity. 

We note that it was not our aim to cover the large number of articles con- 
cerned with the rich and well-developed phenomenology of hadronic atoms within 
the traditional quantum-mechanical setting, based on hadronic potentials, or to 
provide a complete bibliography on this issue. There are excellent textbooks (see, 
e.g., Refs. [66-68]) or review articles (e.g. Ref. [69]), which should be considered 
as complementary to the present work. In addition, the issue of deeply bound 
exotic atoms lies beyond the scope of the present review and will not be discussed 
at all. In order to follow the developments in this field, we refer the interested 
reader e.g. to the EXA05 proceedings [70]. For a short review, see, e.g., Ref. [71] 
and references therein. 



2 Physics background 

As already mentioned, experiments on hadronic atoms can provide stringent in- 
formation on fundamental properties of QCD. In this section we consider in some 
detail the potential of individual experiments in this respect. 



2.1 Pionium 

Pionium is presently investigated by the DIRAC collaboration at CERN [72-78] . 
The experiment is a very challenging one: the lifetime of pionium in the ground 
state is of the order ri ~ 3 x 10~^^ s. Because the ratio of the binding energy with 
the width is of order 10^, one concludes that pionium is a quasi-stable bound state 
as well. At present, it is not possible to measure this lifetime directly. Instead, 
one first measures the ionization probability of pionium into a charged pion pair 
in different targets. From the calculated relation [79, 80] between this probability 
and the lifetime of pionium (which in vacuum almost exclusively decays into a 
7r°7r° pair), one obtains the lifetime. The DIRAC collaboration has reported [77] 
the result 

n = 2.91 1°J^ (stat) t[Ji^ (syst) x 10^^^ s = 2.91 t^^f^ x 10^^^ s. (2.1) 

Efforts to reduce the experimental uncertainty are ongoing [81]. 

The decay width of pionium in the ground state into a 27r° pair is related to 
the S-wave tttt scattering lengths^ a/ with total isospin J = 0, 2, 

ri,2.o = -a3pt(ao-a2)' + --- , (2.2) 

where p\ = (M^ — M^o — | M^a^)^/^ denotes the magnitude of the 3-momentum 
of the TT^ir^ pair in the final state (with higher order terms in a ignored) , and the 
ellipses stand for the isospin breaking corrections. In the derivation of Eq. (2.2) 
one has not used chiral symmetry: this relation is universal and holds - at leading 
order in isospin breaking - as long as the scale of strong interactions is much 
smaller than characteristic atomic distances. The DIRAC collaboration plans to 
measure Fi 27r" with an accuracy of 10%. Once a reliable evaluation of the isospin 
breaking corrections in the relation Eq. (2.2) is carried out, this measurement 
enables one to determine the value of |ao — 02! with 5% accuracy. As we shall see 
later, these isospin breaking corrections are of the order of 6% and therefore not 
at all negligible at the accuracy one is interested here. 

What makes the above enterprise particularly interesting is the fact that the 
difference 00—02 is very sensitive to the value of the quark condensate in QCD [82], 
and thus to the manner in which chiral symmetry is spontaneously broken. The 
so-called "standard" scenario is based on the assumption that the condensate is 
large. This, for example, signifies that in the expansion of the pion mass in terms 
of the quark mass, 

M^ = M^- 32^ ^' + ^^^'^ , M2 = 2m5 , m = ^(m„ + m^) , (2.3) 

the first term is dominant [83]. Here, F stands for the pion decay constant F^^ 
in the chiral limit^, and I3 denotes a particular low-energy constant (LEC) that 



"'We use a convention where tttt scattering lengths a/ are dimensionless. 
^We use F^=92.4 McV. 
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appears in the O(p^) Lagrangian of ChPT, whereas B is related to the quark 
condensate in the chiral hmit [84,85]. Assuming the standard scenario of chiral 
symmetry breaking, a very accurate prediction of scattering lengths was achieved 
by merging two-loop ChPT with Roy equations [13, 14], 

ao = 0.220 ±0.005, oa = -0.0444 ± 0.0010 , ao - 02 = 0.265 ± 0.004 . (2.4) 

If the experimental value of ao — a2 does not agree with this prediction, this 
would unambiguously indicate that chiral symmetry breaking in QCD proceeds 
differently from the standard picture [82]. As a result of this, the first term in 
the expansion Eq. (2.3) would be non-leading, and the chiral expansion must be 
reordered [82]. 

Measuring the energy levels of pionium [86, 87] enables one to extract a dif- 
ferent combination of the S'-wave scattering lengths. Indeed, the strong energy 
shift of S'-states is given by [44] 

A^f = -- ^(2ao + a2) + ---, (2.5) 

where n denotes the principal quantum number, and where the ellipsis stands 
for isospin breaking corrections. Hence, measuring both, Ti2n0 and AEn, and 
assuming ao — 02 > 0, one may extract the values of oq and 02 separately and 
compare them with the prediction Eq. (2.4). 

The present situation concerning the verification of the prediction Eq. (2.4) 
is the following. Lattice results [17-19], the data from DIRAC [77] on pionium 
lifetime and from NA48 on the cusp in K ^ Stt decays [88-93] neatly confirm 
the prediction, although partly with considerably larger error bars. On the other 
hand, preliminary NA48/2 data on i^e4 apparently preferred larger values [94] of 
Oo than the one displayed in Eq. (2.4). This puzzle was very recently resolved by 
the following observation [95,96]. The electromagnetic corrections applied in the 
NA48/2 Ke4: data were calculated using the PHOTOS Monte Carlo code [97], and 
by applying in addition the Sommerfeld factor [98]. However, this does not yet 
take care of all isospin breaking effects in this decay [99]. Indeed, the kaon can 
first decay into a tt^tt" or ir^rj pair, that then re-scatter into the outgoing charged 
pions, through 7r°7r° —>■ 7r~^TT~ and ir^r] — »■ Tr^TT*^ — > tt+tt". Since the charged pion is 
heavier than the neutral one by about 4.6 MeV, these intermediate states generate 
a cusp in the phase of the relevant form factor. As a result of this, the phase is 
pushed upwards by about half a degree, and does not vanish at the threshold for 
-7r"''7r~ production. Once this is taken into account, the i^e4 data from NA48/2 are 
in nice agreement with Eq. (2.4). There is, however, a discrepancy with the E865 
data [100, 101]: the isospin corrections just mentioned spoil the good agreement 
with the prediction Eq. (2.4). We refer to Refs. [18, 19,95, 102] for more details. 
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2.2 ttK atom 

The DIRAC collaboration at CERN plans [78] to measure the nK atoms in con- 
tinuation of the already running experiment on pionium. From the measurements 
of the energy shift and decay width of the ground state, one would in principle 
extract the two independent isospin combinations a^ and a^ of the S-wave tiK 
scattering lengths, because the strong energy shifts and decay widths of t:K 
atoms in their S— states are related to these scattering lengths in a manner very 
analogous to pionium, 

r. = ^p;(ao)^ + ---- (2.6) 

Here the ellipses stand for isospin breaking effects. 

The theoretical interest in these scattering lengths is twofold at least. First, 
we mention that there exists a low-energy theorem [103], which states that the 
expansion of the isospin odd scattering length a^ is of the form 

^-^^ -^0{Ml)). (2.7) 



StiF^{M^ + Mk) 



Here, M^, Mk and F^^ denote the physical meson masses and the pion decay 
constant in the isospin symmetry limit. The theorem Eq. (2.7) states that the 
corrections to the leading current algebra result are proportional to powers of 
the pion mass - they vanish in the chiral limit m„ = m^ = 0.^ Since there is no 
strong final state interaction in this channel, one expects that these corrections 
are modest. On the other hand, the two-loop calculation performed in Ref. [105] 
suggests that the correction at order p^ is substantial, and moreover larger than 
the one at order p^. The numerical result obtained in Ref. [105] for the scattering 
length ttQ agrees with an analysis [106] based on Roy equations^. 

The situation is puzzling - how can it be that higher orders are larger than 
low-order contributions? As was shown by Schweizer [107], in the present case, 
the large two-loop correction stems mostly from counterterms at order p^, esti- 
mated with resonance saturation in Ref. [105], see also Ref. [108] in this connec- 
tion. Is this resonance estimate correct? If yes, ChPT is turned upside down. If 
no, ChPT does not agree with the result from the Roy equation analysis. 

Second, the scattering length do also depend on certain LECs at order p^, 
whose values would be interesting to know in light of the large/small condensate 



^Implicitly, this result is contained in the explicit expressions of the ttA' amplitudes given 
in Ref. [104]. 

^Note, however, that the available experimental data below 1 GeV are, in general, inconsis- 
tent with the solution of Roy equations, where the data above 1 GeV are used as phenomeno- 
logical input. 
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question [82]. So, in case one does have a precise value of the scattering lengths, 
the first point could be clarified in the sense that one knows what is right and 
what is wrong. Concerning the second point, it remains to be seen whether this 
would indeed allow one to pin down those LECs which play a dominant role 
in connection with the large/small condensate issue. A detailed analysis of this 
point remains to be performed. 

2.3 Pionic hydrogen and pionic deuterium 

Measurements of the energy shift and decay width of the ground state of ttH and 
TiD were carried out by the Pionic Hydrogen Collaboration at PSI [1, 109-114]. 
The pionic hydrogen experiment was upgraded in 1998 [115], followed by an 
increased accuracy in particular of the width measurement [116] (for a review, 
see, e.g., Ref. [117]). The preliminary values of the strong shift and width of the 
ground state of tiH obtained in the upgrade [115] are [118] 

AEf = -7.120 ±0.011 eV, 

Fi = 0.823 ± 0.019 eV. (2.8) 

The announced aim [119] of the collaboration is to extract S-wave vrA^ scattering 
lengths from this experiment with 1% precision. This will be a unique experimen- 
tal result for hadron physics. Of course, in order to achieve this goal by measuring 
TiH alone, one should remove possible sources of theoretical uncertainties related 
to isospin breaking effects, at an accuracy that matches the experimental preci- 
sion. At present, this seems to be a very difficult task, as we will show later in 
this report. 

The measured shift and width of ttD are [112] 

^^str,d ^ 2.46 ± 0.048 eV , F'j^ = 1.194 ± 0.105 eV. (2.9) 

This accuracy is expected to improve in the near future [120, 121]. 

We will show in later sections that, using multiple-scattering theory [22,66, 
68, 122-130] or chiral EFT in the two-nucleon sector [47, 131-143], one can relate 
the value of the vrd scattering length, which can be extracted from the experi- 
mental data, to the vrA^ scattering lengths. This results in additional constraints 
on the values of a^j^ and a,^^_. 

Let us assume that one indeed is able to extract the exact values of the nN 
scattering lengths from the experiment with high precision. What information 
about the fundamental properties of QCD are contained in these precise values? 
Of course, iiN scattering lengths are quantities of fundamental importance in 
low-energy hadronic physics by themselves, since they test the QCD symmetries 
and the exact pattern of the chiral symmetry breaking (see, e.g., Refs. [144-146]). 
Moreover, since the knowledge of these scattering lengths places a constraint on 
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the TiN interactions at low energy, it also affects our understanding of more com- 
plicated systems where vrA^ interaction serves as an input, e.g. A^A^ interaction, 
TT-nucleus scattering, three-nucleon forces, etc. 

In addition to this, the high-precision values of the ttN scattering lengths are 
used as an input for the determination of different basic parameters of QCD at 
low energies more accurately. One example is the irNN coupling constant q^nn, 
which is obtained from the Goldberger-Myazawa-Oehme (GMO) sum rule [21], 
where a particular combination of scattering lengths enters as a subtraction con- 
stant, 



M^ /"°° dk 1 

{(T^-p{k) - a^+p{k)) \ , (2.10) 



iJirNN 

Air 



8vr2 J, ,/Mfr¥ 

where rriN is the nucleon mass in the isospin symmetry limit, cr,r±p(^) denotes the 
total elastic cross section for the scattering of 7r+ (tt^) on the proton, and a,r±p 
stand for isospin combinations of S-wave irN scattering lengths, a^Tp = a^_,_±ao+- 
For recent investigations of the GMO sum rule, see Refs. [22,23]. 

Other important quantities, which can be obtained by using the S-wave nN 
scattering lengths as an input, are the so-called ttN cr-term and the strangeness 
content of the nucleon [20]. The a-term, which measures the explicit breaking 
of chiral symmetry in the one-nucleon sector due to the u- and d-quark masses, 
is given by the matrix element of the chiral symmetry-breaking Hamiltonian 
rh{uu + dd) between the one-nucleon states. The latter quantity is expressed 
through the value of the scalar form factor of the nucleon at t = (p' — p)^ = 0, 

u{p' ,s')aT,N{t)u{p,s) = {p's'\'m{uu + dd)\ps) , o'ttTV = o'ttAtIO) , (2.11) 

where the Dirac spinors are normalized to uu = 2m ]^. The a-term is related to 
the strangeness content of the nucleon y and the S'f/(3) symmetry breaking piece 
of the strong Hamiltonian, 

— {ps\uu + dd — 2ss\ps) = ( -^ — 1 ) (1 — y)a.„N , 



2m7v V "^ 

y= 2(P^I^-^N) ^2.12) 

{ps\uu + dd\ps) 

An analysis of the experimental data, carried out some time ago [147], gives 

(T^jv^45MeV, y~0.2. (2.13) 

Note that in this analysis, one has used S-wave tiN scattering lengths as input 
in the dispersion relations which provide the extrapolation of the isospin even 
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pion-nucleon scattering amplitude from threshold down to the Cheng-Dashen 
point. The a-term, which is obtained as a result of the above analysis, is rather 
sensitive to these scattering lengths, see, e.g., Ref. [148, Fig. 1]. Consequently, 
an accurate measurement of the scattering lengths will have a large impact on 
the values of (jjrAf and y. For an analysis that finds a substantially larger value of 
(jTrAT than the one quoted above, see Pavan et al. [20]. 

2.4 Kaonic hydrogen and kaonic deuterium 

The DEAR collaboration at LNF-INFN [149-152] has performed a measurement 
of the energy level shift and width of the KH ground state with a considerably 
better accuracy than the earlier KpX experiment at KEK [153]. The present 
experimental values of these quantities are [154] 

AEf = 193 ± 37 (stat) ± 6 (syst) eV , 

Fi = 249 ±111 (stat) ±30 (syst) eV. (2.14) 

As can be seen from the above result, the uncertainty is still tens of eV in the 
energy shift and more than 100 eV in the width. Now DEAR is being followed by 
the SIDDHARTA experiment that features new silicon drift detectors [154, 155]. 
The plans of the SIDDHARTA collaboration include the measurement of both, 
the energy shift and width of KH , with a precision of several eV, i.e. at the 
few percent level, by 2008. Moreover, SIDDHARTA will attempt the first ever 
measurement of the energy shift of KD with a comparable accuracy and possibly, 
kaonic helium and sigmonic atoms. Eventually, these experiments will allow one 
to extract the values of the K~p and K~d threshold scattering amplitudes by 
using the pertinent DGBT-type formulae. Because these amplitudes depend 
on the other hand on the S'-wave KN scattering lengths Qq and ai, one can in 
principle determine these. 

The necessity to perform measurements of KD energy shift and width of the 
ground state is also due to the fact that, unlike in the case of pionic atoms, the 
measurement of the KH spectrum alone does not allow one - even in principle - 
to extract independently a^ and ai. This happens because the imaginary parts of 
the threshold amplitudes do not vanish in the isospin limit, being determined by 
the decays into inelastic strong channels ttS, tt'^A, ■ ■ ■ - in other words, the separa- 
tion of thresholds is governed by the breaking of S'f/(3) symmetry. Consequently, 
one attempts here to determine four independent quantities (real and imaginary 
parts of ao and ai) that requires performing four independent measurements - 
e.g., the energy level shifts and widths of KH and KD . Note however that, even 
though it is clear that ao and ai cannot be determined separately without mea- 
suring kaonic deuterium, it is still not evident whether it is possible to do so if one 
performs such a measurement. The reason is that the (complex) kaon-deuteron 
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amplitude at threshold, which is directly determined from the experiment and 
which is expressed in terms of oq and ai through a multiple-scattering series, 
is generally plagued by systematic uncertainties due to a poor knowledge of the 
low-energy kaon-nucleon dynamics. Thus, one needs to know whether these un- 
certainties are small enough not to hinder a determination of oq and ai from the 
SIDDHARTA experiment. 

In order to answer this question, a detailed analysis of the problem has been 
performed in Ref. [49]. This investigation revealed that - at least, within the 
lowest-order approximation - a combined analysis of DEAR/SIDDHARTA data 
on KH and KD turns out to be more restrictive than one would a priori expect, 
see also the discussion in section 12 of this report. The combined analysis im- 
poses stringent constraints on the theoretical description of the kaon-deuteron 
interactions at low energies and provides a tool for determining Oq and Oi with 
reasonable accuracy from the forthcoming SIDDHARTA data. It remains to be 
seen, whether this conclusion stays intact, if higher-order corrections are system- 
atically included. 

What fundamental physical information can be gained from the values of Oq 
and Oi? We believe that it could be very useful to carry out a comparison of the 
KN scattering lengths measured in the DEAR/SIDDHARTA experiment with 
different theoretical predictions based on the unitarization of the lowest order 
ChPT amplitude [156-169]. It turns out that even the data from KH alone im- 
pose rather stringent constraints on the values of the KN scattering lengths - 
namely, in certain cases there emerge difficulties to make DEAR data compatible 
with the scattering sector [46, 160-169]. It is clear that imposing additional con- 
straints from KD data makes the issue of compatibility even more pronounced. 
In our opinion, it will be important to check whether the unitarization approach 
passes this test. 

Finally we wish to note that the original physics program of the DEAR/SID- 
DHARTA experiment was largely a direct extension of the vrif experiment to 
the strange quark sector, including, in addition, the measurement of the kaon- 
nucleon cr-term and the strangeness content of the nucleon [150,151]. On the 
other hand, there are significant differences between these systems [170, 171]: 

i) There are much less data on KN scattering than in the nN case. 

ii) There are open strong channels below threshold. 

iii) There is a sub-threshold resonance A(1405). 

iv) The distance between the threshold and the Cheng-Dashen point is much 
larger than in the ttN case. 

It will therefore be very difficult to extract the experimental value of the KN 
a-term [150,151]. 
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3 Non— relativistic effective theories: strong sec- 
tor 

3.1 Introductory remarks 

In relativistic quantum field theories, the number of particles is not a conserved 
quantity. As a result of this, S'-matrix elements have a complicated analytic struc- 
ture. In particular, they contain branch points at each energy that corresponds 
to a threshold for a newly allowed physical process [172]. 

This feature renders the description of bound states very complicated. On the 
other hand, recall that we are concerned with loosely bound states which gener- 
ate poles in the complex energy plane. The poles are located close to the elastic 
threshold and thus far away from any inelastic and crossed channel singularities. 
Further, the magnitudes of the characteristic momenta are much smaller than the 
particle masses, and the characteristic kinetic energies lie far below the inelastic 
thresholds generated by the presence of additional massive particles in the final 
state. This suggests that a framework where all singularities are treated on an 
equal footing is counterproductive and superfluous in this case: the physics at 
low energies can equally well be described by a simpler effective non-relativistic 
field theory, in which the S'-matrix elements possess only the elastic cut and the 
poles that are located close to this cut, whereas the contributions from distant 
singularities can be Taylor expanded in external momenta [37]. This (truncated) 
Taylor series is generated by a finite number of operators in the non-relativistic 
Lagrangian, multiplied with unknown coupling constants. For properly chosen 
couplings, the non-relativistic and the relativistic theories are physically equiv- 
alent at low energy, up to a certain power in the momentum expansion. This 
procedure to fix the couplings in the non-relativistic Lagrangian is called match- 
ing. 

In this non-relativistic theory, bound-state dynamics is described by the 
Schrodinger equation, and the perturbation expansion of the bound-state en- 
ergies is given by the Rayleigh-Schrodinger series [37]. The crucial property of 
the non-relativistic theory which makes this calculation simple, is the fact that 
the production and annihilation of massive particles can be forbidden by con- 
struction: the interaction Hamiltonian has vanishing matrix elements between 
states with different number of massive particles. In the presence of photons, 
only matrix elements between states that contain the same number of massive 
particles and an arbitrary number of photons can be non- vanishing. 

For the Coulomb bound state, momenta are of order a. As a result of this, the 
contribution to the bound-state energy from an operator in the interaction La- 
grangian which contains n space derivatives is suppressed by a"" as compared to 
the contribution from an operator with the same field content and without deriva- 
tives. Stated differently, one counts each space derivative di = 0{v) = 0{a) (here 
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V stands for the velocity of a massive particle). If one wishes to evaluate the bind- 
ing energy at a given order in a, one may truncate the non-relativistic Lagrangian 
at a certain order in v and carry out calculations with the truncated Lagrangian 
- the error in the bound state occurs at a higher order in a. The method works, 
provided that the power counting stays intact in higher-order calculations in- 
volving loops. This can be achieved by applying threshold expansions, in the 
scattering sector as well as in the Rayleigh-Schrodinger perturbation series (see, 
e.g., Refs. [42,64,173-182]). For a review on electromagnetic bound states see, 
e.g., Refs. [183,184]. 

We illustrate in this and in the following two sections the procedure in a 
simple theory: we consider a self-interacting scalar field which is coupled to 
photons - in the absence of the self interaction, the Lagrangian reduces to the 
one of scalar electrodynamics. We determine the widths and energy levels of 
the bound states at next-to-leading order in a, following closely the procedure 
outlined in Refs. [181,185]. As we shall see later, the evaluation of the energy 
levels and widths in hadronic atoms in the framework of ChPT amounts to a 
rather straightforward generalization of the techniques presented here. 

3.2 The Lagrangian and the reduction formulae 

We start from a situation where only one hard scale is present, and where the 
external momenta involved are much smaller than this scale. This situation is 
described e.g. by a theory that contains solely massive charged scalar particles (p^ 
of physical mass M, and where the external momenta are considered to be small 
in comparison to M. In the relativistic case, we introduce a complex massive 
self-interacting scalar field 0. Photons will be included in the following section. 
The Lagrangian is given by 

C = d^4>{d^'4>y - M200t + ^ (^00t)2 ^ counterterms , (3.1) 

where = 75 (0i + ^^2) is the charged scalar field, Mr = Mr{fio) and A^ = Ar(/io) 
stand for the renormalized mass and the quartic coupling, respectively. We use 
the modified minimal subtraction (MS) prescription throughout this paper. The 
scale of dimensional renormalization is denoted by /iq. 

In the following, we restrict the discussion to the sector with vanishing to- 
tal charge and, in particular, consider the scattering of two oppositely charged 
particles 4>'^4>~ -^ 4'~^4>~ ■ The pertinent Green function is given by 

X 2^(O|T0(xi)0t(a;2)0t(^^)0(^2)|O)e, (3.2) 

where the subscript c refers to the connected part. 
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The scattering amplitude 

(pi,P2 out|gi,g2 in) = (pi,p2 in|gi, g2 in) 

+ 2(27r)'^5^(pi+p2 -gi -g2)T(pi,p2;qi,q2) (3.3) 

is obtained from the Green function in a standard manner, through amputating 
external legs and multiplying with appropriate renormalization factors. 

Let us now construct the non-relativistic theory, which reproduces this scat- 
tering amplitude at small 3-momenta Pi,q? <^ M^. We first introduce the 
non-relativistic field operators $±(x) ($±(x)) that annihilate (create) a non- 
relativistic (positive/negative) charged particle from the vacuum. The non- 
relativistic Lagrangian Cnr{x) consists of an infinite tower of operators which 
are constructed from the fields $±(x), its conjugated fields <l>j_(x) and space 
derivatives thereof. Time derivatives at order > 2 can be eliminated by using the 
equations of motion and field redefinitions. The guiding principles for construct- 
ing the Lagrangian are the following. 

i) The effective theory must be rotationally invariant and obey P, T discrete 
symmetries, under which the non-relativistic field transforms as follows, 

P<l>±(x°,x)pt = <l>±(x°,-x), r<l>±(x°,x)Tt = <l>±(-x°,x), (3.4) 

where T is an anti-unitary operator. Lorentz invariance is realized in form 
of relations between different low-energy couplings. 

ii) Since the number of the heavy particles is conserved, all vertices in the non- 
relativistic Lagrangian must contain an equal number of ^±{x) and $j_(x). 
Particle creation and annihilation, which occurs in the relativistic theory, 
is implicitly included in the couplings of the non-relativistic Lagrangian. 
In the case considered in this section, these couplings are real and the 
Lagrangian is a hermitian operator. 

iii) The T-invariance implies 

t/:^r(x°, x)rt = /:^r(-x°, x) . (3.5) 

iv) It is necessary to specify counting rules, which allow one to order different 
operators in the Lagrangian. We introduce a formal small parameter w - a 
velocity of the massive particle in units of the speed of light - and count 
each 3-momentum p and space derivatives V as 0{v). Further, the kinetic 
energy of the massive particles p^ — M (or idt — M) is counted at O(f^). 
This rule amounts to an expansion in inverse powers of the heavy mass M. 
All terms in the Lagrangian can be ordered according to the counting in the 
parameter v. This ordering can be performed separately in each ra-particle 
sector - the sectors with different numbers of massive particles do not mix. 
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According to the above rules, the non-relativistic effective Lagrangian in the 
two-particle sector with Q = is given by 



+ 4M^ V*+ ^ $U(<^+<^-) + h.c. I + ^ ($;$+) A (<|.t_$_ 



(3.6) 



where n A f = n(At') + {/\u)v and gfj denote non-relativistic couplings, which 
are functions of the parameters Aj.(/io)) ^r(/^o) and /iq. Ellipses stand for higher 
derivatives. 

The calculation in perturbation theory with the Lagrangian Cnr proceeds as 
follows. One splits the Lagrangian into a free and an interacting part, Cnr = 
^NR + ^NR with C%ji = J2± ^ii^dt - M + A/2M)$±. The free propagator is 
given by 

z(0|T<l>±(x)<l>U0)|0) = / --4t 1 . (3.7) 

^ ' ^^ ^ ±^ ^1 / J (27r)4M+|i-pO_,o ^ ^ 

It vanishes for negative times, which shows that the field ^± (x) indeed annihilates 
the vacuum. As expected, the free particle has a non-relativistic dispersion law 
p^ = M+p^/2M. In order to arrive at a relativistic dispersion law, one has to sum 
up all higher-order relativistic corrections proportional to A^/8M^ ■ ■ ■ . Taking 
into account that A is of order w^, and counting the propagator in momentum 
space as 0(t;~^), the corresponding series for the propagator is of the form a/v"^ + 
b + cv^ + ■ ■ ■ and thus indeed corresponds to an expansion which makes sense. 

For later use, we note that the one-particle states in the non-relativistic 
theory - with the propagator given by Eq. (3.7) - are normalized according to 

(k|q) = (27r)^ 6^{k - q) [non-relativistic] . (3.8) 

Calculation of the scattering amplitude of two oppositely charged particles in the 
non-relativistic theory proceeds analogously to the relativistic framework: one 
evaluates the connected Green functions 

X z^(0|T<l>+(xi)<l>_(x2)$Uyi)'^U2/2)|0)e. (3.9) 

Here, the non-relativistic self-energy corrections in the external legs are summed 
up, so that the correct relativistic dispersion law p° = tf (p) = (M^ + p^)"*^/^ is 
recovered. Because light (massless) particles are absent here, there are no loop 
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corrections to the self-energy, and the pertinent wave function renormahzation 
factor Zi\ffj = 1. Amputating external legs on the mass shell, we finally obtain 

i{27r)'^S'^{pi +P2-qi- q2)TNBiPi, P2; Qi, q2) 
2 2 

= n n^^"^ ,(^(Pi)-P°)n J^^ iw{qj) - q'^)GNRipi,P2;qi,q2)- (3.10) 

3.3 Matching condition 

The relativistic and non-relativistic theories must be physically equivalent at 
small momenta. This equivalence is enforced by the matching condition, which 
enables one to express the non-relativistic couplings Qi in Eq. (3.6) through the 
parameters A^, Mr and the scale fiQ. 

We first formulate the matching condition for the scattering amplitude Tjsfji 
introduced in Eq. (3.10) and start with the observation that the normalization 
of the asymptotic states is different in the relativistic and in the non-relativistic 
theory: in the relativistic case, one has 

(k|q) = (27r)3 2w(p) 5^{k - q) [relativistic] . (3.11) 

From this we conclude that, in order that the cross sections in the relativistic 
and in the non-relativistic theories are the same, the scattering matrix elements 
must satisfy 

2 

T(pi, P2; qi, q2) = n(2^(PO)'^'(2w(qO)'/'T^fl(pi, P2; qi, q2) , 

2=1 

(3.12) 

i.e., the non-relativistic scattering amplitude acquires an extra factor (2tf;(k))~^/^ 
for each external leg. The equality Eq. (3.12), which goes under the name of 
matching condition, is understood to hold at small 3-momenta, order by order in 
the expansion in Pi/M, q^/M. This corresponds to an expansion in the formal 
parameter v introduced in the previous subsection. We will provide illustrations 
below. 

The matching condition for transition amplitudes with more than 4 legs can 
be formulated in a similar fashion. However, in order to do so explicitly, one 
must introduce the concept of connected S'-matrix elements. As we will not need 
these higher order transition amplitudes in the present work, we do not go into 
more details here. 

Finally, let us mention that the matching condition does not always determine 
all couplings g^ separately. This happens, for instance, when one decides to retain 
a non-minimal set of operators in the non-relativistic Lagrangian. The matching 

21 



condition, which is formulated for the scattering amphtudes, then determines only 
certain combinations of the non-relativistic couplings. The point is that exactly 
the same combinations then enter all observable quantities that can be calculated 
from the non-relativistic Lagrangian. Consequently, the non-relativistic effective 
theory leads to a consistent treatment of the low-energy physics also in this case. 
Explicit examples and further details can be found in Ref. [181], chapter IIIc. 

3.4 Tree level 

As an illustration to the matching procedure, consider the scattering process 
Qi + (i2 ^ Pi + P2 in the sector with Q = 0. The relativistic scattering amplitude 
at tree level is Xr/Q. Using the matching condition Eq. (3.12), and expanding 
the non-relativistic amplitude in the CM frame (pi = — p2 = P, qi = — q2 = Qj 

IpI = |q|) in powers of v, we get 

T(pi,P,;qi,q,)_A. 1 _ A. ^^ ,p^ + of-l,V (3.13) 



2w{p)2w{q) 6 2w{p)2w{q) 2AM^ 2AM^^ \M^ 

On the other hand, the non-relativistic amplitude, evaluated from the Lagrangian 
(3.6) is given by 

^ I \ 31 92P'^ ^3(p-q)^ , ^. 4^ fnTA\ 

TNRiPi, P2; ^1' ^2) = ^ - ^^ j^^ + 0(p ) . (3.14) 

From these relations, we can read off the couplings gi, 92,93 at tree level, 

^i = ^ + 0(A^), g, = ^ + 0{X^.), 93 = 0{\l). (3.15) 

The generalization of the tree-level matching procedure to terms containing more 
field operators, or higher derivatives, is straightforward. Here, we only mention 
that the many-particle scattering amplitudes contain singular parts that stem 
from particles traveling near their mass shell between two interactions. These 
singular parts are exactly reproduced at low energies by the tree graphs in the 
non-relativistic theory. The remaining part of the relativistic amplitudes is a 
polynomial in the small external 3-momenta, and can be reproduced by adjusting 
the couplings in Cj^r. 

3.5 Loops 

At higher order in the low-energy expansion, loops also contribute. In the non- 
relativistic theory, they reproduce the non-analytic behavior of the relativistic 
amplitudes at low energy. The polynomial parts can be fixed with the help of 
the matching condition. 

To illustrate the procedure, we again consider the two-particle amplitude. 
The structure of the perturbation series is particularly simple here, see Fig. 1: at 
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Figure 1: Typical diagrams, contributing to the 2-particle elastic scattering ampli- 
tude in the non-relativistic theory. Filled boxes and crosses denote the derivative 
vertices and the self-energy insertions, respectively. 

any order in the loop expansion, only the chain of S'-wave bubbles plus relativistic 
insertions in the internal and external lines plus insertions of derivative 4-particle 
vertices contribute to the scattering amplitude. The elementary building block to 
calculate a diagram with any number of bubbles in the CM frame P^ 
(P*^, 0) is given by 



Vi- 



-V2 



J{P') 



d^^l 
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1 



(27r)^z 
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iJ-PO + lO-^OM+iJ-P-^0 
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(M(P°-2M)) 



1/2 



atP)^4, P°>2M. (3.16) 



The function J(P°) is analytic in the complex P° plane, cut along the positive real 
axis for P*' > 2M. The scattering amplitude is obtained by putting P^ = 2w{p), 
where p denotes the relative 3-momentum in the CM frame. Diagrams containing 
self-energy insertions, or derivative couplings, are evaluated in a similar manner^. 
The perturbative expansion based on the Lagrangian Eq. (3.6) generates a 
systematic low-energy expansion of the scattering amplitudes. The single loop 
J{P^) is proportional to the magnitude of the small relative momentum |p|. 
Self-energy insertions and derivative vertices contribute additional powers of | p | . 
Hence, evaluating loop diagrams with vertices of higher dimension yields contri- 
butions at higher order in the expansion in momenta. We conclude that, at any 
given order in f", only a finite number of diagrams contributes. For example, 
combining two non-derivative vertices through one loop (Fig. Ic) generates a 
contribution at 0{v), and the contribution from the two-loop diagram in Fig. If 
is O(w^), as is the derivative vertex displayed in Fig. lb. 



^Our theory is perturbative in dimensional regularization. We do not consider here the case 
when some of the non-relativistic couplings are "unnaturally" large, which would lead to the 
necessity of partial re-summation of the perturbative series (example: the "pionless" effective 
field theory in the two-nuclcon sector [186, 187]). 
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Figure 2: One-loop scattering amplitude in the relativistic (0V)^ theory. Self- 
energy insertions are not displayed. 

From the above discussion one concludes that in the CM frame, the non- 
relativistic scattering amplitude takes the form 

Tnr = M- (/o + /i^ + /^^ + /3^ + ■ ■ ■) , (3-17) 

where the expansion coefficients fi are polynomials of finite order in the dimen- 
sionless couplings Qi. The matching condition Eq. (3.12) fixes these couplings 
through Aj., Mr and the scale /io order by order in the perturbative expansion. 

As one observes from Eq. (3.17), the analytic structure of the amplitude is 
particularly simple. The only non-analytic piece of the non-relativistic amplitude 
at low energy is given by terms containing odd powers of |p|, whereas the rest 
is a polynomial in the external momenta. As mentioned above, only s-channel 
bubbles generate non-analytic contributions in the vicinity of the elastic threshold 
(graphs shown in Figs. lc,d,e,- ■ ■ ). 

Let us now verify that the relativistic amplitude has the same non-analytic 
structure at low energy as the non-relativistic one. We consider the relativistic 
scattering amplitude at one loop (see Fig. 2), 

r(s, t,u)=To + ^ (2.1(3) + 2.J(t) + J(u)) + 0(\l) , (3.18) 

where Tq is a constant, s = (pi +P2Y, t = (Pi — QiY, u = (pi — 52)^ are the usual 
Mandelstam variables, and where, for x < 0, 



^0^ =T7r^Kxln^— + 2 , a^ = \l . 3.19 

lOTT^ \ ax + 1 J V X 

The amplitude Eq. (3.18) has a more complicated analytic structure than its 
non-relativistic counterpart Eq. (3.17). For instance, its partial waves contain, 
in addition to a cut at s > 4M^, a left-hand cut at s < 0. However, in the 
vicinity of the elastic threshold, the contribution from these distant singularities 
can be approximated by a polynomial. Expanding the functions J(s), J(t), J(u) 
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near s = 4M^, t = 0, m = 0, we obtain 



16vrV(s) = 2 + ^-^ + 0{v' 



\2 



^'"'■'W = -^w^ + o(„' 



16^ J{u) = - (^±^ + 0(«'') . (3,20) 

Substituting this expansion into Eq. (3.18), one immediately sees that the low- 
energy behavior of the non-relativistic amplitude Eq. (3.17) is reproduced. One 
may further check that the non-analytic part oc A^|p| in the expression for 
the amplitude, which is generated by J(s), is automatically reproduced in the 
non-relativistic calculation with the correct coefficient, provided the matching at 
0{Xr) has been performed. 

We end this subsection with a brief discussion of crossing and charge symme- 
try. In the relativistic theory, the scattering amplitude obeys crossing symme- 
try. The crossed channels, along with (f)~^(j)~ -^ 4''^4'~, include also the reactions 
(f)^(f)'^ -^ 0+0+ and 0~0~ —>■ (t>~ (t>~ ■ These reactions are described by a single 
analytic amplitude, derived from the relativistic Lagrangian Eq. (3.1). 

In the non-relativistic case, crossing symmetry is apparently lost because, 
ffist of all, the non-relativistic expansion is performed in the vicinity of s = 4M^, 
t = M = 0. In difference to the relativistic case, an analytic continuation of the 
amplitude from the vicinity of the s-channel threshold to the t- or w-channel 
thresholds cannot be performed, because the distance between these two regions 
exceeds the radius of the convergence of the pertinent Taylor series (see, e.g. Eq. 
(3.20)). Furthermore, our non-relativistic Lagrangian describes only the sector 
with total charge zero, so that the scattering amplitudes for 0+0"'" -^ 0+0+ and 
0~0~ -^ 0~0~ vanish in this theory (within the non-relativistic approach, this 
choice is a consistent procedure). We therefore conclude that crossing symmetry 
is not present in the non-relativistic approach ah initio. At most, one may 
describe the scattering amplitudes in all crossed channels, including all pertinent 
terms in the non-relativistic Lagrangian. The couplings, which are determined 
from matching to the relativistic theory, will then obey the restrictions imposed 
by crossing symmetry. The same conclusion holds for charge invar iance, which 
in general connects sectors of different total charge. 

3.6 Relation to the effective range expansion 

There is one property of the non-relativistic theories that makes the use of this 
framework extremely effective in many areas of hadron physics. The property is 
related to the expansion parameters in the non-relativistic perturbation series. 
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Following the line of reasoning outlined in the introduction, we consider a 
generic hadronic scattering process at momenta which are much smaller than 
any of hadron masses involved and any of the dynamically generated scales. It 
is clear that any consistent EFT at such momenta should operate only with ob- 
servable characteristics (S'-matrix parameters) of this process (since the distances 
involved are already asymptotic). In other words, for e.g. the elastic two-particle 
scattering process, the expansion parameters in the non-relativistic EFT should 
be expressible through the parameters of the effective range expansion (scattering 
length, effective range, shape parameters), rather than directly through the cou- 
pling constant A^, the running mass M^, etc. It can be immediately seen that the 
non-relativistic EFT, which we have constructed here, passes this test: at a given 
order v^ only a finite number of the non-relativistic loops contribute, and the 
coefficients /o, /i, ■ " ' iii Eq. (3.17) can be expressed through a finite number of 
effective-range parameters of the relativistic theory. For example, the coefficient 
of the diagram shown in Fig. la is proportional to the scattering length a in the 
relativistic theory to all orders in Aj. and not only at 0(Ar), since all other dia- 
grams Fig. 1 give vanishing contributions at threshold. Hence, the contributions 
Fig. Ic and Fig. If are proportional to a^|p| and a^p^, respectively, so that the ex- 
pansion in the scattering length a and in the small momenta is correlated. Note 
that, in order to achieve this convenient ordering of the various contributions, 
the regularization scheme - which is used to calculate loops - should not contain 
a mass scale that destroys power-counting. We use dimensional regularization, 
because it has this property. 

4 Non-relativistic effective theories: including 
photons 

4.1 The Lagrangian 

Photons are included in the theory through minimal coupling, 

C = D^(P{D>'(Py - M^(I)(P^ + ^ {(P^^y - J F^^F"" + counterterms , (4.1) 

where </> is a charged scalar field as before, A^ stands for the photon field, D^(j) = 
dfj^cj) + ieAfj_(j), F^i, = d^Ay — d^A^, and the gauge fixing term is not explicitly 
shown. The theory described by the Lagrangian Eq. (4.1) contains two coupling 
constants A,, and e and thus mimics the hadronic atom problem where an interplay 
of electromagnetic and strong effects is present. For convenience, we refer to the 
self-interaction of the scalar field as "strong interactions" . 

The non-relativistic effective Lagrangian that describes this theory at low 
energy is constructed along the same pattern as before. In the presence of pho- 
tons, a few additional considerations should be taken into account. First, one 
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starts from a gauge-invariant non-relativistic Lagrangian. Appropriate building 
blocks are covariant derivatives of the non-relativistic fields, and the electric and 
magnetic fields E and B. Second, the couplings are expressed in units of the 
physical mass M. The infinitely many terms in the Lagrangian are then gener- 
ated by expanding it in inverse powers of the mass. This amounts to an ordering 
according to the number of space derivatives, and according to the number of 
electric and magnetic fields. Third, one may finally use different gauge fixing 
in the non-relativistic and in the relativistic theories. The Coulomb gauge is a 
convenient choice for the non-relativistic theory. 

The lowest-order terms in the non-relativistic theory are 

T e/^i^^-^ + ---V± + lfe(^l^U($+$-) + ---, (4.2) 



6M2 J ^ M^ 

where ^± denotes the non-relativistic field operator for charged particles, Dt^± = 
dt^± =F ieAQ^±, D<l>-t = V$-|- ± ieA^± are the covariant derivatives, E = 
—VAq — a, B = rot A denote the electric and magnetic fields, and hi,gi are 
non-relativistic effective couplings^. For instance, hi is related to the electromag- 
netic charge radius of the particle, see below. The ellipses stand for higher-order 
derivative terms, as well as for non-minimal terms containing E and B. 

We add several comments. First, the set of operators in the Lagrangian 
Eq. (4.2) is not minimal: the term with hi could be eliminated in favor of four- 
particle local terms by using the equations of motion for the Coulomb photon 
Aq. As a result of this manipulation, e.g. the operator proportional to gi would 
receive a contribution. This means that the same linear combination of hi and 
Qi enter the expressions for the two-particle scattering amplitude and for the 
bound-state energies - one does not need to know hi and gi separately. We pre- 
fer to work with this non-minimal set, because it renders the presentation more 
transparent and simplifies the comparison with results available in the literature. 
Second, calculating loops with Cnr generates ultraviolet (UV) and infrared (IR) 
divergences. Therefore, to obtain UV finite Green functions and S'-matrix ele- 
ments at (i = 3, counterterms should be introduced. On the other hand, because 
the final NR amplitudes will be expressed in terms of the UV finite relativistic 
expressions, such a renormalization is not needed: one simply performs all calcu- 
lations at (i 7^ 3 whenever needed, then does the matching at (i > 3 in order to 



^Unless stated otherwise, we use here and in the following the same symbols for the coupling 
constants, for the running and physical masses and for the running scale as in the previous 
section 3, where e = (both in the relativistic and non-relativistic theories). This avoids an 
unnecessary flooding of the text with symbols. Further, to the order considered in the following, 
we do not need to distinguish between the bare and the renormalized charge, so we keep the 
symbol e throughout. See also the discussion in subsection 5.6 on vacuum polarization. 
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avoid infrared divergences that show up in on-shell amphtudes, and sets d = 3 
at the end. This simphfies the calculation considerably, while the final result is 
the same. For a detailed discussion of this point, and for explicit examples, see, 
e.g., Refs. [42,181]. 

The LECs in the the effective Lagrangian Cnr depend on the fine-structure 
constant. Indeed, they can be considered as functions of A^, e, of the mass Mr 
and of the scale /xq, 

hi = hi + ah'i + 0{a^) , fi'i = ^i + ctg'i + 0{a^) , (4.3) 

such that 



dhi dgi 



. (4.4) 



dfiQ dfiQ 

The barred quantities refer to the limit a = 0. 

Further note that, in general, the couplings of the non-relativistic Lagrangian 
are complex, because they contain contributions from processes of particle cre- 
ation and annihilation in the relativistic theory. The imaginary part arises, if the 
threshold for such a process lies below the threshold for two massive particles. 
In the case which is considered here, the imaginary part of e.g. the coupling gi 
will be determined by the 0+0~-annihilation into intermediate states with two 
or more photons - consequently, Im^fi = O(a^) (a generalized unitarity con- 
dition which holds in the presence of complex couplings is briefiy discussed in 
appendix B). For this reason, the general non-relativistic Lagrangian is not a 
hermitian operator, and T-invariance implies (cf. with Eq. (3.5)) 

r£jy^(x°, x)Tt = Ur{-x\ x) . (4.5) 

4.2 Perturbation theory 

To generate the perturbative expansion, one again splits the Lagrangian into 
a free and an interacting part, Cnr = CPnr + '^tvj? with £^^ = —\F^^F^^'' + 
Xl± *^±(^'9j — M + A/2M)$-|-. In addition to the strong bubbles, relativistic in- 
sertions and derivative vertices already discussed in the previous section, diagrams 
may now contain photon lines. The calculations in the Coulomb gauge are done 
as follows. First, the non-dynamical field Aq is removed from the Lagrangian by 
using the EOM. This procedure generates the non-local operator 

A-^ = -(27r)~=^ f d^ke-'^^^'-y^ /k'' (4.6) 

in the Lagrangian. In the following, for ease of understanding, we keep calling 
the pertinent diagrams as "generated by the exchange of Coulomb photons" . The 
propagator of transverse photons in the Coulomb gauge is given by 

/I / k^k^\ 

dxe'''^0\TA\x)A^{0)\0) = -—i5'^--—\. (4.7) 
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Figure 3: One-photon corrections to the non-relativistic scattering amphtude 
(representative set of diagrams) . Dashed and wiggly hnes correspond to Coulomb 
and to transverse photons, respectively. Diagrams containing non-minimal pho- 
ton vertices, relativistic insertions and derivative couplings are not shown. 

In the following, we restrict ourselves to the first nontrivial order in a, and 
consider 2-particle scattering in different channels. In Fig. 3, several examples of 
diagrams with virtual photons are displayed. These are: one-photon exchange in 
elastic channels (Figs. 3a,b), self-energy corrections in the external and internal 
legs (Figs. 3c, d), vertex corrections to the external legs (Figs. 3e,f,g,h), internal 
corrections due to the virtual photon exchange (Figs. 3i,j). At higher orders in 
the momentum expansion, diagrams containing non-minimal photon couplings 
occur, generated e.g. by hi. These are not shown in Fig. 3. 

4.3 Coulomb photons 

Here, we investigate diagrams containing Coulomb photons at order a: one- 
photon exchange (Fig. 3a), vertex correction, where the photon is attached before 
or after all strong vertices (Figs. 3e,f), and an internal exchange between two 
strong vertices (Fig. 3i). 

The contribution from the graph Fig. 3a to the scattering amplitude of oppo- 
sitely charged particles is given by 



TUP^q) 



Aira 



|p-q| 



(4.8) 



This contribution is non-analytic at p, q — i> 0, and thus very different from the 
one encountered in Eq. (3.17). 

Next, we consider the vertex correction Fig. 3e. Let Vc(p, P°) denote the part 
of the diagram which stands on the left or on the right of the first strong interac- 
tion vertex. After integration over the zeroth component of the loop momentum. 
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the integral to be evaluated is 

The contribution to the scattering amplitudes is obtained by evaluating this ex- 
pression at P° = 2w{p). The result (singular at |p| — > 0) is 

TTttM 

V,{p, 2wip)) = -— — - zae, + 0(|p|, d - 3) , (4.10) 

where 

Oc = ^/-^ I j^ - I [r(l) +ln47r] +ln^| (4.11) 

2|p| [rf — 3 2 /i J 

denotes the (infrared-divergent) Coulomb phase^°, and /i denotes the scale of 
dimensional regularization in the non-relativistic theory (not equal to fiQ in gen- 
eral) . 

Finally, we consider the two-loop diagram Fig. 3i. Let Bc{P^) denote the part 
of the amplitude which corresponds to the diagram with one photon exchange 
between two strong interaction vertices. Integrating over the zeroth components 
of the loop momenta, it is given by 

B(P^) = ^[ "^ '- ^^ (412) 

^^^^ ^ (2vr)2^ y po _ 2M - I 111 - l2|^ po _ 2M - I ■ ^ ^ ^ 

Evaluating this expression at P° = 2w{p), we find a result which is again singular 
at IpI — > 0, 

B,{2w{p)) = -^/A(^)+21n^-l-^7rj+0(|pU-3), 
87T [ n J 

A(/x) = /i2('^"=')|^-ln4vr-r'(l)|. (4.13) 

Here, A(/i) contains an UV singularity. For the consistency of the method it is 
important to note that the contributions from diagrams obtained by adding mass 
insertions and/or using vertices with derivative couplings (not shown explicitly 
in Fig. 3) are suppressed by powers of momenta with respect to the leading terms 
Be and Vc- They are not needed in the following. 
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This phase is identical to the one in the relativistic theory [188] at this order. 
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Figure 4: NR propagator in the presence of transverse photons. 

4.4 Transverse photons and the threshold expansion 

The exchange of one transverse photon (Fig. 3b) contributes with 

This term is suppressed by two powers of momenta as compared to the Coulomb 
contribution Eq. (4.8). This is due to the fact that the transverse photon vertex 
in the Lagrangian Eq. (4.2) contains a derivative, whereas the corresponding 
vertex for the Coulomb photon is point-like. In other words, the ordering in 
the Lagrangian is preserved in the contribution to the amplitude in this case: a 
term which is suppressed in the Lagrangian by powers of the mass M generates a 
term in the amplitude which is suppressed by powers of momenta. This feature 
is called power counting in the following. 

For the non-relativistic theory to make sense, power-counting should also 
persist in loops. This is not, however, the case if one uses standard dimensional 
regularization in the loop calculations. To demonstrate this fact, we consider in 
detail the virtual photon contribution to the self-energy of the massive particle at 
one loop. Summing up one-loop diagrams with a transverse photon to all orders 
(see Fig. 4), one obtains 

where f2 = M + ^— p° = O(f^) and S^(r2, p) denotes the one-loop diagram 

[Coulomb photons do not contribute to the self energy, because the Coulomb 
photon propagator does not depend on the zeroth component /°. Therefore, 
one may close the contour of integration over l^ in that half-plane where the 
propagator of the massive particle has no singularity.] 
The self-energy is obviously of order p^ , 

E,(fi,p) = ^S,(f],p). (4.17) 
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We combine the denominators by use of the Feynman-Schwinger trick and obtain 
at p = 

1 VL f 1 \ °° f VL\'^ 

S.(aO) = -^ln— (l --==)+ 5: a„(-j .a„^0, (4,18) 

The result consists of two parts: a non-analytic piece, which starts with a term 
proportional to filnl^, and a polynomial part, which starts with a constant. The 
non-analytic part contributes all in all at order v^ to the self-energy and therefore 
respects power counting: its contribution is suppressed by v^ with respect to the 
tree contribution Vt in the denominator of the last term in Eq. 4.15. On the other 
hand, the polynomial part breaks power counting: it generates a term which is 
of the same order as the tree level one. Considering derivative vertices in Fig. 4, 
it is seen that the corresponding diagrams contribute at the same order in w, so 
that the low-energy effective theory is inconsistent in this sense. 

It is easy to identify the reason for this failure. First, we note that the 
integration momenta /°,1 can vary in wildly different regions. For example, one 
or both of them can be of order of the hard scale M, be "soft" (~ v) or "ultrasoft" 
(~ f^). Consequently, this integral is given by a sum of different terms, where 
each one corresponds to an integration in a different regime (see, e.g., Ref. [178]. 
For a detailed investigation of this statement for a specific one-loop diagram in 
the non-relativistic effective theory, see also Ref. [45]). Therefore, in order to 
construct a consistent low-energy theory, one should modify the Feynman rules 
to get rid of the high-energy contribution. In a most general way, this can be done 
by using the so-called "threshold expansion" [42,175,176,178-182]^^. A simple 
algorithm for the calculation of Feynman diagrams that occur in the present work 
can be formulated as follows: 

i) Perform contour integrals over the zeroth components of the momenta. 

ii) Assume all 3-momenta (external as well as the integration momenta) to be 
much smaller than the hard scale M; expand the integrands in the small 
3-momenta. 

iii) Interchange the order of integration and expansion, then perform integra- 
tions in dimensional regularization. 

After applying the threshold expansion, the integrand becomes a homogeneous 
function in small kinematic variables, and the counting rules are restored. On 
the other hand, the high-energy contribution, which is polynomial at low energy, 
and which upsets power counting, is removed from the result. 



^^In case of TriV scattering at one loop in ChPT, the threshold expansion is obviously equiv- 
alent to the infrared regularization introduced in Ref. [146]. 
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Let us see how this procedure works in the present case. The contour inte- 
gration in Eq. (4.16) gives 

It is seen that 1 should scale like Q ^ v^, otherwise the expansion of the inte- 
grand leads to no-scale integrals that vanish in dimensional regularization. With 
such scaling assumed, the individual terms in the denominator of the last term in 
Eq. (4.16) scale like v'^,v'^,v^,v^. Therefore, this denominator is not a homoge- 
neous function in v, leading to a breakdown of the counting rules. Applying the 
threshold expansion, we obtain (we put a hat on threshold-expanded quantities) 

T (O ^ e' f d'^l f , {pir\ 1 / 1 



. (^-41..^. -V (4.20) 



M 2MJ (1^+|1|)2 
Performing the remaining integration, we obtain 



o2 



67r2M2 



p^n |L(/i)+ln— -U + 0(M"3,rf-3) ,(4.21) 



Hfi) = /-^(^-^(r'(l) + ln4vr + l)). (4.22) 

As expected, the non-analytic piece ~ filnfi stays unaffected. On the other 
hand, the part of the polynomial piece which scales like ~ v"^ has disappeared 
from the result. Since the threshold-expanded self-energy is proportional to ^2, 
the position of the particle pole is not affected by radiative corrections. It is 
furthermore seen that in the vicinity of the mass shell, the propagator Eq. (4.15) 
is unaffected, 

D{p) -> 1, n^O (rf>3), (4.23) 

as a result of which the external lines are not renormalized (wave function renor- 
malization constant Z = 1.). Further, it is seen from Eq. (4.21) that the intro- 
duction of a counterterm rendering the two-point function finite at d = 3 would 
add a contribution of order p^ to Z, which diverges as (i — > 3. 

We investigate a second example, and consider a diagram that describes the 
exchange of a transverse photon between the initial pair of charged particles 
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with momenta p^ = (p°,±p), see Fig. 3g. The contribution of this diagram is 
proportional to the integral 









X — fT^ . (4.24) 

We put the external particles on the mass shell, p^ = M + p^/(2M) +0(p^), and 
perform the threshold expansion in the integral. Note that with this procedure, 
the integrands should be expanded in the O(p^) remainder of p°. The threshold- 
expanded integral can be rewritten in the following manner, 

J n l^ e^ f d^l I [ ^ (pl)^ 



M J {2'kY P V 12 y P - 2pl 

e^ Ipl ie^ \v\ f , ^ , 2|p|\ 



16M SvrM V /i 

The divergence at c? = 3 is an infrared one. We see that power-counting is at work 
also here: in comparison to the exchange of a Coulomb photon, this contribution 
is suppressed by p^. This particular contribution to the scattering amplitude 
vanishes at threshold as well. 

4.5 Matching 

The matching condition given in Eq. (3.12) is universal and holds in the presence 
of photons as well (as photons are relativistic particles, the corresponding states 
have the same normalization as in the relativistic theory). However, as we have 
explicitly seen, the singularity structure of the amplitudes near threshold is dif- 
ferent from one given in Eq. (3.17). Further, in order to determine all couplings in 
the non-relativistic Lagrangian Eq. (4.2) separately (e.g. the coupling constant 
/ii), it is convenient to consider amplitudes with external photon legs as well. 

We start with the coupling hi which is related to the charge radius of the scalar 
particle. We closely follow the method described in Ref. [173]. We consider the 
transition amplitude S/j of the charged particle in an external field A'f* and 
define the relativistic form factor F{t) through the linear term in the expansion 
in A-*, 

S/i = 2w(p) {2^f5^{^ - q) + ieA''^\p -q){p + qYF{t) + 0{e^) [relativistic] , 

(4.26) 

where A^^^{p — q) denotes the Fourier-transform of the external field, and where 
t = {p — qY- At a small t, 

F{t) = l + l:{r')t + 0{t'), (4.27) 

o 
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where (r^) stands for the mean square radius in the hmit a = 0. 

The definition of the form factor in the non-relativistic case is very similar, 

S/, = i27iYS'\p - q) + «ei;f*(p - q)Fj^j^{p, q) + 0{e^) [non-relativistic] . 

(4.28) 

The matching condition Eq. (3.12) yields in our case 

{p + qrFit) = i2wip)2w{q)Y/'Fi:,j,{p, q) . (4.29) 

The zeroth component of the non-relativistic form factor at the lowest order in a 
can be directly read off from the Lagrangian Eq. (4.2), considering the coefficient 
of the term linear in the field A^, 

and Eq. (4.29) finally gives 

hi = M^{r^)+0{a). (4.31) 

This relation is an example of matching to threshold parameters, discussed in 
subsection 3.6. Namely, the quantity (r^) is the charge radius of the relativistic 
scalar particle at a = and to all orders in the strong coupling constant A^. In 
order to perform matching in terms of A,., all what one has to do is to calculate 
(r^) in the relativistic theory at a given order in A^ and substitute this result into 
Eq. (4.31). Additional calculations that would invoke non-relativistic EFT are 
not needed. 

We now perform the matching of the coupling gi. For this, we consider the 
elastic scattering amplitude of two oppositely charged particles. Let us restrict 
to diagrams of order a. At this order, the relevant diagrams can be divided in 
two groups: those which can be made disconnected by cutting one photon line, 
and those which cannot, 

T = Ti^ + f, TNR = T],l + f^R. (4.32) 

For the relativistic theory, this splitting is unambiguous. On the other hand, 
the non-relativistic Lagrangian is not unique - e.g., depending on whether one 
keeps the coupling hi in the Lagrangian or removes it using the equations of 
motion, the corresponding term in Tj^j^ will be absent. Here, we stick to Cnr 
in Eq. (4.2), which generates the relativistic one-photon exchange amplitude to 
the relevant order in the momentum expansion. The matching condition then 
holds separately for the 1-photon exchange contribution, and for the remaining 
part of the amplitude. For this reason, we retain only the one-particle irreducible 
amplitudes T and T^ji in the matching condition. 
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We now first discuss the singularity structure of tlie non-relativistic scattering 
amplitude at threshold, and combine information about the strong sector from 
Eq. (3.17) with the known threshold behavior of the various diagrams with virtual 
photons, given in Eqs. (4.10), (4.13), (4.21), (4.23) and (4.25). It can be easily 
seen that at order a 

i) diagrams with transverse photons do not contribute at threshold; 

ii) each strong loop gives rise to a suppression factor |p|, so that only diagrams 
with not more than one strong bubble contribute at threshold; 

iii) in the expression for T/vr, one may combine the strong diagram shown in 
Fig. la together with the Coulomb correction to this diagram in Fig. 3e 
(there is second diagram of this type, where the Coulomb photon is ex- 
changed after the strong vertex). One may further check that at order a 
the (infrared-divergent) Coulomb phase 1 + 2ia9c = e^*"^'' + 0{a^) can be 
factorized in the whole non-relativistic amplitude. 

Using Eqs. (4.10), (4.13), (4.23), (4.25), we find that the threshold behavior 
Eq. (3.17) is modified - at order a - by virtual photon contributions in the 
following manner, 

e-'^^'^T^n = ^+A2\n^ + A, + 0{\p\). (4.33) 

IpI ^^^ 

A straightforward calculation of the coefficients gives 

Tragi . _ agf 
^1 — ^zirr 5 ^2 
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Suppose now that one calculates the relativistic amplitude at order a, and to any 
order in A,.. In order to be consistent with the non-analytic behavior predicted 
by the non-relativistic theory, the threshold behavior of the relativistic 1-particle 
irreducible amplitude must be given by 

e-2-^^f = ;^ + ^2 In %1 + T + 0(|p|) . (4.35) 

|P| J-^ 

The infrared-finite quantity T will be referred to hereafter as the "relativistic 
threshold amplitude." It can be considered the 0{a) generalization of the stan- 
dard definition of the scattering amplitude, valid in the presence of real and 
virtual photons. [ Remark: While defining the threshold amplitude, one has to 
discard in particular terms which diverge logarithmically at threshold. This pro- 
cedure is ambiguous, and depends on the choice of the scale at which the log- 
arithm is set to vanish. In Eqs. (4.33) and in (4.35), it was chosen to be the 
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reduced mass of the system. We are free, however, to choose instead a different 
scale, as a result of which the threshold amplitude will change accordingly. On 
the other hand, observable quantities, which are expressed through the threshold 
amplitude (e.g. bound-state energies of hadronic atoms), are independent of this 
scale. In the present model, this can be verified order by order in the perturbative 
expansion.] 

Finally, we express the coupling gi in terms of T, 

r = Ag^h - ^ ('a(^) +ln^ - 1 - 2n^l\ \ . (4.36) 

This completes the matching of gi and of hi. 

5 Bound states 

5.1 Introductory remarks 

In the previous sections, we constructed a systematic non-relativistic field theory, 
which is equivalent to the underlying relativistic theory at small 3-momenta. The 
equivalence is achieved by performing a matching procedure order by order in 
the momentum expansion. The non-relativistic approach does not provide new 
information about the low-energy behavior of the scattering amplitudes, because 
this low-energy behavior is the information that enters the matching condition. 
The non-relativistic approach becomes useful in the description of the shallow 
bound states in the theory. To be specific, we consider bound states in the system 
described by the Lagrangian Eq. (4.1). In this model, in the vicinity of the elastic 
threshold, there exists a tower of nearly Coulombic bound states, whose energies 
are approximately given by 

E„ = 2M-^-^, n=l,2---. (5.1) 

Due to the combined effect of the strong and the residual electromagnetic inter- 
action, the energy levels are displaced and acquire a finite width. Our aim is to 
find the pertinent corrections to the leading order formula Eq. (5.1). As already 
mentioned, the non-relativistic Lagrangian constructed in the previous section 
offers a simple, elegant and very efficient framework to solve the problem [37]. 
We illustrate in this section the procedure. 

5.2 Coulomb problem 

We start from the unperturbed solution which corresponds to a pure Coulomb 
potential. The unperturbed situation is described by the Lagrangian which is 
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obtained from the original non-relativistic Lagrangian Eq. (4.2) by discarding ev- 
erything but the minimal coupling of the Coulomb photons to the non-relativistic 
massive particle, 



Eliminating the field A^ through the EOM, the Hamiltonian becomes 



(5.2) 



2 

($:^$+ - $L$-)A-^(<l>i^<l>+ - $L$-) } . (5.3) 



^ ' "st dS _ rfit dS ^A-l('<TitdS _ rfit 



2 

We introduce creation and annihilation operators, 

$±(0,x) =y"^e^P''a±(p), [a±(p),al(k)] = (27r)353(p - k) , (5.4) 

and construct the bound state of two charged particles in Fock space, 

|^nim,P)= / T^^^„«„(k)|P,k), 

|P,k) = aU^P + k)aU^P-k)|0), (5.5) 

where \l/„;m(k) denotes the Schrodinger wave function in momentum space, and 
n,l,m stand for the principal quantum number, the angular momentum and its 
projection in the z— direction, respectively. See appendix A for further notation, 
in particular, for an explicit expression of the wave functions \l/„;m(k). The state 
vectors (5.5) satisfy 

(Ho + Hc)|^„,,n,P)= r^n + ^VnZ^n,?). (5.6) 

To proceed, we introduce the resolvent for the Coulomb problem, 

Gc{z) = J^ „ = Go(^) + Go(^)HcGc(^) , (5.7) 

2 - Ho - He 

whose matrix element between the states |P,k) develops poles at z = £"„. To 
remove the CM momentum of the matrix elements, we introduce the notation 

(q|r(^)|p) = yi^ (P,q|R(^)|0,p) , (5.8) 
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where R(-z) is any operator in Fock space, and where r{z) denotes the pertinent 
operator in the CM system. The matrix element of the resolvent G^ is related 
to Schwinger's Green function [189], 

( \ t ^\ ^ (27rmq-p) 1 Ana 1 

(q|gc(2;)|p) = — Z2 






47iar]I{E- q,p) -^ , (5.9) 



^- M ^ M 



with 



f X ^ dx 

I{E;q,p)= -, —, (5.10) 

Jo [(q-p)2x + W«2(i_^)2(^_2!)(i^_E!)]' 

where r] = ^a {—E/M)~^^'^ and E = z — 2M. This function develops poles at 
rj = 1, 2, ■ ■ • . The Coulomb wave functions in momentum space can be read from 
the residues of these poles [189]. 

5.3 Feshbach formalism and the Rayleigh— Schrodinger 
perturbation theory 

We now construct a systematic perturbation theory to include the remaining 
strong and electromagnetic interactions which are contained in the Lagrangian 
Eq. (4.2). It is convenient to apply the so-called Feshbach formalism [38,39]. 
The Hamiltonian of the system takes the form 

H = Ho + Hc + V, (5.11) 

where V stands for all interactions other that the static Coulomb potential. Con- 
sider the full resolvent 

G{z) = j^. (5.12) 

It satisfies the equation 

G = Gc + GcrGc, r = V + VGct . (5.13) 

The spectral representation of this resolvent contains a sum over all states. In 
order to calculate the energy shift of the n, I level, it is convenient to remove the 
corresponding unperturbed pole contribution {z — En)~^ (in the CM frame) from 
Gc{z) by defining 

G^' = Gc I 1 -J2 J -^s\'^nlm,P){^nlm,P\ J • (5.14) 
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We further introduce 

r"' = V + VGj'r"' , (5.15) 

and find for G the representation 

G = G^' + G^'r"'G^' + (1 + G^'T"')n„i(l + r"'G^') , (5.16) 



where 



t^^ (2-)^--S-^n-(v&„z|r"'(.;P)|vI...)' ^ ^ 

and 

(vl/,;|T"'(z;P)|vl/„0 = y — -^(vI>„,„,P'|t"'(^)|v1/„;^,P). (5.18) 

Here we used the fact that the right-hand side does not depend on the magnetic 
quantum number m - here and in the following, we therefore omit the subscript 
m, whenever no ambiguity can occur. The singularity at z = En is absent in 
the barred quantities. Therefore, the pertinent pole must occur through a zero 
in the denominator of the expression Eq. (5.17). In the CM frame, the relevant 
eigenvalue equation to be solved is 

Znl - En - {^nl\T''\z)\^nl) = , (5.19) 

where the matrix element denotes the quantity on the left-hand side of Eq. (5.18), 
evaluated at P = 0. To ease notation, we often will omit the indexes in z^i- 

The master equation (5.19) is a compact form of the conventional Rayleigh- 
Schrodinger perturbation theory, valid for unstable systems as well. It fixes the 
convergence domain of perturbation theory: the theory is applicable as long as the 
level shift does not become comparable to the distance between the ground-state 
and the first radial-excited Coulomb poles. Equation (5.19) is valid for a general 
potential - containing e.g. the interaction with the transverse photons - since in 
the derivation, we did not use the explicit form of the interaction Hamiltonian. 

Finally, the level shifts and the widths can be obtained from the solution of 
the master equation, 

Z = En + AEnl , (5.20) 

where AEni is complex. As we shall see later, the poles of the resolvent G{z) 
occur on the second Riemann sheet in the complex 2;-plane, in accordance with 
analyticity requirements. 
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5.4 Energy levels 

The energy shifts admit a Taylor series expansion (up to logarithms), 

AEni = ttnia^ + hnia^ In a + Cnia^ + O(a^) . (5.21) 

We now show how to calculate the coefficients a^/, hni and c„;. 

We expand the matrix element in the master equation (5.19) in a Taylor series 
in [z — En) and obtain 

(*„z|t"'(EJ|^„;) 

z-En = ^F^ , ' "' + • • • . 5.22 

Because the numerator is a quantity of order a^, the second term in the denom- 
inator starts to contribute at order a^ and can thus be dropped [42], and one 
finds 

z-En = (^„z|r'^'(^n)|^nO + 0(«^) . (5.23) 

It remains to evaluate the pertinent matrix elements of the operator T"'(i?„). 

We start with the determination of the perturbation V from the Lagrangian 
Eq. (4.2). One ffist eliminates the non-dynamical field A^ by using EOM. At 
the order of accuracy required by Eq. (5.21), it suffices to retain only following 
terms: 



9i 

M2 



Hs = -^ /rf3x$t^$L$+$ 



lb ^ 

2 7 /• 

e^Hf,n = -^y rf=^x($V$+-$^_$-)(<$+-$^_$-). (5.24) 

Additional terms generated by the Lagrangian Eq. (4.2) do not contribute to 
the energy shift at O(a^), and are hence omitted. Evaluating the energy shift is 
now straightforward. Solving Eq. (5.15) by iterations, at the order of accuracy 
required we get 

t^'{z) = H5 + H5G^'(z)H5 + H« + e^Uf^n + e^Il^G^'{z)Il^ + ■■■ . (5.25) 
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Figure 5: Diagrams contributing to the electromagnetic energy shift: a) relativis- 
tic insertions; b) finite-size effect; c) self-energy corrections due to transverse 
photons; d) transverse photon exchange. 

Then, using Eq. (5.23), we find 

+ e2H^Gj'(E„)H^|^„0 + O(a^). (5.26) 

In order to evaluate this quantity, one inserts a complete set of states between the 
various operators - only a finite number of terms is non-zero in each sum - and 
calculates the integrals over the intermediate momenta. It is convenient to further 
split the energy shift Eq. (5.26) into so-called "strong" and "electromagnetic" 
parts, 

Tistr -p I I A ZT'cm 



AEni = 6io iAEf - 2 Tn 1 + ^KT + 0{a') [scalar QED] . (5.27) 

The (real) electromagnetic shifts AE^ are evaluated below. The above equation 
should be taken as a definition of the bracketed term on the right-hand side. 
Since the strong Hamiltonian is local, the strong corrections to the energy levels 
with angular momentum / ^ are suppressed by additional powers of a. As a 
result of this, they vanish at the accuracy considered here, as is indicated by the 
Kronecker symbol in Eq. (5.27). 

The naming scheme for the various corrections should not be understood lit- 
erally. For example, the "electromagnetic" contribution contains the coupling hi 
which, according to Eq. (4.31), depends on the quantity (r^) which is a func- 
tion of the "strong" coupling constant A^. Vice versa, there are electromagnetic 
corrections to the coupling gi which enters AE^'^, see Eq. (4.3). 

Let us first consider the electromagnetic shift. At O(a^), it has contributions 
from the relativistic insertions, finite-size effects, and transverse photons (see 

42 



Fig. 5), 



/ ^^P l,Tr ^„M2 P^ , e^^l,:^, ,„,|2 






/(l;^(l^^"'(p)^'^^P'^^^"'^^^- ('-2^ 



The contribution from the self-energy diagram Fig. 5c is of order a^ and is 
therefore not displayed. The diagram with one transverse photon exchange (Fig. 
5d) gives 



4M^|p - q| V- ■ - |p-qp;^ + jj + |l + |p_q| 

= / ((p + qy^-lP-iqyU.... (5.29) 

4:M^\p-q\^\ |p-q| / 

The evaluation of the pertinent integrals in Eq. (5.28) for arbitrary n, I is made 
easier by eliminating the terms of the type p^^„;(p) by use of the Schrodinger 
equation in momentum space [44]. We obtain 

.™ „4.w «m ,11 1 \ , a'Af'{0 



A£;r = «-M^^ + — -^-5^^^^j+^-^^*,„. (5.30) 

Note that in the calculations of the electromagnetic shift, we have replaced Gc{z) 
by Go (2). It can be verified by direct calculations, that this does not affect the 
result at O(a^). In the limit of point particles, the above result coincides with 
the one from Ref. [190]. Next, we turn to the calculation of the strong shift. 
According to Eq. (5.26), one has 

^ET -\t^ = -^(1 - |,(gc°(i^.))) + 0{a^) , (5.31) 

where 

^^^^^"^^ =J J2^^j2^^ (P|g£°(^n) |q) = ^^ (^A(/i) + 1^ ^ - 1 + ^n(«) j , 

s„(a) =2(^(n)-V^(l)-- + lna-lnn), ip(x) = T'(x)/T(x) . (5.32) 

n 

Finally, we substitute the value of the coupling constant gi from Eq. (4.36) into 
the expression for the strong energy shift Eq. (5.31). In this manner, one obtains 
an expression which contains only quantities defined in the relativistic theory, 
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In order to calculate the strong shift of a given level, we take the real part of 
the this expression and obtain 

Equations (5.30), (5.33) and (5.34) are the main results of this section. They 
provide, at this order, a complete expression for the total (electromagnetic and 
strong) shift of the energy level with a given quantum numbers n, I. Furthermore, 
Eq. (5.34) is the generalization of the DGBT formula for the strong shift of the 
energy level of the bound state [8, 191-193] to next-to-leading order in a. 

5.5 Decay into 2 photons 

In the present theory, the ground state can only decay in two or more photons. 
This is a hard process, with a mass gap of order of the heavy particle mass M. 
Consequently, within the non-relativistic theory, the decay is described through 
the imaginary parts of the coupling constants. At leading order, this is the 
coupling gi in Eq. (5.31), with Iragi = O(a^). According to Eq. (5.22), we may 
then use Eq. (5.31) to calculate this width up to and including terms of order a^. 
The result is 

Tl' = ^^lmgi + 0{a^). (5.35) 

An 

At the order of accuracy we are working, the imaginary part of gi is given by (see 
also appendix B) 



Im^i = -^ J2^2nY6\p,+p2 - q, - g;)|T<^+(Pi)^-feH2,|2 



16 
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threshold 



(5.36) 

Evaluating the matrix element at tree level, we find 

lmgi = 7ra'^ + 0{a^Xr,a^), (5.37) 

see also subsection 8.4. We finally arrive at the following result for the two-photon 
decay width of the ground state [54] , 

r?^ = ^ + 0(a^A„«^). (5.38) 

In the real world, additional channels are open for decay (strong as well as elec- 
tromagnetic). These will be treated in later sections. 
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Figure 6: Modification of tlie static Coulomb potential by an electron loop. 

5.6 Vacuum polarization due to electrons 

In this subsection we separately consider a particular correction to the hadronic 
atom observables, which does not emerge in the model described by the La- 
grangian Eq. (4.1), but is present in the case of real hadronic atoms - namely, 
the correction due to vacuum polarization induced by an electron loop. In or- 
der to take this effect into account, one would have to use the formulation of 
the low-energy effective theory of the Standard Model not in terms of photons 
and hadrons only [194-197], but to consider the explicit inclusion of the leptonic 
sector of the Standard Model as well [198] . Further, although in the treatment 
of the hadronic atoms weak interactions can be safely neglected, one has still 
to take into account the presence of electrons, which couple to photons. The 
reason for this is that the electron mass rrie ^ 0.5 MeV is numerically of order 
or smaller than aM - a scale which is "resolved" in the non-relativistic theory 
(M stands now for a typical reduced mass in the hadronic bound state). Exam- 
ining, however, the possible contributions from the explicit electron degrees of 
freedom to the observables of hadronic atoms, one may easily check that these 
all vanish at the next-to-leading order in a with one exception. Modification 
of the static Coulomb potential by an electron loop (vacuum polarization) gives 
rise to corrections in the hadronic atom observables which are amplified by large 
numerical factors - powers of M/rrie, which emerge from the calculation of the 
matrix elements in Eq. (5.19). If one would count the quantity aM/mf, at 0(1) 
in order to get rid of these large factors, the vacuum polarization effect starts 
to contribute to the bound-state energy at 0{a^), i.e. at the same order as the 
leading-order strong shift which is calculated by using the DGBT formula. To 
the decay width, the contribution of the vacuum polarization effect leads to a 
0{a) correction^^ . 

Since the effect from the presence of electrons appears only at one place, it 
would be counterproductive to carry out a full-fledged inclusion of the explicit 
electron degrees of freedom into the framework. More easily, the same goal may 
be achieved, using e.g. the methods described in Ref. [200] (see also Refs. [41, 201- 
203]). Here one first integrates out the electron field in the generating functional. 



""^^For a systematic discussion of vacuum polarization contributions within a potential model 
framework, see, e.g., Ref. [199]. 
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At the order of accuracy we are working, this amounts to a modification of 
the kinetic Lagrangian —^F^^F^'^ by the electron loop. Furthermore, at this 
order one may safely neglect the corrections to the transverse photon propagator, 
performing the static limit k"^ — »■ — k^ in the photon self-energy part. Finally, 
the whole effect reduces to the well-known modification of the static Coulomb 
potential in Eq. (5.3) through the electron loop (see Fig. 6). In momentum space, 
this modification reads 



Ana Ana Aa^ f^ ds 1 / 2ml\ / Ami 

= -^ + V^"'(k). (5.39) 

[Here, one has removed the UV-divergence in the propagator by charge renormal- 
ization.] Considering the second term in Eq. (5.39) as a perturbation, one may 
calculate the corrections to the bound-state observables in a standard manner. 
Namely, adding V™'^ to the Hamiltonian in Eq. (5.24) and iterating, we arrive at 
the following correction terms to the right-hand side of Eq. (5.26), 

(^„;|V™^ + V^-=GJ'(E„)H5 + H5G^'(E„)V™'= + V™^GJ'(E„)V™= + ■ ■ ■ |*„0 . 

(5.40) 

Starting with the first term, we note that the relevant matrix elements 

AE^r = (^„,|V^^=|^„0 (5.41) 

have been evaluated analytically for any n, / in Ref. [200]. For completeness, here 
we reproduce the integral representation provided in Ref. [200, Eq. (B.3)] which 
is valid in the generic case of a bound state of two oppositely charged particles 
with arbitrary masses, 

Ap.vac ^ /icQ^ "y^ f n-l-l\ f n + l \ .2in^l~k^l) 

37m2 A^V ^ J \n-l-k-l J"^'^ 

/•I 2Z+2fc+l „™ 

X dx-^^—^VT^l^i2 + x'); in = ^^. (5.42) 

Here fic denotes the reduced mass of the bound system, and me is the electron 
mass. 

The second term in Eq. (5.40) contains both V™'^ and H5. Because the strong 
Hamiltonian Hs is local, the corrections in the states with / 7^ start at higher 
order in a than in the S'-states, see also the remark after Eq. (5.27). We therefore 
stick to the corrections in S'-states. It is convenient to introduce the quantity 

5*no(p) = / (0? (0J (Plgc (^n) |k) V^-(k - q)vl/„o(q) , (5.43) 

46 



which is called "modification of the Coulomb wave function at the origin, due to 
vacuum polarization." The corrections are then most easily expressed in terms 
of the ratio 

Sr = 25§„,o(0)/§„o(0) . (5.44) 

The next-to-leading correction to the energy shift and width emerges through the 
modification of the value of the wave function at the origin (cf. with Eq. (5.31)), 

|^no(0)r^|^no(0)|'(l + 5r). (5.45) 

The quantity 6^^^ is explicitly evaluated for n = 1 in Ref. [200]. At the present 
experimental accuracy, the correction 6^^^ is not negligible for ttH , where it 
amounts to approximately half a percent. 

The notion of a "modified Coulomb wave function" appears here for the sake 
of convenience only - in this manner, one may easily parameterize the correction 
term and carry out a comparison with other approaches. The Feshbach formalism 
does not refer to the exact wave function of the bound system at all - this might 
even fail to be a well-defined quantity in the case of meta-stable bound states. 
All calculations in the Feshbach formalism are performed in terms of the resolvent 
G{z), which is a well defined quantity. 

We expect that the remaining terms in Eq. (5.40) are very small, and we 
discard them in the following. 

5.7 Energy shift and width: summary 

Amazingly enough, the above results carry over nearly unchanged to the real 
world of hadronic atoms, described in the framework of ChPT. For this reason, 
we collect the results of this and the previous section for later reference. 
The relevant master equation to be solved is 

^-E„-(*„z|r"'(^)|^„0 = O, (5.46) 

where the matrix element on the right-hand side is the one in Eq. (5.18), evalu- 
ated at P = 0. We also introduced several energy shifts, 

z = En + AEni , (5.47) 

AEr^i = AEJ + AE:T + Sio(^AEf-'-Tr}j+0{a'). (5.48) 

The understanding of these relations is as follows, i) AEni is a complex quantity, 
defined through Eq. (5.47). ii) The real shifts AE^ and AE^ are defined in 
Eqs. (5.30) and (5.41), respectively, iii) The remainder AE^i — AE^ — AE^"^ is 
split into the real components AE^'^ and r„, which are related to the threshold 
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amplitude in the underlying relativistic theory. Including vacuum polarization 
effects, 

where T is the threshold amplitude defined in Eq. (4.35). 



6 On DGBT formulae in ChPT 

In this section, we investigate in some detail the procedure to extract informa- 
tion on the hadronic scattering lengths from the energy spectrum of hadronic 
atoms. To start with, we perform a Gedankenexperiment: let us assume that the 
fundamental relativistic Lagrangian is given by Eq. (4.1) (the effects of vacuum 
polarization do not play any role here and are ignored). Assume further that one 
measures the energy levels of the two-body bound state and extracts the real part 
of the threshold amplitude Re T via the expressions Eqs. (5.27), (5.30) and (5.34). 
Can we purify Re T from electromagnetic interactions and determine the "purely 
strong" threshold amplitude? The question is tailored to quantum-mechanical 
models, where one considers bound states formed by the sum of Coulomb and 
short-range strong potentials, and where the Coulomb potential may be easily 
switched on or off. After this investigation, we discuss the analogous question 
in the framework of ChPT. We follow the ideas outlined in Ref. [204], see also 
Refs. [205-210]. 

6.1 Electromagnetic corrections in scalar QED 

The amplitude T defined in Eq. (4.35) depends on two coupling constants, A^ 
and e. We expand T in powers of the fine-structure constant, 

ReT = T + aT^ + 0{a^). (6.1) 

Here, T denotes the "purely strong" threshold amplitude. Consider now the 
energy shift of the bound state, determined by Eq. (5.34). At leading order we 
have the standard DGBT formula [8], 

which may be used to measure T, 

T = £^ + 0{a) , £:„ = — . 6.3 
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To include higher order corrections, we rewrite Eq. (5.34) as 

Sr = f(l-^|^T)+aTi + 0(a^). (6.4) 

This formula leads - provided that one can calculate T^ - to a refined determi- 
nation of the threshold amplitude, 



T = £T (l + ^^ C J - "^1 + 0{a^) . (6.5) 

It turns out, however, that T and T^ cannot be uniquely defined. Indeed, let us 
imagine a perturbative expansion of the amplitude, and let us retain the terms 
up to and including one loop. The result will have the structure 



ai 2 M2 a2 M2 

— A,, In -^ + — Ka In -^ 

2 /io 2 /i^ 



Re T = Ar + — Aj, In — y + ir ^rCt In — ^ + ^i-^r + ^2Artt + • ■ ■ , (6.6) 



where Oj, 6, are pure numbers, and where the ellipsis stands for the contributions 
beyond one loop. This amplitude is scale independent, because A^ satisfies the 
RG equation 

jjiQ- — \r = aiX"^ + a2\rO: + ■ ■ ■ . (6.7) 

d/io 

What is the meaning of the splitting in Eq. (6.1)? The amplitude T is given by 



T = A^ + ^A^ln^ + feiA^ +... , (6.8) 






where A^ and M are the couplings and masses in pure theory, compare the 
footnote at the beginning of section 4. The coupling constant A^ satisfies the RG 
equation with a = 0, 

fXo^Xr = aiXl + 0{\l)- (6.9) 

d/io 

The second term in the representation Eq. (6.1) is then simply the difference 
T — T. In order to actually calculate it, we need to know the relation between 
Xr and Xr [204]. For this, we may specify the scale /ii at which the two couplings 
coincide. We denote the corresponding coupling by Xr{fio;fii), 

KifJ'O] /ii) = Xrifio) ( 1 - aaa In — H ) , (6.10) 

where the arbitrary scale Hi is referred to as the matching scale. As A^ is inde- 
pendent of /ii, we have 

1^1-: — Ar(/io;/^i) = "OsAr H . (6.11) 
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We express the coupling A,, through A,- and arrive at the representation Eq. (6.1), 
with 

ri = A,(/io;^i)|y In^ + fea} +■■■• (6.12) 

The difference between M and M does not matter at the order of the perturbative 
expansion considered here. 

We now come to the main point: because T^ depends on the scale /xi, 

/ii- — T^ = -a2Ar H , (6.13) 

a/ii 

the amplitude T in Eq. (6.5) becomes convention dependent as well [whereas 
ReT = T + aT^ + ■ ■ ■ is independent of Hi] The relation between the amplitudes 
evaluated from Eq. (6.5) with matching scales /xi and /X2 is 

%, = %,{l + aa^ln^) + ■ ■ ■ . (6.14) 

The ambiguity shows up at order a and does not affect the interpretation of the 
"purely strong" amplitude at leading order, which is extracted by using Eq. (6.2). 

6.2 Isospin breaking effects in ChPT 

In the real world, the dynamics of hadronic atoms can be analyzed with ChPT. 
As this is a quantum field theory, the issue of purification of the amplitude from 
electromagnetic contributions is affected with the same ambiguity as the toy 
example discussed above. Hadronic atoms are not an exception in this sense: 
whenever an attempt is made to split electromagnetic effects in ChPT, one is 
faced with this problem. We now explain how this can be handled. 

We start from the observation that the Lagrangian in ChPT contains, aside 
from the particle fields $, a set of LECs, which are not determined by chiral sym- 
metry alone. We consider ChPT in the hadronic sector, with photons included. 
Generically, one has 

Cc^^T = C{g-lC-M,e-^). (6.15) 

We have classified the LECs in two groups: the ones in Q [strong LECs), which 
stand for those which survive at e = 0, like the pion decay constant in the chiral 
limit, the LECs Lj at order p'^, etc. The group /C stands for the so called electro- 
magnetic LECs that one has to introduce while incorporating electromagnetic 
interactions in the theory, like the pion mass difference in the chiral limit, the Ki 
introduced by Urech [194], etc. In addition, there are the quark masses collected 
in M., and the electromagnetic coupling e. 
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Given the structure of this Lagrangian, Green functions can be evaluated in 
a straightforward manner. The Lagrangian is so constructed that all UV diver- 
gences cancel - the result for any quantity is independent of the renormalization 
scale /io and can symbolically be written in the form 

dX 

^ = 0. 6.16 

d/io 

We have indicated the dependence on the renormalized LECs Q^ , K7 and on 
external momenta p,. In particular, one can determine the algebraic form of the 
masses tt^, 7r°, K^, K^, r], . . .. We now define the isospin symmetry limit in the 
following manner. 

• Set e = 0, niu = trid- The quantities X so obtained depend on the renor- 
malized parameters Q^ , on the quark masses and on the momenta. 

• In the isospin symmetry limit so defined, physical masses are grouped into 
mass degenerate (isospin) multiplets. Assign numerical values to these. In 
particular, for the pion, kaon and proton mass, choose the physical values 

for M7r+, MK+,mp, and for the pion decay constant take F,^ = 92.4 MeV. 

• Isospin breaking terms are defined to be the difference X—X. These depend 
on the full set of renormalized parameters Q^; e; KJ , and, in addition, on the 
physical masses, on quark mass ratios, on F^r and on the momenta. 

To numerically calculate the isospin breaking corrections, one uses measurements 
or estimates to assign numerical values to the needed renormalized LECs, and 
on quark mass ratios, as a result of which the isospin breaking terms can be 
calculated for any quantity. 
We comment the procedure. 

• It is algebraically well defined and internally consistent. As far as we can 
judge, it agrees with calculations of isospin breaking corrections performed 
in recent years by many authors. 

• The problem with the ambiguity in purifying quantities from electromag- 
netic interactions is hidden in the numerical values assigned to the LECs. 
Here, we assumed that all LECs are known, to within reasonable error 
bars. In other words, fixing the values of all LECs unambiguously defines 
the splitting in physical observables. 

• The relation of the chiral Lagrangian to the underlying theory, QCD+QED, 
will not be discussed here [see Refs. [204-210] where certain aspects of the 
problem are considered]. It is however clear that the mentioned ambi- 
guity will show up again once the matching of the chiral Lagrangian to 
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QCD+QED is performed. We take it that the uncertainties in the electro- 
magnetic LECs are chosen in such a manner that this ambiguity is taken 
care of appropriately. 

• Assigning, in the isospin symmetry limit, physical values to F^,M^+,M/^+ 
and to the proton mass overconstrains the available parameters in QCD 
i^QCD^'ftiifns)'- fixing the first three quantities allows one to calculate the 
proton mass, and there is no reason why this mass should coincide with the 
physical value rUp. However, we expect the difference to be small, with a 
negligible effect on the quantities considered here. 

• The manner in which electromagnetic corrections are treated in Ref. [206] 
differs from the present prescription in at least one important aspect: in 
that treatment, the leading-order meson mass terms that enter the strong 
chiral Lagrangian become scale dependent and run with a in the full theory. 
This is not the case in the procedure advocated here. 

We add a note concerning the bookkeeping of isospin breaking corrections. 
Isospin breaking contributions can be evaluated as a power series in a and rrid—mu 
(modulo logarithms). It is useful to define a small parameter 5 as a bookkeeping 
device, 

Sr^a,md — mu. (6-17) 

A particular calculation is then carried out to a specific order in 6. We note that 
the above counting is not the only possible choice. Indeed, in the treatment of 
pionium [36, 40, 42, 44], the bookkeeping was performed differently: 5 ~ a , {md — 
rriu)'^. The reason for this is that, both in the pion mass difference and in the tttt 
scattering amplitudes, the terms linear in m^ — ruu are absent - isospin breaking 
generated by the quark mass difference start at 0{{md — muf'). On the other 
hand, if one still wants to count {nid — niu) like a in pionium, this amounts to 
discarding all strong isospin breaking corrections in the final result, because the 
calculation was carried out at next-to-leading order in 6. We have indeed found 
[40, 42] that numerically, these corrections are tiny, and nothing changes in the 
answer if these are left out. 

This concludes our review on the formulation of the general approach to 
hadronic atoms. We have used the lagrangian in Eq. (4.1) to illustrate the main 
ideas, which do depend neither on the details of the underlying interaction, nor 
on the choice of a particular bound state. In the remaining part of this article 
we consider the application of this general theory to several hadronic atoms, in 
the framework of chiral perturbation theory. 
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Figure 7: Decay channels of various hadronic atoms. Different thresholds are 
separated due to electromagnetic/strong isospin breaking [a ^ 0, m„ ^ m^], 
and/or ^[/(S) breaking [ttt.^ 7^ m], and/or a non-zero value of the pion mass 
[m 7^ 0]. Multi-photon states are not explicitly indicated: e.g., the decay channel 
of pionium into 71*^7?° also implies decays into states with any number of additional 
photons. The scales in the Figure are chosen arbitrarily. 

7 Hadronic atoms and scattering lengths: 
General observations 

Before starting the investigation of the individual atoms listed in Table 1 , we find 
it useful to provide already here some general information on their spectra, and 
on the relation of these to the pertinent threshold amplitudes. We start with an 
investigation of the counting rules in the parameter 5 introduced at the end of 
the last section. 

7.1 Counting rules for the widths 

Strong interaction energy shifts start at order b^ . On the other hand, the counting 
rules are slightly more complicated for the decay widths, because the bookkeeping 
of the phase space factor (in powers of 5) depends on the final state involved. 

The hadronic atoms decay into states whose thresholds lie below the bound- 
state energy. In the cases that we consider here, the separation of the thresholds 
is due to (one or several of) the following mechanisms (see Fig. 7): 
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i) Level splitting in the atoms, which is O(a^); 

ii) Isospin breaking [a 7^ and/or rriu 7^ rrid] ; 

iii) SU{3) breaking [m^ 7^ m]; 

iv) Chiral symmetry breaking [m 7^ 0]. 

The phase space factor in the S'-wave counts as 0(6^^"^) for ii) and as 0(1) for 
iii) and iv). Furthermore, up to order S^, only the decay width of the S'-states is 
non-zero, 

lmAEra = -lTJio + 0{6'). (7.1) 

In the case of particles with spin, AEni denotes the averaged energy shift, see 
appendix A. Transition widths between different energy levels also start at 0{6^) 
(see, e.g., Ref. [211]). 

Taking into account the fact that the square of the Coulomb wave function 
at the origin counts at 6^, the decay widths into different final states are seen to 
obey the following counting rules: 

Pionium 

In the iS-wave, pionium decays overwhelmingly into tt^tt". At leading order, the 
decay width counts at 0(5^/^). Next-to-leading order isospin breaking correc- 
tions to this partial decay are of order S^^^. The anomaly-induced decay into 
7r°7 cannot proceed in the S'-wave, and the decay into 77 starts at 0{6^) [54]. 
At order 0(5^/^), the multi-photon final states (e.g. 7r°7r°7) do not contribute 
either, see Ref. [44] for a detailed discussion. At order S^^'^, only the decay into 
the 7r°7r° final state needs to be taken into account. 

ttK atom 

The decay width of the iiK atom into ii^K^ starts at order 5^/^, and the next- 
to-leading order corrections count as 0(5^/^). Up to 0{6^) there are no other 
contributions, see Ref. [44] for a detailed discussion. At order 5^'^, only the decay 
into the tt^K^ final state needs to be taken into account. 

Pionic hydrogen 

In difference to pionium and to the ttK atom, there are two decay channels 
in ttH with comparable widths: the mr^ and the nj states, with ~ 60% and 
~ 40% decay probabilities, respectively. The ratio of these two quantities gives 
the Panofsky ratio P ~ 1.546 [212,213]. The decay width into nn^ starts at 
0(5^/^) and next-to-leading isospin breaking corrections to it count as 0(5^^^). 
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The decay width into 727 starts at 0(5^) - there are no corrections to the leading 
order at the accuracy we are working. The Panofsky ratio at leading order counts 
as 0(5~^'^). Again, there is no need to consider isospin breaking corrections. 
Further, at the accuracy we are working, the decays into final states with more 
than one photon can be neglected. At order S^''^, only the decays into the n7r° 
and wf final states need to be taken into account. 

Kaonic hydrogen 

Since Mj^o + m^ > M^- + rrip, the decays of KH proceed overwhelmingly into 
the strong channels vrA and ttS. The decay width at leading order counts as 
0{S^). The decay width into S7, A7, as well as next-to-leading order isospin 
breaking corrections to the decay into strong channels start at 0(5^). Other 
decay channels are suppressed. 

Remark: The couplings 7r"'"7r~7r'^7 and K'^TT~K^'~f in the relativistic theory are 
proportional to the antisymmetric tensor gt^'^'^f^. For this reason, the pertinent 
transition amplitudes necessarily contain at least one soft momentum and are 
suppressed as compared to the n^p -^ wj amplitude. There is therefore no need 
to introduce a Panofsky ratio for these channels at the order considered here. 

Pionic deuterium 

Pionic hydrogen can decay e.g. into nn. Therefore, the phase space factor is of 
order 5°, and the width is of the order 5^. 

Kaonic deuterium 

Kaonic deuterium can decay e.g. into Envr, AriTr. The corresponding phase space 
factor is of order (5°, and the width is of order 5^. 



7.2 Spectrum and scattering lengths 

We are concerned in this review to a large extent with the connection between 
the hadronic atom spectra and the underlying hadronic scattering amplitudes at 
threshold (scattering lengths). Four elastic reactions are involved: tttt -^ tttt, 
ttK -^ ttK, ttN -^ ttN and KN -^ KN . In the isospin symmetry limit a = 
0, rriu = rrid, each of them is described by two independent isospin amplitudes. 
These are real in the first three reactions, and complex in the last reaction. The 
atoms considered fall into two distinct groups: 

Group 1: To this group belong the first three atoms displayed in Figure 7: 
A2n,A.,rK, and irH . At the order considered here, these atoms can decay into 
states that are separated in energy from the bound state through isospin breaking 
effects, which are calculable in ChPT for which can be taken into account with 
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the Panofsky ratio in case of ttH ). For each of these three atoms, there are two 
real threshold amplitudes in the isospin symmetry limit. Therefore, it suffices to 
have two experimental numbers to pin them down: the strong energy level shift 
and width of the ground state. 

Group 2: To this group belong KH , KD and vrZ) . These atoms can decay 
into states that are separated in energy from the bound state also in the absence 
of isospin breaking, because the splitting is due to m 7^ and(or) m ^ rris. As a 
result of this, the widths are enhanced. We already mentioned in subsection 2.4 
that the KN scattering lengths cannot be determined from the energy level shift 
and width of the ground state of KH even in principle, because one has to fix two 
complex numbers. As we will see later, it turns out that KD provides the lacking 
information, see subsection 12.4. The case of irD plays a special role here. From a 
principle point of view, information on this compound is not needed: to measure 
the two real pion-nucleon scattering lengths, knowledge of the strong energy level 
shift and width of the ground state of nH suffices. On the other hand, due to an 
insufficient knowledge of one of the LECs involved, information aside from irH is 
presently needed. It is provided e.g. by the strong energy shift in nD , which 
can be related to the pion-nucleon scattering lengths. On the other hand, its 
width cannot be used for this purpose. The physical reason for this is the fact 
that ttD can decay e.g. into nn final states, a reaction which cannot be related 
to the elastic pion-nucleon scattering amplitude. The technical reason behind 
this is the fact that the ratio between the strong energy level shift and the width 
in ttD is of order S^, while the imaginary part of the pion-nucleon amplitude is 
of order S^^^: there can be no DGBT formula that relates the energy shift and 
width directly with the pion-nucleon scattering amplitude, because there would 
be a mismatch in counting powers of 6. 



8 Pionium 

8.1 DIRAC experiment at CERN 

As already mentioned, the aim of the DIRAC collaboration at CERN is to mea- 
sure the lifetime of the tt^tt^ atom in the ground state with 10% precision [72- 
78] . This allows one to determine the difference | oq — 02 1 of 5*- wave vrvr scattering 
lengths at 5% accuracy. Measurements of the lifetime of tt+tt" atoms have also 
been proposed at J- PARC and at GSI [214,215]. For an earlier attempt to mea- 
sure pionium production, see Ref. [216]. 

Details of the set-up of the DIRAC experiment at the CERN Proton Syn- 
chrotron can be found in Ref. [73]. The underlying idea is the following. High- 
energy proton-nucleus collisions produce pairs of oppositely charged pions via 
strong decays of intermediate hadrons. Some of these pairs form tt+tt" atoms 
due to Coulomb final state interaction. Once produced, the 7i~^n~ atoms prop- 
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agate with relativistic velocity. Before they decay into pairs of neutral pions, 
the atoms interact with the target atoms. This interaction excites/de-excites or 
breaks them up. The vr'^vr" pairs from the break-up exhibit specific kinematic fea- 
tures, which allow one to identify them experimentally. Excitation/de-excitation 
and break-up of the atom compete with its decay. Solving the transport equa- 
tions for excitation/de-excitation and break-up leads to a target-specific relation 
between break-up probability and lifetime, which is believed to be known at the 
1% level [77]. Measuring the break-up probability then allows one to determine 
the lifetime of pionium [79, 80, 217]. 

The first observation of vt+tt^ atoms [218] has set a lower limit on its lifetime, 
r > 1.8 X 10^^^ s (90% CL). Recently [77], the DIRAC collaboration reported a 
measurement based on a large sample of data taken in 2001 with Ni targets, see 
Eq. (2.1). 

As pointed out Refs. [86, 87], a measurement of the energy splitting between 
2s- and 2p-states of pionium allows one to determine the combination 2ao + 02 
of scattering lengths, see also section 2 of this report. Assuming 09 — 02 > 0, 
knowledge of the width and energy shift thus allows one to determine separately 
Oq and 02. 

8.2 Two— channel problem 

In contrast to scalar QED considered in previous sections, one is concerned here 
with two coupled channels of non-relativistic particles, tt+tt" and 7r°7r°, which 
are separated by a mass gap M^ — Mt,o -C M^, see Figure 7. As a result of this, 
the ground state of pionium dominantly decays into 7r°7r°, whereas its decay into 
77 is suppressed by a factor ~ 10^. 

We investigate the decay process in the non-relativistic framework developed 
in previous sections. In particular, the master equation (5.46) and the decompo- 
sitions Eqs. (5.47, 5.48) apply also here - we simply need to adapt the potential 
V that occurs in the quantity f"', see Eq. (5.15). This can be achieved by con- 
structing the relevant Hamiltonian. We follow Refs. [36,40,42,44], which use a 
two-channel formalism. An alternative method, based on a one-channel frame- 
work, will be invoked in the description of pionic hydrogen later in this report. 

The Lagrangian of the system contains the non-relativistic charged and neu- 
tral pion fields 'k±{x), vro(x) and the electromagnetic field A^{x). At the accuracy 
we are working, it suffices to retain only the following terms (cf. Eq. (4.2)), 



T I (0(DE - ED) + ■ ■ -^vTi + 4[idt - ^.° + ^ 



A^ 



TTO 



8M3o 

+Ci7ri7rl7r+7r_ + C2{-k\tiItiI + h.c.) + C3(4)V^ + ■ ■ ■ , (8.1 
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where ci, C2, C3 denote non-relativistic four-pion couplings, and (r^) is the charge 
radius of the pion. As is shown in Ref. [44], the second-order derivative term 
retained in Ref. [36,42] is in fact redundant and can be ehminated by using 
the equations of motion. This procedure, which merely amounts to choosing 
the coupling constant C4 = in Refs. [36,42], does not affect any of observable 
quantities. 

Without performing any calculation, the electromagnetic energy shift in the 
bound state with angular momentum / and principal quantum number n up to 
terms of order a^ can be directly read off from Eq. (5.30). One finds [44] 

^^-^ -"^-(,8;^ + 64;^"2n3(2/ + l)J+ 6n3 ^'° " ^^'^^ 

Further, the contribution due to vacuum polarization AE^^'^ is defined by 
Eq. (5.41) and is calculated in Ref. [200], see Eq. (5.42). 

In order to evaluate the strong shift of a given level, we return to Eq. (5.46) 
and perform the calculations by explicitly assuming e = everywhere, except 
in the static Coulomb interaction (this is a perfectly legitimate procedure up to 
terms of order 6^). Suppressing everything but four-pion local interactions and 
the relativistic mass insertions, one has 

r"0(E„) = V + VGg.°V + VG^°VG^°V + VGg.°VG^°VG^°V + • ■ ■ , (8.3) 

where V = H5 + Hji, and 



He 



/ c/^x < Ci7r^7rl7r+7r_ + C2(7r|7rl7ro + h.c.) + C3(7rJ)^7ro '> 



Hh = -/'''-{E^I^-.+^g^-o}^ (8.4) 

Note that at the accuracy we are working, three interactions in Eq. (8.3) are 
sufficient (each iteration is suppressed by one power of S) . The general expression 
for the strong energy shift was already worked out in Eq. (5.19). Here, it remains 
to simply calculate the pertinent matrix elements. 

As an illustration, let us calculate the (complex) level shift at lowest order in 
6. It is clear that one should obtain the DGBT formula for the real and imaginary 
parts of the shift. In order to demonstrate that this is indeed the case, it suffices 
to consider only those contributions to the quantity t"'^{z), which are shown in 
Fig. 8 (the tree diagram and the neutral pion loop). We obtain 

{z - EX' = -|^no(0)P(ci + 2clUz)) + o{5'/^) , (8.5) 

where the quantity Jo(^) is the loop integral defined in Eq. (3.16), with the mass 
M replaced by M^^o. This function has a branch point at 2 = 2M^o and its 
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Figure 8: Lowest-order contributions to the complex energy shift of the tt+tt" 
atom, see Eq. (8.5). The labels "+", "— " and "0" (solid and dashed lines) denote 
the charged and neutral pions, respectively. 

imaginary part has the same sign as the imaginary part of z throughout the cut 
plane. Therefore, Eq. (8.5) has no solution on the first Riemann sheet. On the 
other hand, if we analytically continue Jo{z) from the upper rim of the cut to the 
second Riemann sheet, we find a zero at 

Re z = En T- ci H , 

Imz = _^M^,l + ..., (8.6) 

where p„ = Mn-o(£'„ — 2Mn-o). Further, matching in the isospin limit yields 

SM^ci = 47r(2ao + as) + ■ ■ ■ , 

3M^C2 = 47r(a2 - ao) + ■ ■ ■ , 

SM^cs = 27r(ao + 203) + ■ ■ ■ , (8.7) 

where the ellipses stand for terms of order a and {m^ — m^)^. Using |\E'„o(0)P = 
a^ M^ / [Sirn^) , we finally conclude that i) Eq. (8.6) reproduces the leading order 
expressions Eqs. (2.2) and (2.5), and ii) the bound-state pole is indeed located on 
the second Riemann sheet. Note that the real and imaginary parts of the energy 
shift are of order S^ and 5^/^ , respectively. 

Calculations up to and including first non-leading order in S are straight- 
forward and proceed in complete analogy to the case of scalar QED, considered 
earlier - details may be found in Refs. [36,40,42,44]. The result is 

ry^M 9 
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with 

A a 

Kn = ^7^(ao + 2a2)^- -s„(a) (2ao + a2), 

Ct 

K'n = -- Sn{a) {2ao + a2) , 
p: = \'/'{ElM^o,M^o)/i2En), /^^ = Ml-Mlo, (8.9) 

where Sn(a) is given in Eq. (5.32). The quantities Ax and Ac coincide with the 
real parts of the relativistic threshold amplitudes in the channels vr+vr" — >■ 7r°7r° 
and TT+TT" -^ TT+vr", respectively, up to kinematic normalization factors. The 
precise relations are given in appendix A. With the normalization chosen there, 
one has 

Ax = aQ-a2 + 6Ax + o{5), A = 77 (2ao + ^2) + 5A + o(5) , (8.10) 

where 5Ax and 5 Ac denote isospin breaking corrections of order 5. Finally, 
including vacuum polarization correction in the above expressions amounts to 
Kn-^ Kn + C and K'^ -^ K'^ + 5™^ where C^ is defined by Eq. (5.44). 

For completeness, we also display the strong shift of the P-wave states at 
leading order [44], 

AE- = -^4^«^M,X, (8.11) 

where a\ denotes the P-wave vrvr scattering length. The strong decay into n^n^ 
final states from P-wave states is forbidden by C-invariance. 

The equations (8.8), (8.9), (8.10), (8.11) provide the framework for the nu- 
merical analysis of pionium observables. 

8.3 DGBT formula and numerical analysis 

The relations Eqs. (8.8), (8.9) and (8.10) can be rewritten as follows, 

Ai^f = -^(2ao + a2)(l + + o(^'),^; = /^ + <, 
bn-^ IttQ + 02 

9n'^ Oq — a2 



(8.12) 

At this stage one may invoke ChPT and evaluate the quantities 5 Ax and 5Ac- 
This procedure leads to a systematic calculation of the isospin breaking correc- 
tions. 
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We now demonstrate how the calculations at this concluding step are carried 
out. In general, the following representation for the quantities 6Ax and 6Ac holds 
in ChPT, 

5Ax = hi{md-muf' + h2a, 5Ac = h'i{md- muf + h^a, (8.13) 

where hi,h[ are functions of the quark mass m and of the RG-invariant scale 
of QCD. According to our definition of the isospin limit, the parameter rh is 
adjusted so that in this limit, the pion mass is equal to the charged pion mass 

Consider, for instance, the amplitude Ax at tree level in ChPT, 

A = ^^(4M.' - M^o) + 0{p\ eV) . (8.14) 

In the isospin symmetry limit this expression reduces to 

ao-a2 = ^|^ + 0(p^). (8.15) 

Comparing this with Eq. (8.14), we find 

3A 

A = ao - a2 + ^1^ + 0(p^ eV) . (8.16) 

From this result, we may read off the coefficient /12 at leading order in the chiral 
expansion, 

3(M|-M|^ (8.17) 

[To be precise, the first term on the right-hand side of this equation should be 
evaluated at a = 0. To ease notation, we omit this request here and in the 
following]. On the other hand, the above calculation is not accurate enough to 
determine hi at leading order, because for this purpose, the amplitude is needed 
at order p^. This procedure may obviously be carried out order by order in the 
chiral expansion - all that is needed is the chiral expansion of the scattering 
amplitude at threshold, at m^ 7^ m^, a 7^ 0. As a result of this, the quantities hi 
are represented as a power series in the quark mass rh (up to logarithms). 

The evaluation of the amplitude for tt+tt^ -^ tt^tt^ was carried out at 
0(p^, e^p^) in Refs. [40,42,195], and a detailed numerical analysis of the pio- 
nium decay width was performed in Refs. [40,42]. This analysis results in the 
following value for the ground-state correction. 



(5A = (0.61 ±0.16) X 10^^ i^i = (1.15 ±0.03) X 10"^ (8.1^ 
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AE^f [eV] 


AE:r [eV] 


AE^r [eV] 


lO^'rni [s] 


n=l, /=0 
n=2, /=0 
n=2, 1=1 


-0.065 
-0.012 
-0.004 


-0.942 
-0.111 
-0.004 


-3.8 ±0.1 

-0.47 ±0.01 

~ -1 ■ 10~6 


2.9 ±0.1 
23.3 ±0.7 
~ 1.2-10^ 



Table 2: Numerical values for the energy shift and the lifetime of the tt+vt atom, 
taken from Ref. [44]. The quantity AE^^J stands for AE.^*'', and r„o = (r„)"^ 

Using input scattering lengths from Refs. [13, 14], we obtain for the isospin break- 
ing correction 6i in Eq. (8.12) and for the lifetime of the ground state. 



r = r^i = (2.9 ±0.1) X 10 



-15 



6i = (5.8 ±1.2) X 10" 



^.19) 



We note that ^i amounts to a 6% correction to the leading-order DGBT for- 
mula [8]. On the other hand, the preliminary result of DIRAC, displayed in 
Eq. (2.1), corresponds to the following value of the difference ao — 02, 

-0.033 /M—l 



ao - 02 = 0.264 



-0.020 



^.20) 



In Ref. [44], the energy shift of the ground state of pionium was worked out in 
an analogous manner (see also Ref. [219]), with the result 



AEf 



-3.8±0.1)eV, 5; = (6.2 ±1.2) X 10" 



(8.21) 



As was discussed in the previous sections, the width and strong energy shift are 
modified by vacuum polarization effects, (5„ -^ 5„ ± 5™'^, 6'^ -^ 5[^ ± 5^^, where 
^.^'^ is defined in Eq. (5.44). For example, in the ground state, the correction 
amounts to 5\^ = 0.31 ■ 10~^ [200]. We conclude that it is safe to neglect, in this 
system, the 5™^ altogether: the uncertainties in 5n and 5'^ are larger than 5^^ 
itself. 

Finally, in Table 2 we list different contributions to the energy shift and 
the lifetime for the first few levels in pionium [44], calculated by using the input 
scattering lengths from Refs. [13, 14]. The calculations for the S'-wave states were 
carried out at next-to-leading order in isospin symmetry breaking. The bulk of 
the uncertainty in these quantities is due to the uncertainties in the pertinent 
scattering lengths. The lifetime of the 2p-state is calculated at leading order only 
and is determined by the 2p— Is radiative transition [87]. Finally, the theoretical 
value for the 2p — 2s energy splitting is given by [44] 



AE' 



-2p 



AEf 



AEl^ 



AE^T 



AEl^ 



AEJf 



-0.59 ± 0.01 eV. (8.22) 



The uncertainty displayed is the one in AEf''' only. To the accuracy we are 
working, we may neglect the strong shift in the 2p state, because it is of order 



a 
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8.4 Two— photon decay of pionium 

The decay width of the pionium ground state in two photons is given by 

r?^ = ^a^|A(4M2,-M2,-M2)|2 + o(5^), (8.23) 

where the invariant amphtude A[s, t, u) for the process 77 — »• tt+tt" is defined 
in Eq. (2.4) of Ref. [220]. This amphtude was evaluated at NLO in chiral 
SU{3) X SU{3) by Bijnens and Cornet [221], and to two loops in SU{2) x SU{2) 
in Refs. [220, 222]. Here, we use the 5t/(2) x SU{2) version displayed in Eq. (5.1) 
of Ref. [220], 

2 M^ / - J — J \ 

Ml A{4Ml -Ml -Ml) = 1 + ^;f G^{AMl) + ^—^ + 0{Mt) , (8.24) 



487r2 



where 



1 /, 2Ml r^ dx 



s 



''-W = -l6^(^ + ^XT'"(l-Mj^(l-^»j' 

G.(«f.') = ^- (8-25) 

The terms of order M^ in Eq. (8.24) denote two-loop contributions, and /s, Iq are 
O(p^) chiral LECs [84]. These are related to the pion polarizabilities Ott and (3.,^, 
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The two-photon decay width of pionium was evaluated already earlier [54, 
192]. Here, we comment on the expression given by Hammer and Ng [54], who 
were the first to incorporate contributions from the pion polarizabilities. Their 
expression does not contain the contribution from G^r, which is, in the chiral ex- 
pansion, of the same order as the one from the polarizabilities. Numerically, this 
difference amounts to a large effect. Using Tg — Ts = 3.0 ± 0.3, which corresponds 
to (a - 13)^ = 6.0 • 10^^ fm^ [220], we obtain from Eq. (8.24) 

Ml A{4Ml -Ml, -Ml) = 1 + 4.2 ■ lO'^ + 2.9 ■ 10"^ + 0{M^) , (8.27) 

i.e., the contribution from G7r(4M^) is larger than the one from the pion polar- 
izabilities. Note also that in Ref. [54] much larger values of the polarizabilities 
a^ = — /3^ = (6.8 ± 1.4 ± 1.2) ■ 10~^fm have been used, tending to mask the 
absence of G^{AMl). 

The reason why the contribution from Gtj{s) was missed in Ref. [54] is the 
following. In that article, the annihilation amplitude 77 -^ tt+tt^ is obtained from 
the Compton amplitude through analytic continuation. Moreover, the authors 
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expand the amplitude in small photon momenta and retain only terms which 
are at most quadratic in this expansion. In the annihilation channel, however, 
the photon momenta are of the order of the pion mass. Therefore, the analytic 
continuation of the truncated Compton amplitude can be justified if and only if 
the expansion of this amplitude in photon momenta has a convergence radius 
of the order of the heavy scale in ChPT. In this case, the coefficients of such an 
expansion are regular in the chiral limit. However, as can be seen from Eq. (8.25), 
this is not what happens: at the Compton threshold, the expansion is carried out 
in powers of s/M^ (with s ^ t). The coefficients of the higher-order terms in 
the expansion of G{s) become more and more singular in the chiral limit. As a 
result of this, they contribute at the same order at s = 4M^ and must thus be 
kept. 

We summarize this section with the observation that the vt+tt" atom is now 
completely understood in the framework of QCD+QED, on a conceptual as well 
as on a quantitative level. For bound-state observables, the theoretical predic- 
tions are made with percent accuracy. If the experiments are performed with the 
planned accuracy, a precise measurement of the tttt scattering lengths is indeed 
feasible. 



9 ttK atom 

As far as the description of the bound state within the non-relativistic effective 
theory is concerned, the ttK atom problem is completely analogous to the case 
of pionium considered in the previous section. For this reason, we shall display 
only the final results here. Details may be found in Ref. [44]. 

For definiteness, we consider the bound state of K~^ and tt^. The strong 
energy shifts and widths of the ttK atom are given by expressions similar to 
Eq. (8.8), 






Tn = ^pMl(l + ir„) + 0(5^/2), (9.1) 



n^ 



where Ax and Ac are now related to the relativistic threshold amplitudes for 
the scattering processes n~K~^ -^ -k'^K'^ and tt^K^ -^ 7r~K~^, respectively, see 
appendix A. In particular. 

Ax = a-o + SAx + o{S) , Ac = a^ + Uq + SAc + o{S) , 

«o = o ' "o = o ' (9.2) 
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where a'l'^^c?!'^ are the -nK scattering lengths with total isospin / = 1/2 and 
/ = 3/2, respectively, and bAx.^Ac denote isospin breaking terms of order b. 
Furthermore, 






e; ^"° 



2a/ic(a^ 



% )SnKa) 



(9.3) 



where A^ = M"^ — M^o, S+ = M^, + M^, ^ic is the reduced mass of the vr K^ 
system, and 



En 



a'^ldc 



+ 2n2 ' 
The equations (9.1) can be rewritten as 

2aV^ 



p: = A^/^(^^,M^o,M^o)/(2E„). 



(9.4) 



AEl 



n 



n 



3 (a+ + ao)(l + <5;) + o(5^) 
^(l + 5„) + o(5^/^). 



5' 



5A, 



«o + «o 



+ X 



3, ,2 



5n 



2 (5^,, 



+ ifn 



(9.5) 



The electromagnetic shift is given by 



Ai?:r 




^. 3/ie\ 

" 1 , 1 


"3 1 " 

_8n 2/ + l_ 

3 


4aV?A 
+ n3 




n3(2/ + l)_ 


5 



5W 



(9.6) 



where A = | ((?"^) + {^k))- Vacuum polarization introduces a small change in 
Eq. (9.5): 5„ — *> 5„ + 5™'^, (5^ — > ^^ + ^™'^- As in pionium, this correction is very 
small: for example, in the ground state it amounts to 61^^^ = 0.45 ■ 10"^ [200]. In 
the numerical analysis given below this tiny contribution is neglected. 

Finally, the strong shift in the P-wave state and the P-wave decay width into 
ir^K^ at leading order are given by [44] 



str 



AP- 



■■- nl 



2(n' 



n'' 



in' -I 



a Vc [o-t + «! ) 



n 



5 «v^(p:)'(«r), 



(9.7) 



where a^ denote P-wave tiK scattering lengths. 

A comprehensive numerical analysis of the ttK atom observables in ChPT 
has been carried out in Ref . [44] . Below we give a short summary of the results 
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102^1 


102^2 


102(51 


W5'^ 


4.0 ±2.2 


3.8 ±2.2 


1.7 ±2.2 


1.5 ±2.2 



Table 3: Isospin breaking corrections in the DGBT-type formulae for the tiK 
atom in ns states, taken from Ref. [44]. See Eq. (9.5) for the definition of 5i and 

of this analysis. The isospin breaking corrections to the ttK threshold ampli- 
tudes Eq. (9.2) have been worked out by several authors [223-227] at 0(p^, e^p^). 
Whereas the analytic expressions for 5Ax and 5Ac obtained in Refs. [223-227] 
are not identical, the numerical values agree within the uncertainties. In the 
calculations of the the ttK atom observables carried out in Ref. [44] , the values 



5A = (0.1 ±0.1) X lO-^M; 



5Ac = (0.1 ±0.3) X 10~2m^ 



(9.8) 



were used. Table 3 contains the final results of these calculations for the isospin 
breaking corrections to the strong shift and width of the ns state with n = 1,2. 

Next, we consider the calculation of the lifetime of the tcK atom ground 
state, given the input for the S-wave tiK scattering lengths. As discussed in 
section 2, different values for scattering lengths are available in the literature. In 
the analysis carried out in Ref. [44], the results of O(p^) calculations in ChPT [104, 
228], as well as the solutions of the Roy-Steiner equations [106] have been used 
as an input. The prediction for the lifetime is^'^ 



r= (4.8 ±0.2)- 10"^^ s 
r= (3.7 ±0.4) -10"^^ s 



ChPT O(p^), 
Roy-Steiner equations [106] . (9.9) 



Finally, in Table 4 we display various contributions to the energy shift and width 
of the low-lying levels of the ttK atom, taken from Schweizer [44]. Moreover, 
because the strong shift in the P-wave is small, one obtains for the 2s — 2p level 
splitting [44] 

-2p 



AE^ 



AEf ± AE: 



cm 
20 



AE^T 



± AE™'^ - AEJf = -1.4 ± 0.1 eV . (9.10) 

We summarize that, as a result of the investigation carried out in Ref. [44], 
the next-to-leading order isospin breaking corrections to the observables of the 
TiK atom in the ground state are now known with a precision that is comparable 
to the one in pionium. In our opinion, this is a completely satisfactory situation 
from the point of view of the data analysis of possible future experiment. We 
refer the reader back to subsection 2.2 for comments concerning the relevance of 
measuring nK scattering lengths. 



""^•^ Schweizer [44] does not quote an uncertainty for the hfetime which result from ChPT 
at order p^. For completeness, we re-evaluated the first line of Eq. (9.9), using the value 
a^ = 0.079 ± 0.001 Af-\ given in Ref. [44]. 
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AE^f [eV] 


AE:r [eV] 


A^:r [ev] 


lO^V^z [s] 


n=l, /=0 
n=2, /=0 
n=2, 1=1 


-0.095 
-0.019 
-0.006 


-2.56 
-0.29 
-0.02 


-9.0 ±1.1 

-1.1 ±0.1 

~ -3 X 10^6 


3.7 ±0.4 

29.4 ±3.3 

~0.7x 10^ 



Table 4: Numerical values for the energy shift and the lifetime of the nK atom, 
taken from Ref. [44]. The quantity AE^^J stands for AE^'^, and r„o = (rn)"^- 

10 Pionic hydrogen 

10.1 Pionic hydrogen experiments at PSI 

During the last decade, the Pionic Hydrogen Collaboration at PSI performed 
high-precision measurements of strong interaction parameters of both, pionic 
hydrogen (ttH) and pionic deuterium (vrZ) ) [1,109-118] in order to extract in- 
dependent information about the S'-wave vrA^ scattering lengths. In this section, 
we mainly discuss pionic hydrogen. 

This exotic atom is formed, when the kinetic energy of a negative pion is of 
the order of a few eV. The incoming pion is captured by the Coulomb field in 
highly excited states and a de-excitation cascade starts, which proceeds through 
X-ray emission as well as through other (non-radiative) mechanisms. The atomic 
cascade ends in the ground state, which then decays, mainly in wy and nn^ states. 
The experimental setup uses the high-intensity low-energy pion beam irEd at 
PSI. The energy shifts are extracted from the measured X-ray lines, see Fig. 9. 
In particular, the strong interaction shift in the ground state is obtained from 
the measured mean transition energy E^p^ig. 

The theoretical framework is analogous to the one summarized in subsection 
5.7, see also below. The quantity -Esp-is is related to the energy levels by -Esp-is = 
Re (2:31 — 2:10), where 



Znl 



En + AEnl, X., 



1 



nl 



2(2/ ±1) 






i)x, 



nlj 



X = z,AE. (10.11 



Here, En stands for the pure Coulomb energy, j = / ± | denotes the total angular 
momentum of a given eigenstate, and AEnij is the energy shift of a generic 
energy level Znij of irH , labeled by the quantum numbers n,l,j. Further, in 
analogy with Eqs. (5. 46), (5. 47) and (5.48), this energy shift is further split into 
an electromagnetic piece, a contribution from vacuum polarization and from the 
strong shift. 



AE, 



nlj 



ae:t, + AE:r 



5,o(AEf --r„,)±o(5^ 



;io.2) 
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Figure 9: The scheme of energy levels of txH with the 3p — Is X-ray transition. 
The displacement and the finite width of the ground state are explicitly shown. 

The strong shift of the ground state is determined from 

A^f = E'?^^^^ + Ell'Lis - -Esp-is , (10.3) 

where -E'lJ^i^ and E^^tis are the (theoretical values for) the electromagnetic tran- 
sition energy and the vacuum polarization contribution, 

^sT-i. = E,-E, + AE3T - AE,7 , EllUs = Ai?3T - Ai?iT • (10.4) 

Here, AE'^™ is defined through AE'^, according to Eq. (10.1). In addition, one 
makes use of the fact that the strong-interaction shift of the 3p state is suppressed 
by additional powers of a and is therefore negligible (see, e.g., Ref. [229]). 

The following values for the strong energy shift and width of the ground state 
have been reported back in 2001 [114], 

AEf = -7.108 ± 0.013 (stat) ± 0.034 (syst) eV , 

Ti = 0.868 ±0.040 (stat) ±0.038 (syst) eV. (10.5) 

One determines the S-wave ttN scattering lengths aj^_ and a^^. from 
AEf = -2«V'(4+±ao+)(l + ^'i) + o(5'), 

Ti = 8a=^/i,M(l + ^)[«o+(l±5i)]'±o(5^/'), (10.6) 

where P denotes the Panofsky ratio, and fic is the reduced mass of the ir'p 
system. Further, p* = X-^^'^{Ef, m^, M^o)/(2-E'i) is the CM momentum of the ir^n 
pair after decay. 
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Figure 10: Determination of the vrA^ S-wave scattering lengths from the energy 
shift and width of ttH : i) Old analysis (bands 1 and 2), which is using the data 
of Ref. [114] given in Eq. (10.5) and applying the isospin breaking corrections, 
calculated in the potential model [230], see Eq. (10.7); ii) New analysis (bands 3 
and 4), which is using preliminary data of Ref. [118] given in Eq. (2.8), together 
with the isospin breaking corrections calculated in ChPT, see Eqs. (10.42) and 
(10.44). The scattering lengths evaluated [132] from the ttD energy shift Eq. (2.9) 
are shown in band 5. The cross indicates the solution Eq. (10.8) for the scattering 
lengths in the old analysis. 

The quantities 6[ and 6i stand for the isospin breaking corrections, which 
in ChPT are of order a or m^ — m„. The lowest-order DGBT formula has 
S'l = 6i = 0. In the analysis of the data given in Ref. [114], one has used the 
values obtained in a potential scattering model [230], 



5[ 



-2.1 ±0.5) -10' 



5i 



-1.3 ±0.5) -10' 



;io.7) 



Recently, a new calculation of the isospin breaking corrections within a potential 
model has been performed in Ref. [231]. The pertinent scattering lengths still 
contain residual electromagnetic effects that cannot be removed with the potential 
model used in Ref. [231] (of course, the same is true for the analysis of Ref. [230]). 
We refer the reader to Ref. [231] for more details. 

The analysis of the experimental data proceeds as follows. Equations (10.5) 
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and (10.6) generate two constraints on the possible values of the scattering lengths 
aQ_^ and a^^_. They are displayed with the two linear bands 1,2 in Fig. 10 (here, one 
assumes that Cq^ is positive). The widths of the bands are determined by both, 
the experimental error in AE'f'" and in Fi, and by the theoretical uncertainty 
in Eq. (10.7)^^. In addition, in the same Figure we also plot the band which 
emerges [132] from the measurement of the strong energy shift Eq. (2.9) of irD, 
see comments in subsection 2.3. 

In the analysis carried out in Ref. [132], all three bands happen to intersect 
in a very small domain, suggesting a coherent overall picture. The S'-wave vrA^ 
scattering lengths can then be determined very well from the intersection domain 
of these three bands: 

a+^ = (-0.0034 ± 0.0007)M^^ , ao+ = (0.0918 ± 0.0013)M-^ Ref. [132]. 

(10.8) 

This picture, however, has undergone substantial changes in the last few years, 
on the experimental as well as on the theory side. First, upgrading of the experi- 
ment has resulted in an increased accuracy of the measured nH width, as well as 
in substantially decreasing its central value, compare Eqs. (10.5) and (2.8). Sec- 
ond, as already stated in the introduction, results like those given in Eq. (10.7) 
cannot be trusted ab initio: the potential model, which was used to derive these 
results, does not include the full content of isospin breaking in QCD+QED (see 
also Ref. [231]). We will discuss this in a quantitative manner later in this section. 
Here, we simply note that, if one applies the isospin breaking corrections calcu- 
lated in ChPT [41,43,45] to the iiH energy and width given in Ref. [118] and 
Eq. (2.8), the central values of the scattering lengths are shifted, see the bands 3 
and 4 in the Figure. There is now no more a common intersection area between 
the TiH and the irD bands, hinting at an internal inconsistency of the theoretical 
methods which were used to analyze the data. The most radical proposal for 
circumventing the problem is to exclude the ttD data from the global analysis, 
thus avoiding the potentially largest source of the systematic error which stems 
from a poor control of the multiple-scattering series for the ird scattering length. 
It was argued that, if the planned 2% accuracy for the width is achieved, it must 
be possible to determine the vrA^ scattering lengths from ttH alone [119]. We will 
investigate this possibility later in this section. 

10.2 Effective theory and counting rules 

The case of nH which is considered here, differs in two aspect from pionium and 
from the irK atom: first and most importantly, it decays with a ~ 40% prob- 
ability into the channel wf, where the CM momentum of the decay products is 



^''We treat the uncertainties in the isospin breaking corrections Eq. (10.7) as systematic ones, 
combine them quadratically with the systematic uncertainty in the energy shift and width in 
Eq. (10.5), and add to the result the statistical uncertainties in Eq. (10.5) linearly. 
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of order of the pion mass. It is clear that, if one tries to include 727 interme- 
diate states explicitly, this will upset the whole power-counting scheme of the 
non-relativistic approach, where M^ plays the role of a hard scale [45]. For this 
reason, we apply a method already used in section 5 for the calculation of the 
two-photon decay width: we construct a non-relativistic theory in which the 77,7 
channel is "integrated out" , and where its contribution appears only through the 
low-energy effective couplings. One could then study tiH in the non-relativistic 
theory with two channels ir'p and ir^n [45], which closely resembles the pio- 
nium case. However, we choose a different strategy here. Namely, in order to 
demonstrate the flexibility of the approach, we find it instructive to consider the 
inclusion of the rni^ channel into the list of "shielded" channels as well and work 
with a one-channel theory. The final result must of course agree with the one 
obtained in the two-channel setting [45]. 

The second generalization is related to the fact that the nucleons are fermions. 
However, the inclusion of particles with spin into the non-relativistic framework 
is nearly trivial and proceeds straightforwardly. 

In order to establish the Lagrangian of the system, one has to merely write 
down all possible operators with a minimal number of derivatives, which are 
consistent with all symmetries (see also Refs. [37,173,174]). The result is 

1 J p2 p4 

C-NR = --F^,F^'- + ijUiVt-mp + ^^ + —^ + --- 
4 [ Znip 8mj^ 

p eaB ^ e(PE - EP) ^ ie(T{V x E - E x P) 1 

v- t f .. .. D2 D^ pe(DE-EDl 



^' mi +-^"^ 



(10,9) 

where if) stands for the non-relativistic proton field, D and D denote the co- 
variant derivatives acting on the proton and the charged pion, with ip'^ V"^ ip = 
ip'^{V'^ip) + {V'^ip^)ip (similar for the pions) and Cp,Cp,Cp,c^,gi,ei are various 
non-relativistic couplings. In the 4-particle sector, we discard the vertices that 
correspond to the P-wave interactions, as well as spin-flip terms, since these do 
not contribute at the accuracy we are working. 

The effective theory based on the Lagrangian Eq. (10.9) enables one to eval- 
uate the energy shift and the decay width of pionic hydrogen perturbatively in 
the isospin breaking parameter 5 ~ a ~ (m^ — m^j). The power counting rules in 
this parameter are as follows (see also subsection 7). The leading-order strong 
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energy shift in Eq. (10.6) is of order S^, and the correction term S'l counts at 0{S). 
In the formula for the decay width, the CM momentum p^ corresponding to the 
decay into the ir^n final state, counts at 0(6^''^). Further, since the cross section 
for 7r~p — > 77-7 starts at 0{a), one gets 1/P = 0(5^'"^). The correction term 5i in 
Eq. (10.6) counts again at 0{6). Our final aim is to evaluate the energy shift 
up to and including terms of order 6'^, and the decay width up to and including 
terms of order 5^^'^. It can be checked that at this accuracy no other terms than 
already displayed in Eq. (10.9) are needed. 

Next we note that Eq. (10.9) now includes the effective-range term with the 
coupling ei. This term has been discarded previously for scalar QED, as well 
as for pionium and the irK atom. The reason for this is that, if the coupling 
ei stays finite in the limit 5^0, the effective-range term does not contribute 
at the accuracy we are working. It can be however shown that in our case ei is 
singular as 5 — *> 0. In general, the counting of the effective couplings gi,ei,- ■ ■ 
in the parameter 6 changes, when one integrates out the vr^n state. To illustrate 
this, we compare the threshold momenta for the 'yn and ir^n intermediate states 
in the CM frame, denoted by p*(0) and p*(0) respectively (see appendix C): 



_ AV^((m, + M.)^m^,0) _ ^^^,^ 



^-p^n^ : p;(0)= " U"P^-'-^.;."-n.^; ^Q(^o 



TT p ^ 7r°n 






As it will be demonstrated below, in the presence of a shielded vr^'n intermediate 
state, the dependence of the constants g\ and e\ on the parameter b becomes 
non-analytic. For example, unlike the previous cases, unitarity now gives Im^fi = 
0(5^/^) (see appendix C). Moreover, we shall see that Imei = 0(5~^/^) and, due 
to this fact, the effective-range contribution survives at next-to-leading order 
in the decay width of pionic hydrogen. To summarize, more complicated power 
counting rules in the parameter b is the price one pays for using a one-channel 
formalism. 

10.3 Matching 

Let us next consider matching of the parameters of the effective non-relativis- 
tic Lagrangian Eq. (10.9) to ChPT. The matching in the one-particle sector, 
which determines the couplings c^ , c^, c^ and c^, is completely analogous to one 
carried out in subsection 4.5. The matrix elements of the electromagnetic current 
between two proton / two charged pion states at leading order in the coupling 
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constant e are given by the standard expressions 

{p's'\j;^{0)\ps) = eu{p',s')L,F,{Q') + ^ a,,Q''F2{Q')\u{p.s) , 

{^^{p')\J7m7rHp)) = ±e{p'+P),FAQ'), (10.11) 

where Q^ = {p' — p) n- At a small momenta, the form factors can be expanded, 
yielding 

FiiQ') = 1 + ^ (rj) + 0{Q') , F,{Q') = n, + 0{Q') , 

F^{Q') = l + ^{rl) + 0{Q^). (10.12) 

Here FilQ"^) , F2{Q'^) and Ft^{Q'^) are the Dirac and Pauli form factors of the 
proton and the electromagnetic form factor of the pion, respectively (see ap- 
pendix A). Further, (r^) and Kp denote the charge radius squared and the anoma- 
lous magnetic moment of the proton. Expanding the form factors in powers of 
3-momenta and performing the matching enables one to read off the values of the 
low-energy constants at leading order in the parameter 5 [41], 

Cp = 1 + fi^p , Cp = 1 + iKp + 2 ^^p\fp} , 

(10.13) 
c^ = l + 2«:p, c^ = M^{rl). 

The four-particle couplings gi,ei are determined through matching to the 7r~p 
elastic amplitude in the vicinity of threshold, calculated in ChPT. The procedure 
is described in detail in the appendices A and C Below, we merely recall the 
notation. In total, we shall need the pion-nucleon amplitudes in 4 different 
physical channels labeled "c/x/O/n", respectively: 

pvr" -^ pir^ (c) / pn' -^ mr^ (x) / mr^ -^ mr^ (0) / mT~ -^ nvr" (n) . 

Channel n shows up only in the ttD case. One calculates these amplitudes at 
0{a, {rrid — rriu)) in ChPT and obtains the threshold amplitudes %^x,o,n- It is 
convenient to define the real threshold amplitudes Ac^x,o,n, which are proportional 
to Iie%,^x,o,n, see Eq. (A. 29). The normalization is chosen so that, approaching 
the isospin limit, one has 

Ao = a+^ + 6Ao + --- , An = a^+ - a^+ + 6An + ■ ■ ■ , (10.14) 

where SAc,SAx,5An = 0{5), 5Ao = 0(5^^'^) are leading isospin breaking correc- 
tions and the ellipses stand for higher-order terms. The corrections at order 5^^"^ 
emerge due to the unitary cusp. At the order of accuracy we are working, the 
couplings gi and ei are finally expressed in terms of A^ Ax and Aq. 
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10.4 Bound state — electromagnetic shift 

A generic solution of the unperturbed Schrodinger equation with pure Coulomb 
potential, which describes the bound state of spin (charged pion) and spin | 
(proton) particles is characterized by the following set of quantum numbers: the 
principal quantum number n = 1, 2, ■ ■ ■ , the angular momentum / = 0,1,2,---, 
the total angular momentum j = |/ — i|, /+i, audits z-component m = — j, ■ ■ ■ ,j. 
The explicit expression for the wave function is given by 

/d'^a 1 

— -- (jm|/(m - (r)-a) *„i(^_.^)(q) |P, q, a) , 

|P, q, a) = b\7]iP + q, a)a\r]2P - q)|0) . (10.15) 

Here, r/i = mp/T,+, 772 = M^/S_|_, S+ = (m^ + M^), fic = mpMT,/T,+, a"!", b'^ stand 
for the creation operators of the non-relativistic pion and proton, {jm\l{m—a)ja) 
are the pertinent Clebsch-Gordan coefficients, and \E'nim(q) denotes the Coulomb 
wave function in the momentum space. 

The energy shift is given by the counterpart of Eqs. (5.46,5.47), 

AE^i, = (*„,,|t"'^(EO|*„;,) + o{6^) , (10.16) 

where E^, = S_|_ — ^ iicOp' is the Coulomb energy for this state, and where we 
have again used the fact that the matrix element on the right-hand side does not 
depend on m. The quantity T^^^{En) is defined through the equation 

r-'J(2) = V + VG^'^'(z)r"'^'(z) , (10.17) 

where, as before, V denotes the perturbation Hamiltonian (everything except the 
static Coulomb interaction) and the pole-subtracted Coulomb resolvent is now 
given by 

^ J 1,-^vrj z - En- ;t?t- 



2S-I 

According to Eq. (10.2), the total energy shift AE^ij is given by a sum of several 
terms. The electromagnetic shift is [41,46] 

2a^lj,l f Cp c^ \ j n n n 

^ n^ \mpM^ ^ 2^7 \ 2/ + 1 ~ 2j + 1 ~ 2 '° 

. ^^ii..^V (10.10) 



n3 '"V8m2 QMl 
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Corrections 


Notation 


Ref. [230] 


Ref. [41] 


Point Coulomb, KG equation 


pKG 
-^10 


-3235.156 


-3235.156 


Finite size effect (proton, pion) 


771 fin 


0.102 


0.100 


Vacuum polarization, order o? 


pvac 


-3.246 


-3.241 


Relativistic recoil, proton spin and 
anomalous magnetic moment 


771 rel 
^10 


0.047 


0.047 


Vacuum polarization, order a^ 




-0.018 




Vertex correction 




0.007 





Table 5: Contributions to the electromagnetic binding energy of the vr^p atom 
ground state and the correction due to vacuum polarization, in eV. Vacuum 
polarization at order a^ and the vertex correction have not yet been calculated 
within the non-relativistic approach. 

The above analytic calculation neatly reproduces the numerical results of 
Ref. [230] for the ground-state energy. In order to facilitate the comparison, 
we split the electromagnetic shift by introducing the same naming scheme as in 
Ref. [230] (see Table 5), 



-g/^cft 



AE^o" 



nKG 



^in^ + E 



fin 
10 



E- 



rel 
10 ' 



(10.20) 



where 



E- 



KG 
10 



-2/^c« 



1 + -a^ 



E 



fin 
10 



«V^ (0 + (<) 



E 



rel 
10 



Ilea 



l^c 



4S. 



+ 






1 



2k„MI 



(10.21) 



The comparison is displayed^^ in Table 5. As is seen, the agreement is very good 
up to some higher-order contributions, which have not been yet calculated in 



i^We used AU = 139.56995 McV, nip = 938.27231 MeV, rUe = 0.51099907 McV, (r^) = 
(0.849 fm)2, (rl) = (0.657 imf, Kp = 1.79284739, a"^ = 137.035989561, see also Ref. [230]. 
The above value for (r^) somewhat differs from the recent update that can be found in the 
hterature [232]. To make the comparison easy, we however do not change this value. 
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the effective field theory framework. Further, we have checked that the analytic 
expressions for E^ and for EIq agree with the analytic result given out in 
Ref. [233]. [These authors did not determine -Eio]- 

We conclude this subsection the with a few remarks. 



1. Eq. (10.19) provides a compact analytic expression for the electromagnetic 
shifts at order a^, for any n, l,j. 

2. To include the contribution AE^^'^ from vacuum polarization at order a^, 
one may use the integral representation worked out in Ref. [200]. It is 
reproduced for convenience in Eq. (5.42). 

3. One may wish to include still vacuum polarization contributions from order 
a^, and vertex corrections, as displayed in Table 5. Note, however, that the 
uncertainty generated by the LEG /i that occurs in the strong energy shift is 
an order of magnitude larger than vacuum polarization at order a^. Indeed, 
one has [see Eq. (10.31) below] 

Ai^flf, = ^/i. (10.22) 

For /i = —1 GeV~^, this amount to a shift of —0.15 eV. An uncertainty in 
/i of the order of —100 MeV~^ induces therefore a shift in AEf^ which is 
of the same size as the contribution from vacuum polarization at order a^. 
We see no way to pin down /i to this precision in the foreseeable future. 

4. For these reasons, the electromagnetic shifts Eqs. (10.19) and (5.42) provide 
a convenient representation, of sufficient accuracy, that allows one to easily 
incorporate changes in the values of the pertinent parameters (r^)^, ■ ■ • , 
whenever needed. Together with Eqs. (10.2), (10.3) and (10.4), they allow 
one to translate measured energy shifts into strong shifts in an easy manner. 

10.5 DBGT formula for ttH 

Finally, at the accuracy we are working, the strong shifts and widths are (cf. with 
Eq. (5.31)) 



Ai^f - ^ r„ = -^ (g, + 47^ei - ^?(gS°(i?„))) + • ■ ■ , (10.23) 



2 TTTl-^ 

where 7„ = 7/ra, and 



'^"■0/ 



l^(^(^)-'"f 



^ciEn)) = -£r A(/.) + In -^ - 1 + sr^ia) . (10.24) 
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Moreover, in order to arrive at Eq. (10.23), we have used that 

/ W)" P'^^^^^P^ = -7n^n0(0) . (10.25) 

To arrive at this result, we used the Schrodinger equation in momentum space for 
the Coulomb wave function \E'„o(p); and took into account that no-scale integrals 
vanish in dimensional regularization. 

Expressing now the couplings gi and ei through the threshold scattering am- 
plitudes (see appendix C for details), we arrive at our final result 



2aV' 

i3 






Tn = ^^^ (l + ^) A (1 + K^) + o{5''^) , (10.26) 

where p* stands for the relative 3-momentum of the nvr^ pair after the decay of 
the Coulombic bound state, 

Pn = ^ ;V • (10-27) 

Finally, the correction terms at this order are given by 

K'n = -afJ'c (a(j+ + %^)sn{a) + 61^" , 

Kn = -2a/i,(a++ + a^+)s„(«) + 2^,Am(a++)2 + 5r, (10.28) 

where s„(a) is given in Eq. (5.32). For n = 1, this result coincides with the one 
given in Ref. [45]. 

10.6 Isospin breaking corrections and data analysis 

Chiral expansion of the threshold amplitude 

The threshold amplitudes which enter Eq. (10.26) contain itself isospin breaking 
corrections, that we identify in the manner displayed in Eq. (10.14). The cor- 
rection to the bound state energy shift in the notation Eq. (10.6) is then given 
by 

S'n= /;f^- +K- (10.29) 

0-0+ + '^0+ 

The correction 6Ac can be calculated systematically, order by order in ChPT. 
Further, K'^ can be evaluated by using the values of the scattering lengths from 
e.g. Ref. [114] and the vacuum polarization correction from Ref. [200]. For 
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example, for the ground state, this gives A^( = 0.66-10~^ + 0.48-10"^ = 1.14-10"^, 
where SJ^^ = 0.48 ■ 10"^. For the chiral expansion of the isospin breaking part in 
the threshold amplitudes, we write 

6Ac = 6Ac,2 + SAc,3 + 0{p^) . (10.30) 

The first two terms have been evaluated by Refs. [41, 43, 234, 235] (Ref. [235] does 
not provide a complete analytic expression of the scattering amplitude). At order 
p^, the result is 

For 6Ac,3, see Ref. [43]. Unlike in the case of pionium, this correction contains 
nontrivial strong and electromagnetic LECs ci, /i, /2 [236,237] already at lowest 
order. This introduces a sizable systematic error in the theoretical prediction of 
the corresponding contributions to 6'^. In order to proceed further, we have first 
to specify the values of these LECs. 

Strong LECs: q 

There are different options for fixing the values of the LECs in the strong sector. 
In the context of the present problem, it is consistent to determine these con- 
stants from the threshold data on irN scattering, similar to Refs. [43,238,239]. 
Below, we closely follow the method of Ref. [238], using the same conventions 
and notations. We denote the isospin symmetric spin-nonflip vrA^ amplitude 
by D~^{q'^,t), where g^ = A(s,mp, M^)/4s is the square of the CM momentum 
and s, t, u are the usual Mandelstam variables. In the vicinity of threshold the 
quantity ReL'"'"(g^, t) can be expanded in Taylor series. 

Re D+(g2, t) = D+o + D+ g' + D+^t + ■ ■ ■ . (10.32) 

The amplitudes Df- contain the LECs Ci,C2,C3 [out of which we need Ci]. On 
the other hand, these amplitudes can be expressed in terms of the threshold 
parameters (see, e.g., Ref. [238]), 



^o'o = 


= 47r(l + r/) a++ , 


^^0 = 


At(1 I 'n^ > '^°+ 1 /)+ \ n+ 1 ^n+ >■ 


- ^TTi^i 1 r/j < ^2 1 "o+ 1 '^1- 1 -'^1+ r ' 


^0^1 = 


= 27r(^ + a+_ + <(2 + 3r/)l, 



(10.33) 

where r/ = M.,^/mp and a^"^, hf^ stand for the scattering lengths and for the 
effective ranges. Consequently, given an algebraic expression for the amplitude 
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LECs 


0{p') 


0(p3) 


0{p') 


Cl 


-0.9 


-1.2 


-1.2 


C2 


2.6 


4.0 


2.6 


C3 


-4.4 


-6.1 


-6.1 



Table 6: The central values of the LECs Ci, C2, C3, extracted from the experimental 
data at different chiral orders (in units of GeV~^). 

D+(g^, t), one may solve for the constants ci, C2, C3 in terms of the experimentally 
measured threshold parameters and other strong LECs, whose contribution enters 
D^{q'^,t) first at order p^ and is small. Moreover, using D~^{q'^,t) calculated at 
O(p^), O(p^), 0{p'^), ■ ■ ■ , one may extract the couplings Cj from data at the same 
order - the differences between the numerical values corresponds to the residual 
quark mass dependence in these constants^^ (note that, while doing this, we treat 
Dqq, D^q and Dq^ as a fixed input and do not expand in the variable 77). 

In Table 6 we summarize the central values of ci, 02,03 at different chiral 
orders. As an experimental input, we have used Koch's values for the thresh- 
old parameters [5]. Other LECs which contribute to Z}+(g^,t) at order p^ are: 
0{p^) LEC /a from the pion sector and the fourth-order pion-nucleon LECs 
Cl ■ ■ -eg [238]. In the calculations we use I3 = 2.9 and set the finite part of the 
constants Cj to at the scale /i = nip. We have checked that the dependence of 
Cl, C2, C3 on these LECs is rather weak, so even a large uncertainty here does not 
affect the final result. From this Table one readily sees that, going from 0{p^) 
to 0{p^) leads to substantial quark mass effects in all q, see also Ref. [146]. On 
the other hand, ci and C3 become stable already at 0{p^), suggesting that the 
procedure is convergent. 

Table 6 contains only central values. The estimate of the uncertainties is more 
subtle. The quoted experimental uncertainties on the threshold parameters are 
rather small. However, since in actual calculations one is truncating amplitudes 
at a given order in ChPT, the uncertainty in the LECs is set to reflect the 
uncertainty due to the higher orders as well. Bearing this in mind, we assign the 
following error estimate to our central value of Ci : 



Cl = -1.2 ± 0.3 GeV*^ [order /] 

At order 0{p'^) we assign an asymmetric error: 

Cl = -0.9+°J GeV"^ [order p^] . 



;i0.34) 



;i0.35) 



Here we take the same uncertainty as in the recent analysis of LECs of the O(p^) 
Lagrangian [240]. In this paper, the author compiles various determinations of 



^®We thank T. Becher and H. Leutwyler for providing us with the exphcit analytic expression 
for the amphtudc D^{q^ , t) up to and inchiding the fourth order terms. 
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the couplings known in the hterature, including the information about the sub- 
threshold vrA^ amplitudes and the A^A^ potential [146, 236, 241-245]. The resulting 
values of ci range from —0.7 GeV~^ to —1.4 GeV~^, which is in reasonable 
agreement with our finding (the sign of the uncertainty for ci in Ref. [240] must 
be reversed). On the other hand, our error bars are generous enough to include 
the shift of the central value when going from O(p^) to 0{p^). We therefore 
expect that the result remains stable with respect to higher-order quark mass 
effects. 

Finally, we note that our former result Ci = —(0.93 ± 0.07) GeV~^ at O(p^) 
[43] includes the uncertainty of the scattering lengths only. Even though that 
result is compatible with the above one within the error bars, we prefer to use a 
more conservative error estimate here. 

Electromagnetic LECs: /i and /2 

Next, we turn to the electromagnetic constants /i, /2- The quantity /i occurs 
in the chiral expansion of the nucleon mass and in elastic pion-nucleon scatter- 
ing 7r^p{n) — »■ 7T^p{n). The electromagnetic part of the proton-neutron mass 
difference is given by the constant /2 at leading order in the chiral expansion, 

-e^F^f2 = {mp-mnr'^. (10.36) 

Here we disagree with Ref. [246, Eq. (12)] by a factor of 2. Numerically, we use 
{mp - m„)<='" = (0.76 ± 0.3) MeV [247], or 

/2 = -(0.97 ± 0.38) GeV"^ (10.37) 

We are left with the determination of /i. The sum rrip + rrin contains the com- 
bination e^(/i -|- /a) - the constants /i and /s can thus not be disentangled from 
information on the nucleon masses. We may consider rup + rUn as a quantity 
that fixes /s, once /i is known. Therefore, elastic pion-nucleon scattering is the 
only realistically accessible source of information on /i. In principle, one may 
consider combinations of amplitudes that vanish in the isospin symmetry limit, 
and determine /i from those. The combination 

has this property. The tree graphs of X start at order p^ and contain /i - that 
one may try to determine hence from here. Of course, one is faced with a problem 
of accuracy: in order to determine X, one needs to consider the difference of two 
large numbers, quite aside from the fact that the cross section n'^p -^ ix^p is not 
known experimentally^^. It remains also to be seen whether a combination of 
experimental data and lattice calculations can provide a reliable estimate of /i. 



^^For a proposal to measure the elastic scattering of neutral pions, see, e.g., Ref. [248]. 
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In the absence of precise experimental information, we can i) rely on order- 
of-magnitude estimates, or ii) consider model calculations. As to order-of- 
magnitude estimates, we follow Fettes and MeiBner [249] and write 



FVlM^^^m,, 



a 

2^' 



or 

|/i| ~1.4GeV-\ (10.39) 

because /i is due to a genuine photon loop at the quark level (we divide by 
2tt rather than by 4-7r [249] to be on the conservative side). This estimate also 
confirms the expectation [41] that |/i| has the same size as I/2I, see Eq. (10.37). 
In Refs. [250,251], O(p^) LECs have been estimated in the framework of a 
quark model, with the result (in our notation) 

ci = -1.2 GeV~S 
(/i,/2,/3) = (-2.3 ±0.2, -1.0 ±0.1, 2.1 ± 0.2) GeV^i. (10.40) 

For a recent attempt to calculate /i by using a method related to resonance 
saturation, see Ref. [252]. 

Isospin breaking correction 5i,5[ and data analysis 

Using the values of LECs at O^p^), which are given in Eqs. (10.35), (10.37) and 
(10.39), we find that 

S[ = (-4+1) . 10-2 [order p^] , (10.4I) 

where the central value corresponds to /i =0. We add the uncertainties quadrat- 
ically. The bulk of the uncertainty is due to the one in /i and ci. 

As mentioned above, the calculation of the isospin breaking correction at 
0{p^) in ChPT has been carried out in Ref. [43], see also Ref. [235]. Using the 
estimate Eq. (10.34) for ci, these calculations lead to 

S[ = (-9.0 ± 3.5) • 10-^ [order p^] . (10.42) 

It turns out that, albeit the loop contributions at 0{p^) are sizable due to the 
presence of chiral logarithms, the contribution of LECs at this order is quite small. 
The uncertainty is dominated by the uncertainty in /i and, to a certain extent, 
by the one in ci. We expect that the contributions of the higher-order LECs, 
which are multiplied by additional powers of M^, will be even more suppressed 
and do not alter significantly the uncertainty displayed. At this order, isospin 
breaking in the 7r~p elastic scattering amplitude is still of purely electromagnetic 
origin: the terms with {rrid — my) appear first at 0{p^). 



The isospin breaking corrections in the charge-exchange and to the vr n 
7r~n elastic threshold amplitudes at 0{p'^) are [45,48] 



One may use the above expression to evaluate the correction term to the decay 
width in Eq. (10.6). At this order, one has [45]^^ 

5i = (0.6 ± 0.2) ■ 10^2 [order p^] , (10.44) 

where we take qa = 1.27, and where the error comes from the uncertainty in 
/2. Note that, unlike the amplitude Ac, the isospin breaking part of the charge- 
exchange amplitude at 0{p^) does not contain the quantities /i and ci. In this 
channel, one is, therefore, able to determine the isospin breaking corrections 
to a much better accuracy. Still, achieving an accuracy that is comparable to 
the present experimental precision is not possible without the evaluation of the 
corrections of order 0{p^) (at least). This task should urgently be addressed. 

The results given in Eqs. (10.42) and (10.44) should be contrasted with the 
ones obtained in the framework of potential models, see Eq. (10.7). As antic- 
ipated, these two sets are rather different, because potential model do not, in 
general, fully reflect QCD+QED. 

We now come back to the discussion of Fig. 10, which was presented in the 
beginning of this section. It is seen that presently, the theoretical uncertainty 
in the scattering length a^j^ is much larger than assumed in the potential model 
approach. For this reason, a precise determination of a^l^ from hydrogen data 
alone is not possible, unless one finds a way to pin down /i more precisely. On 
the other hand, as we will show later in this article, a simultaneous analysis of 
T[H and ttD data may allow one to pin down the scattering lengths in a more 
reliable manner. 



11 Kaonic hydrogen 

11.1 The kaonic hydrogen and kaonic deuterium experi- 
ments at DA$NE 

An unique source of negative kaons, which provides important conditions for 
the study of the low-energy kaon-nucleon interaction, has been made available 



Note that in Rcf. [45], 0{p ) expressions for the scattering lengths were used in numerical 
calculations, instead of their experimental values. If we use the aj^_ and Oq^. from e.g. Ref. [114], 
the result changes insignificantly to 5i = (0.7 ± 0.2) • 10~^ [at order p^]. 
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by the DA$NE electron-positron collider in Frascati. The DEAR experiment 
at DA$NE and its successor SIDDHARTA aim at a precision measurement of 
the strong interaction shifts and widths of KH and KD . The final aim is to 
extract the antikaon-nucleon scattering lengths from the measured characteristics 
of these atoms [149-152, 154, 155]. 

In the DEAR experiment the low-momentum negative kaons, produced in the 
decay of the (/)-mesons at DA$NE, leave the thin-wall beam pipe, are degraded 
in energy to a few MeV, enter a gaseous target through a thin window and are 
finally stopped in the gas. The stopped kaons are captured in an outer orbit 
of the gaseous atoms, thus forming the exotic kaonic atoms. The kaons cascade 
down and some of them will reach the ground state emitting X-rays. 

As mentioned before, recent results of the DEAR collaboration [154] for the 
shift and the width of the kaon hydrogen in its Is ground state (see Eq. (2.14)) 
considerably improve the accuracy of the earlier KpX experiment at KEK [153] 
and confirm the repulsive character of the K~p scattering at threshold. It should 
be also pointed out that the DEAR central result deviates from all earlier exper- 
iments [153, 253-255], see also Fig. 12 below. Below we shall discuss the implica- 
tions of these beautiful new measurements for establishing the precise values of 
the KN scattering lengths oq and ai. The discussion is based on Ref. [46]. 

11.2 DGBT— type formulae for kaonic hydrogen 

The non-relativistic theory of KH [46] is almost a carbon copy of the nH case. It 
is, on the other hand, amusing that the existing small differences at the end result 
in a strikingly different picture (this, in our opinion, is another demonstration 
of the power and flexibility of the NRQFT approach). We do not present many 
explicit formulae, because they often look identical in the ttH and KH case - 
we rather concentrate on those properties of these two systems which are not 
the same. For example, there is no need to carry out again the calculation of 
the electromagnetic shift - it is given by the corresponding expression derived in 
the pionic hydrogen case Eq. (10.19) after obvious replacements M^r -^ Mk and 

{r-D - {rl). 

We start to list those characteristic features of the KN system which cause a 
difference to the case of pionic hydrogen. 

i) The only states that are degenerate in mass with the K~p state in the 
isospin limit 5 — > 0, are the states K~p -\- rwy, K^n + m'7, with m^m! = 
0, 1, ■ ■ ■ . It is convenient to use a framework which explicitly "resolves" all 
these states in the non-relativistic theory. On the contrary, the effect of 
other intermediate states, whose mass is not degenerate with that of the 
K~p state in the isospin limit, could be included in the couplings of the 
non-relativistic effective Lagrangian. 
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ii) Breaking of SU{3) symmetry, which is proportional to the quark mass dif- 
ference TUs —rh, is much larger than isospin breaking effects. For this reason, 
we count rris — m as 0(1) in S. Because there are open strong channels be- 
low K~p threshold - e.g. vrS, vr'^A, the real and imaginary parts of strong 
kaon-nucleon couplings count as 0(1). This leads to a completely different 
power counting as compared to the nH case (see also subsection 7): the 
strong shift and width are quantities of order 6^. 

In the non-relativistic theory, where the neutral {K^n) channel is explicitly 
resolved, the couplings are analytic in 5 (there are no corrections that go like 
a/^, cf. with the discussion in section 10 and in appendix C). The expansion 
of a typical 4-particle coupling constant (similar as in Eq. (10.9)) is given 
by 

g. = c/f^ + ac/^ + (m, - m.)^f + 0{5') , (11.1) 

where gl \gl \gl , ■ ■ ■ are functions of the strange quark mass m^, along 
with ?Ti = I (m„ -|- md). 

iii) The Panofsky ratio Eq. (A. 34) in -kH is of order (5~^'^. Numerically, the 
branching ratio into the 727 channel amounts to a contribution of ~ 40 % 
in the total decay width. In contrast to this, in the case of KH , this 
branching ratio counts as 0{5~^). The measured branching ratios into 
the leading A7, S7 channels are of the order of a per mille [256] (the 
theoretical description of this quantity by using chiral Lagrangians [57, 
257] gives a result which is consistent with the experiment by order of 
magnitude). Consequently, the perturbative treatment of the effects due to 
these channels is justified. At the numerical precision one is working, and 
bearing in mind present experimental and theoretical uncertainties, they 
may even be neglected altogether. 

iv) The K^n intermediate state lies in the vicinity the K~p threshold. Due to 
this fact, the well-known unitary cusp emerges in the K~p elastic scattering 
amplitude (see, e.g., Refs. [258-260] and appendix C) which is the source 
of huge isospin breaking corrections (note that the cusp effect is also the 
dominant isospin breaking effect in some other low-energy processes, e.g. 
in neutral pion photo-production off nucleons [261, 262] as well as K ^ Svr 
decays [88-93].). 

v) The A(1405) never appears explicitly in this approach, that is correlated 
with counting (m^ — m^) as a hard scale. The influence of the resonance 
is only indirect and results in a large KN threshold scattering amplitude, 
which then leads to a significant increase of the isospin breaking corrections. 

Due to the presence of the K'^n states, the relative weight of the isospin 
breaking correction in the K~p threshold amplitude is dramatically increased. 
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Figure 11: K^n bubbles in the threshold amplitude of elastic K p scattering. 
Note that the bubbles with K~p in the intermediate state vanish at threshold. 

Consider, for instance, the real part of this amplitude, which contains any num- 
ber of the s-channel bubbles with K~p and K^n intermediate states. Further, 
the contribution of the former vanishes at threshold, and only the neutral loops 
contribute to the threshold amplitude (see Fig. 11). 

Consider now the series of the bubble diagrams shown in this Figure. Since 
the K^n state lies higher than K~p state, the contribution from a single K^n 
bubble, which is shown in Fig. lib, at threshold is real and proportional to 
go = (2/io("^n + Mf^o — nip — Mk))^^"^ = 0{v6), where yUo is the reduced mass of 
the K^n pair. Thus, the isospin breaking correction to the threshold amplitude, 
which stems from the diagram Fig. lib, starts at 0{vS) (note that go = in 
the isospin limit). Moreover, since both real and imaginary parts of the KN 
couplings Qi, which appear in the vertices, count as 0(1), the contribution of the 
diagram in Fig. lib to the real part of the threshold amplitude counts always as 
0(vo), irrespective of the fact whether the neutral loop is real or imaginary (in 
the latter case, the unitary cusp is below threshold and influences the real part 
of the amplitude through coupling to the inelastic channels). 

For comparison, let us consider the counterpart of the diagram Fig. lib in the 
vrA^ case, in the non-relativistic theory where the ir^n state is explicitly resolved. 
It is purely imaginary at threshold, since m„ + M^o < rup + M^. Furthermore, 
the imaginary part of the non-derivative ttN couplings counts as 0{S) since there 
are no open strong inelastic channels in this case (the channel jn contributes at 
0{5)). For this reason, there is no contribution from this diagram to the real part 
of the threshold amplitude at O(v^). Only the product of two neutral bubbles, 
which is shown in Fig. lie can contribute to the real part, but this contribution 
starts at 0{S). To summarize, the isospin breaking corrections in the KN and 
vrA^ amplitudes emerge at 0{V6) and 0(5), respectively - no wonder that the 
former are much larger in magnitude than the latter. 

The crucial observation which enables one to find the way out is the following: 
albeit the isospin breaking corrections which are non-analytic in the parameter S 
are large, these can be expressed solely in terms of those scattering lengths which 
one tries to extract from the experiment. Consequently, the presence of these 
large corrections does not affect the accuracy of the extraction procedure. On 
the other hand, the isospin breaking corrections at 0{6) cannot be expressed it 
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terms of the scattering lengths alone - these in particular contain the LECs from 
the isospin breaking sector of the ChPT Lagrangian. One expects, however, that 
this error, introduced by the 0{6) terms, should be of the order of magnitude of 
a few percent, as in the pionic hydrogen. 

To illustrate the above statement, let us consider the bubble sum shown in 
Fig. 11 in the non-relativistic effective theory, using only non-derivative 4-particle 
vertices. The corresponding approximate non-relativistic K~p — > K~p scattering 
amplitude, which we denote as T^j!]^-, is given by the following expression, 

(0) _ 27r \ (gp + Qi) + goaoai ^ 27r 

where /ic is the reduced mass of the K~p system, and where we have used the 
matching in the absence of isospin breaking to relate various non-derivative KN 
couplings to the pertinent combinations of two scattering lengths oq and ai (the 
error introduced is of order 5). Note also that for our purposes, it would be enough 
to include the contributions from the first two diagrams only. We shall however 
sum up the whole series shown in Fig. 11 - the difference which we introduce 
is again of order 6. The crucial observation is that, albeit the isospin breaking 
corrections in T^ jvfl ^^^ °f order v^, the difference Tc^nr — Tcnr — 0{5), where 
Tc^NR denotes the full non-relativistic amplitude. In other words, the leading- 
order isospin breaking effect due to the neutral loop has been explicitly included 
in the quantity T^^j^- 

Finally, using the matching condition, one may define the relativistic ampli- 
tude, corresponding to the relativistic bubble sum given by Eq. (11.2), 

Ti°) = 2mp 2Mk r}°jin = Svr (m^ + Mk) a, , (11.3) 

which has the same property as its non-relativistic counterpart, namely, the 
difference %. — % = 0{5), where % is the relativistic elastic amplitude for 
K~p -^ K~p. In general, one has 

Tc = rj°) + ^^ (ri°))^ + s% + 0(6) , (11.4) 

which is nothing but the definition of the correction term 6% (the second term 
in the above definition with the imaginary coefficient has been added for conve- 
nience. It starts at order 6.). Note also, that Eq. (11.4) is the generalization of 
the relation 

% = Svr {rup + Mk) - (ao + Oi) + O(v^) . (11.5) 

The main message here is: in contrast to Eq. (11.5), the quantity 5% in Eq. (11.4) 
is of order 6 and therefore can be assumed to be not too large. Note also that 
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Eq. (11.2) exactly coincides with the amphtude introduced in Refs. [258, 259]. In 
our approach we however in addition demonstrate that this modification accounts 
for a// potentially large 0(^/6) corrections to the threshold amplitude. 

The results of our findings can be finally summarized an a modified DGBT- 
type formula for generic ns energy levels of KH , 

A^f-^r„ = - "'f; 3 (rw + ^T.){i-g/^I^TW + -..| + o(^^), 

(11.6) 

where 5Tc = 0{S) and ellipses stand for (tiny) vacuum polarization contribution, 
which we do not display explicitly. 

One expects that the equation (11.6) is much better suited for the analysis of 
the DEAR experimental data than the original DGBT formula. In this equation, 
potentially large (parametrically enhanced) isospin breaking corrections at 0(v5) 
and^^ 0{6 In 6) are explicitly separated from the rest, which is analytic in 6 and is 
assumed to be small. Indeed, in Ref. [46] an estimate of this 0{6) term has been 
carried out at tree level in the SU{3) version of ChPT. The calculations result in 
the effect at a percent level that supports the above conjecture. In the numerical 
analysis of the DEAR data that follow, we shall always set 6Tc = 0. 

Finally, it is interesting to note that Eq. (11.2) enables one to independently 
test the limits of applicability of the method for a given values of the scattering 
lengths a posteriori. Namely, in order to be consistent, the term of order of 6 in 
the expansion of the above amplitude should not exceed a few percent, and the 
following terms must be negligible (see Table 7). This is, however, not the case for 
all input values of scattering lengths available in the literature, see Refs. [46,49] 
for a detailed discussion. 



11.3 Analysis of the DEAR data 

The Figure 12 summarizes the analysis of the DEAR data with the help of 
Eq. (11.6). Here we display the predicted value of the energy shift and width 
in the ground state, using different scattering lengths as an input, see Table 8. 
These predictions are compared with the DEAR measurement, as well as earlier 
experimental data. The comparison enables one to immediately conclude that 
the scattering data, to which the parameters of Refs. [159, 161, 163, 166] are fitted, 
are, in general, not consistent with the DEAR data [46]. This conclusion is valid 
for all input scattering lengths shown in Fig. 12, except those from Ref. [163]. 
For a further discussion on this issue, we refer to Refs. [162, 164]. 
The comparison allows one to conclude that: 



^^The corrections at 0{S\nS), which are referred to as "Coulomb corrections" hereafter, 
emerge ft-om the second term in curly brackets in Eq. (11.6). 





Ref. [159] 


Ref. [263] 




ao = -1.31 + 1.24i 


ao = -1.70 + 0.682 




ai = 0.26 + 0.662 


a^ = 0.37 + 0.60i 


O-cfi 


-0.52 + 0.95i 


-0.66 + 0.64i 


ac,i 


-0.68 + 1.09Z 


-0.98 + 0.66Z 


ac,2 


-0.67+1.15i 


-1.04 + 0.73i 


ac,3 


-0.65 + 1.16i 


-1.04 + 0.75i 


^C^OO ^c 


-0.65 + 1.15i 


-1.03 + 0.76i 



Table 7: Expansion of the K~p — ^ K~p scattering length Oc in powers of go 
(bubble approximation, see Eq. (11.2)). The index n in a^n corresponds to the 
n-th iteration. The results are given in fm. 



Source 


Oo 


ai 


MeiBner and Oiler [159] 


-1.31+Z1.24 


0.26 + Z0.66 


Borasoy, NiBler and Weise, fit "n" [161] 


-1.48 + Z0.86 


0.57 + Z0.83 


Oiler, Prades and Verbeni, fit "Aj" [163] 


-1.23 + Z0.45 


0.98 + Z0.35 


Martin [263] 


-1.70 + Z0.68 


0.37 + Z0.60 


Borasoy, Meifiner and NiBler, fit "full" [166] 


-1.64 + Z0.75 


-0.06 + Z0.57 



Table 8: KN scattering lengths ao and ai (in fm) from the literature. These 
scattering lengths are used as an input in the calculations of the kaon-deuteron 
scattering length. 



i) The corrections due to the unitary cusp, which start at O(v^), are indeed 
huge. Even at the present accuracy level, it is absolutely necessary to take 
them into account. 

ii) The Coulomb corrections that are amplified by In a, are also quite sizable. 
For example, choosing scattering lengths from Refs. [159,263], we obtain 
that the real part of the correction term in the ground state amounts up to 
9% and 15%, respectively. 

iii) The key point here is that one does not indeed need to know the numerical 
values of these large corrections very accurately. Since both the unitary and 
Coulomb corrections depend on the scattering lengths ao, ai only, in the 
numerical analysis, which aims to extract exactly these scattering lengths 
from the data, it suffices to know the algebraic form of this dependence. 
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Figure 12: Isospin breaking corrections in the energy shift and width of kaonic 
hydrogen: different input for scattering lengths Oq and Oi from Refs. [159,161, 
163,166,263], see Table 8. Filled dots correspond to the DGBT formula (no 
isospin breaking), empty dots include the effect of the unitary cusp and filled 
squares, in addition, take the Coulomb corrections into account. 

One may further refine the analysis of KH data, fully including the restric- 
tions from unitarity [49]. We first rewrite Eq. (11.2), 



ao + ai 



2go 



2ar 



1 - q^a. 



flgOl 



1 - goOc 



0. 



;ii.7) 



The quantity Qc is fixed from experiment, by using Eqs. (11.2), (11.6), and ne- 
glecting corrections of order 5. We decompose the scattering lengths into their 
real and imaginary parts. 



ai = xj + iyi , / = 0, 1 ; ac = u + iv . 



:ii.8) 



Solving Eq. (11.7) for qq and requiring that its imaginary part is positive leads 
to a quadratic constraints on Xi,yi: in the {xi,yi) - plane, Oi is expelled from a 
circle determined by u,v and Qq, and similarly for 1 <-;> 2. Explicitly, 



(x/ + -)2 + (y, - r)2 > r^ ; / = 0, 1 
r = [(1 + qouf + qy] {^0^ ■ 



;ii.9) 
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Figure 13: Restrictions set by the DEAR data on the values of the scattering 
lengths Co and a\ (thick solid line). For comparison, we give the scattering 
lengths from Table 8. The dashed line corresponds to the restrictions obtained 
by using KpX data instead of DEAR data. 



Part of this circle is shown in Fig. 13 (note that, bearing in mind the preliminary 
character of the DEAR data [154], we use only central values in order to illustrate 
the construction of the plot. We do not provide a full error analysis). In order 
to be consistent with the DEAR data, both a^ and a\ should be on the outside 
of this universal DEAR circle^''. For comparison, on the same Figure we also 
indicate (much milder) restrictions, which arise, when KpX data are used instead 
of DEAR data. The values of a^ and ai, plotted in this Figure, are again taken 
from Table 8. As before, we see that in most of the approaches it is rather 
problematic to get a value for a^ which is compatible with DEAR (the price 
to pay for this is the presence of a very narrow 7 = 1 pole in the scattering 
amplitudes close to K~]) threshold for the solutions from Ref. [163]). This kind 
of analysis may prove useful in the near future, when the accuracy of the DEAR 
is increased that might stir efforts on the theoretical side, aimed at a systematic 
quantitative description of the KN interactions within unitarized ChPT. 

In conclusion we note that, in our opinion, the present experimental data 
are still not precise enough to clearly distinguish between different low-energy 
approaches to KN interactions which are based on the unitarization of ChPT. 



^"Note that an analysis of the KH data alone docs not allow one to determine both scattering 
lengths ao and a\ separately, even at leading order in b. The reason for this lies in the fact that 
both real and imaginary parts of a/, 
channels below threshold. 



/ = 0, 1 count at 0(1), due to the presence of the open 
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In the future, however, DEAR/SIDDHARTA can pose very stringent constraints 
on theoretical models, provided the measurements are carried out with the an- 
nounced accuracy. 

12 Pionic and kaonic deuterium 

12.1 Introductory remarks and DGBT formulae 

Pionic (kaonic) deuterium is a hadronic atom, made up from 7r~ {K~) and the 
deuteron, which in turn is a shallow bound state of the proton and the neutron, 
held together by strong interactions. Constructing a theory of these bound states 
constitutes one step further in the sophistication of the NRQFT approach. 

The characteristic momentum for the deuteron is 7 = y/rriNEd — 46 MeV, with 
Ed = 2.22 MeV its binding energy. This quantity is still three times smaller than 
the pion mass, but much larger than the typical momenta in hadronic atoms 
~ 1 MeV. Consequently, in a very good approximation one may consider the 
pionic deuterium as a 2-body bound state and apply the machinery that has 
been developed in sections 3, 4 and 5. This program has been carried out in 
Refs. [47, 61] for irD. It can be generalized to the case of KD without any change 
- the interested reader may consult these article for details. The analysis of pionic 
(kaonic) deuterium data directly determines the pion-deuteron (kaon-deuteron) 
threshold amplitude. For illustration, we give the expression for the ttD energy 
shift at next-to-leading order in isospin breaking [47], 

AEf''' -l^i = -2«VdXd |l - 2«/idXd (Ina - 1) + ■ ■ ■ I , (12.1) 

where /x^ denotes the reduced mass of the 7r~d system and (cf. with Eq. (11.4)), 
and 

And = Ac- 27Tiafid{Acf = a^d + A^^^ . (12.2) 

Here, Ac denotes the pion-deuteron threshold amplitude, defined in Eq. (A.31), 
a^d is the pion-deuteron scattering length in the isospin limit and A^^^ stands 
for the isospin breaking correction. The energy shift of KD at next-to-leading 
order is also given by Eq. (12.1), after an obvious replacement A^d —* Ap^d ^^^ 
using an appropriate reduced mass. However, the case of KD inherently differs 
from TT-D in one important aspect. Namely, the estimated ratio of the binding 
energy to the decay width in the ground state of KD amounts up only to ~ 8.6 
(see Ref. [49]) and is much larger in the pionic deuterium. This value can be 
still considered as large enough [49] to justify using the machinery based on the 
Rayleigh-Schrodinger perturbation theory (see section 5), but the corrections 
at next-to-next-to-leading order might be not completely negligible, when the 
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accuracy of SIDDHARTA data is close to the planned one [154, 155]. We therefore 
conclude that it could be interesting to perform - at some point in the future 
- an estimate of these corrections within the non-relativistic EFT in a manner 
described in section 5. 

At the next step of the investigation one has to "resolve" the scattering length 
(threshold amplitude) on the composite object (deuteron) in terms of the under- 
lying hadron dynamics, since for us this scattering length is primarily interesting 
as an additional source of information about the pion-nucleon (kaon-nucleon) 
scattering lengths. This is the most difficult part of the problem. 

Let us first ignore the isospin breaking corrections altogether and start with 
TiD . In lowest order in 7/M^ there exists an universal relation, which relates 
the pion-deuteron threshold amplitude to the threshold parameters of the pion- 
nucleon interactions, 



l + r//2 



'^^d = , , ,^,^ (a^p + a^n) + 0(7/M^) , r] = M^/rriN , (12.3) 



where a,rp and a-i^n denote pertinent linear combinations of the S-wave tiN scat- 
tering lengths, see, e.g., Ref. [124]. The relation for the Kd scattering length is 
completely similar. However the correction term in Eq. (12.3), albeit down by 
the small factor "^/M^ ~ 1/3, turns out to be even larger numerically than the 
first term. The reason for this in the case of vrD is that the first term is chirally 
suppressed. A similar situation emerges also in the case of KD , however for a 
different reason. Here, the kaon-nucleon scattering lengths are so large that the 
multiple scattering expansion seems not to converge anymore (see, e.g., Ref. [264] 
for a nice discussion of this issue). We conclude that, in order to extract useful in- 
formation from experiments on pionic (kaonic) deuterium, an accurate evaluation 
of the correction term in Eq. (12.3) is necessary. 

Existing calculations of pion-deuteron and kaon-deuteron scattering lengths 
have been carried out in different settings. The description of low-energy meson 
scattering on the deuteron is one of the central topics of the potential scattering 
theory and has been thoroughly investigated during decades. It would be abso- 
lutely impossible to cover all this very interesting work here, or even to provide 
a fairly complete bibliography. We cite here only few sources [22,66,68,122- 
130,265-269], which we have consulted on the subject. 

In recent years, calculations based on low-energy effective theories of QCD 
in the two-nucleon sector have started to appear. These calculations enable one 
to extract (in principle) the pion-nucleon and kaon-nucleon scattering lengths 
in QCD directly from the experimental data, without any additional model- 
dependent input. Below, we give a very condensed review of these calculations. 
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12.2 Pion— deuteron scattering 

Recently, there has been a considerable activity [131-138] in the study of nd 
scattering on the basis of the low-energy effective field theory of QCD with 
non-perturbative pions, where multi-pion exchanges are included into the NN 
potential that is later iterated to all orders. The approach was first formulated in 
Weinberg's seminal paper [139] and has recently reached a level of sophistication 
that allows one to perform systematic precision calculations in two-nucleon as 
well as in many-nucleon systems. For a recent review of this approach, see 
Ref. [141]. Alternatively, the calculations of the vrrf scattering length have been 
also performed in the EFT with heavy /perturbative pions [47, 142, 143]. 

It is instructive to briefly compare calculations carried out in the two differ- 
ent settings. The couplings of the EFT with heavy pions include the threshold 
parameters of the ttN scattering (scattering length, effective range, etc). Thus, 
the perturbative expansion in this EFT, which is carried in powers of 'j/M.,^, pro- 
duces the multiple-scattering expansion exactly in a form one is looking for. 
Unfortunately, individual terms in the multiple-scattering series get strongly 
scale-dependent after the renormalization. The scale dependence can be only 
canceled by a (large) three-body contribution, coming from the 6-particle LECs. 
And, since the value of these LECs is completely unknown, this leads to a large 
uncertainty in the multiple-scattering series [47, 142, 143]. 

On the other hand, it has been shown [131, 132, 135, 136, 270] that the scale 
dependence in the EFT with perturbative pions is rather mild and therefore 
the pertinent three-body LECs need not be large. This conjecture is indepen- 
dently supported by using dimensional estimate and resonance saturation for 
these LECs [47, 132] . Physically, the strong scale dependence in the heavy pion 
EFT (large LECs) reflects the increased uncertainty caused by working with a 
smaller energy cutoff, than for the EFT with non-perturbative pions. The price 
to pay for suppressing the size of the three-body force is however that the ex- 
pansion parameters in the latter approach are initially quark masses and not the 
nN scattering lengths. 

The conversion of the chiral expansion for the ird scattering length into the 
form of the multiple-scattering series can be achieved on the basis of the following 
heuristic observation (see Ref. [132]): it turns out that those contributions in the 
chiral expansion, which do not have a counterpart in multiple-scattering series 
(these are the diagrams that correspond to virtual annihilation and creation of 
pions, see e.g. Fig. 14b, c), are numerically suppressed with respect to the di- 
agrams with no mass gap (Fig. 14a), although both emerge at the same chiral 
order. In Ref. [132] this property was formalized by introducing the so-called 
"modifled power-counting." Now, re-grouping the contributions in accordance 
to this new counting, in Ref. [132] it has been shown that up-to-and-including 
fourth order the pion-deuteron scattering length can be expressed in terms of 
the pion-nucleon scattering lengths (calculated at 0{p^) in ChPT) and a remain- 
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Figure 14: Illustration of the modified power counting in the theory with non- 
perturbative pions. The shaded blobs denote the wave function of the deuteron, 
and the solid and dashed lines correspond to the nucleons and pions, respec- 
tively. Diagrams b) and c), albeit being of the same order in ChPT, are strongly 
suppressed as compared to the diagram a), see Ref. [132]. 

der (the boost, or Fermi-motion correction), which turns out to be numerically 
not very large [132] albeit strongly scale-dependent. Moreover, unknown LECs 
emerge first at fifth order in modified counting and are expected to give a small 
contribution. 

If one neglects the nucleon recoil (static approximation), the final expression 
for the pion-deuteron scattering length can be rewritten in a remarkably simple 
form [132], 

+ 2 ^j^ ((a+J' - 2(a„-j2(a++ - a,-J) i- (^-^^ + a,,„, + ■■■, (12.4) 

where (■ ■ ■ )wf stands for various wave function averages (we remind the reader 
that the above expression is obtained under the assumption of exact isospin sym- 
metry). Note that Eq. (12.4) has indeed the form of usual multiple-scattering 
series, known from EFT with heavy pions [47, 143] (or the potential scattering 
theory). The difference however is that the wave functions, which are used to 
calculate the averages, are the wave functions in the EFT with non-perturbative 
pions. In the calculations in Ref. [132] the NLO wave functions with the cut- 
off mass in the interval A = (500 • ■ ■ 600) MeV [271] have been used, yielding 
(l/q2)^^ = (12.3 ± 0.3)M^ and (l/|q|)„f = (7.2 ± 1.0)M2 . The boost correction 
in Eq. (12.4) is aboost = (0.00369 ■ ■ ■0.00511)M-i (the strong scale dependence of 
the boost correction is related to neglecting A-resonance contribution, see below). 
Further, the first term in the series is chirally suppressed and is anomalously 
small. The contribution of the second term (which is quadratic in the scattering 
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lengths) is large, but higher-order contributions become again smaller, so that the 
series are likely to converge. The measured value of the pion-deuteron scattering 
length is [112] 

ReA^J = -(0.0261 ± 0.0005)M-^ . (12.5) 

From Eqs. (12.2) and (12.4) one obtains a relation between aQ_^_ and ao+, provided 
that the isospin-breaking part AA^d in Eq- (12.2) is dropped. This relation is 
shown in Fig. 10, see the band named "Deuteron, Beane et al." In the same 
Figure we plot the bands which correspond to the a,|^_ and ao+, determined from 
the 7iH energy shift and width [118]. The combined analysis of the ttH and the 
nD data finally yields the result displayed in Eq. (10.8), see Ref. [132]. 

We would like to mention that there were additional approximations made in 
Eq. (12.4). In particular, as already stated, the static approximation for the pion 
propagator was used, resulting in averages (l/|q|")wf- The possibility of lifting 
this approximation has been considered in Refs. [47, 137,272], where it has been 
demonstrated that different corrections to the static limit largely cancel each 
other, so that Eq. (12.4) indeed describes the exact result quite accurately. What 
is important, corrections to the static limit can be calculated perturbatively, be- 
cause the multiple-scattering series converges in the pion-deuteron case. Further, 
higher-order contributions (most notably, the so-called dispersive contribution) 
have been calculated recently and shown to be rather small [133] (for another 
estimate of the size of this correction in a different setting, see, e.g., Ref. [140]). 
The boost corrections were recently re-calculated taking into account the explicit 
A-resonance [138]. It has been in particular shown that the inclusion of the A- 
resonance removes the large scale-dependence of the boost correction. However, 
since our main aim here was to demonstrate the general framework for the numer- 
ical analysis of the ttD data, we have refrained from including the results of the 
ongoing work into our final plot shown in Fig. 10. Note also that the error bars in 
Fig. 10 reflect the uncertainty in the coefficients of Eq. (12.4), but not these addi- 
tional approximations or the higher-order terms. In other words, we believe that 
there is still some room left for systematically improving the numerical precision 
in Eq. (12.4) in the framework with non-perturbative pions. 

12.3 Isospin breaking in the ird scattering length 

As already mentioned in section 10, the pion-deuteron scattering length, given 
by Eq. (12.4), is not consistent with new data on pionic hydrogen. Namely, from 
Fig. 10 one immediately observes that the intersection area for new ttH bands 
moves far away from the irD band. However, the correction terms, which were 
mentioned in the previous subsection, are too small to be responsible for the large 
shift of the ttD band, which is needed to reconcile the new irH and irD data. 
Consequently, new mechanisms should be sought that could be a possible source 
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of large contributions to Eq. (12.4). It should be also understood that the only 
loophole left in the multiple-scattering series Eq. (12.4) is assuming isospin sym- 
metry. It is therefore interesting to check whether isospin breaking corrections to 
Eq. (12.4) can bridge the gap with the new hydrogen data [48]. 

Already in 1977, Weinberg has pointed out [273] that the isospin breaking 
corrections to certain pion-nucleon scattering amplitudes could become large, if 
the leading isospin-symmetric contributions to these amplitudes are chirally sup- 
pressed. Unfortunately, Weinberg's statement refers to the scattering processes 
with neutral pions that makes it difficult to verify with present experimental 
techniques [see, however, footnote 17]. 

It turns out, however, that this large isospin breaking correction emerges in the 
quantity Aj^d as well, since the leading-order term is proportional to the isospin 
even pion-nucleon scattering length a^^_ and is thus very small (see Eq. (12.4)). 
Quite surprisingly, such (a rather obvious) phenomenon has not been explored 
so far until very recently [48] . Studies of isospin breaking in the nd system (see, 
e.g., Refs. [126, 134, 140,274]) include effects coming from the virtual photons at 
low energy and/or the particle mass differences in the loops. Numerically these 
effects, which emerge at higher orders in ChPT, indeed turn out to be moderate. 
However, as it is well known (see section 10), isospin breaking in ChPT at leading 
order emerges through the direct quark-photon coupling encoded in the LECs 
/i, /2 as well as due to the explicit quark mass dependence of the pion-nucleon 
amplitudes, which have not been taken into account in these investigations. 

Quite obviously, the leading-order isospin breaking correction to the pion- 
deuteron scattering length is given by (cf. with Eq. (12.3)) 

Re AAd = T^^ {SA, + 6 An) + 0{p^) , (12.6) 

1 +?7/2 

where the quantities 6Ac and 6An are defined by Eqs. (10.31) and (10.43), re- 
spectively. 

In the numerical calculations the same parameters as in section 10 are used. 
Namely, since the isospin breaking correction for A-^d is evaluated at 0{p^) only, 
for consistency we use the value ci = — 0.9lo;5 GeV~^ at order p^, see Eq. (10.35). 
For the same reason, here we have applied 0{p'^) isospin breaking corrections ev- 
erywhere, although for the energy shift 0{p^) result is also known (see section 10). 

At the leading order, the isospin breaking correction to the nd scattering 
length is independent on the deuteron structure and is extremely large 

Re AXd = -(0.0110t°:°°^?) M~' , (12.7) 

or. Re A^7rd/Re a^^^ = 0.42 (central values). Plotting again the bands in the 
(ao+, ao+)-plane (see Fig. 15), it is seen that the correction in Eq. (12.7) moves 
the deuteron band in in the right direction: the isospin breaking corrections 
amount for the bulk of the apparent discrepancy between the experimental data 
on TiH and tiD in Figure 10. 
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Figure 15: Determination of the ttN S'-wave scattering lengths Oq^ and ao+ from 
the combined analysis of the experimental data on the ttH energy shift and width, 
as well as the ttD energy shift. Isospin breaking corrections have been evaluated 
in the EFT framework: see Eqs. (10.41), (10.44) and (12.6) for bands 1), 2) and 
3), respectively. The filled circle denotes the current algebra prediction [275], and 
the dashed box corresponds to the 0{p^) calculation in Heavy Baryon ChPT [236]. 
See main text for more comments. 
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Figure 16: The quantity a+ defined in Eq. (12.8), plotted against Cq^, see 
Ref. [276]. It is seen that, with the isospin breaking corrections evaluated at 
O(p^), the three bands 1), 3) and 4) have no common intercept [276]. However, 
one also concludes that the corrections from higher orders are larger than the un- 
certainty at O(p^) coming from the LECs alone. The distance between the bands 
1) and 2), corresponding to the same nH energy shift and to the isospin-breaking 
corrections evaluated at O(p^) and at 0{p^), respectively, may serve as a rough 
estimate of the higher-order terms. 
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It is useful to comment here on the pecuhar algebraic structure of these isospin 
breaking terms, which allows one set up a more convenient framework for the 
analysis of ttH and ttD data. First, we note that the widths of the two bands 1) 
and 3) in Figure 15 are mainly due to the uncertainties in the LECs ci, /i that 
occur in 5Ac,n- On the other hand, as pointed out in Refs. [48,276,277], these 
two LECs occur in the same combination in the energy shifts of tiH and tiD at 
order p^ in isospin breaking. Introducing the quantity [276] 



ci - 2e2/i , (12.8) 



"+-"°- + 4vr(l + ,)U.^ 

the expressions for the energy shifts in ttH and tiD then contain only a+ and a^^ 
(up to small contributions from the scattering length a^^). One may then display 
the available information in the (a+,a(j'^)-plane [276], see Figure 16, from where 
it is seen that the data on nH and ttD do not lead to a common intersection^^ . 

Taking into account higher order terms in the multiple scattering theory of 
ttD improves the situation [276]. On the other hand, up to now, isospin breaking 
effects have been taken into account only at O(p^) in vrD, albeit certain contri- 
butions are already calculated [133,274]. We illustrate in Figure 16 that these 
higher order terms do indeed matter: the blue band displays the band from the 
energy shift in vrif , provided that terms at order p^ in isospin breaking are taken 
into account (that band corresponds to the blue band in Figure 10, translated 
here to the variable a_|_). We conclude that a similar investigation for tiD is 
urgently needed. A more complete analysis may then lead to an improved de- 
termination of the scattering length a^^ and of the quantity a+, because there 
is an additional consistency constraint between tiH and nD data. On the other 
hand, the problem with the measurement of a^^_ persists: unless one finds means 
to pin down /i more reliably, the uncertainty in a^j"^, which must ultimately be 
determined from a+, is of the order of aQ^ itself. 

For completeness, we now display results of independent analyses of the 
T[H data performed by other groups. 



In Ref. [23] an analysis of the tiH data has been performed and the GMO 
sum rule has been used for a precise determination of the ttNN coupling 
constant. The energy shift and width of irH were taken from Ref. [116] and 
Ref. [114], respectively, and the isospin breaking corrections 

S[ = (-7.2 ± 2.9) ■ 10-2 [43] , ^i = (0.6 ± 0.2) ■ 10^^ [45] (12.9) 

have been applied. The final results for the scattering lengths are 

a^+ + ao+ = 0.0933 ± 0.0029 M-\ ao+ = 0.0888 ± 0.0040 M'^ (12.10) 



^"'^This seems to be in contradiction to the common intersection of the bands 1) and 3) in 
Fig. 15. However, that intersection corresponds to a different choice of the LECs in the ttH and 
ttD energy shifts, and signals an apparent consistency only [278]. 
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• In the second work of Ref. [118], the isospin breaking corrections from 
Eq. (12.9) have been apphed to the recently measured values of the energy 
shift and width given in Eq. (2.8). The result is 

a(j+ = 0.0069 ± 0.0031 M-\ ao+ = 0.0864 ± 0.0012 M'^ (12.11) 

12.4 Kaon-deuteron scattering 

It is expected that the forthcoming measurement of the (complex) energy shift 
of KD by the SIDDHARTA collaboration at DA$NE [154, 155] enables one 
to extract - in combination with KH data - independently both S-wave KN 
scattering lengths ao and ai. We show in this subsection that this indeed appears 
to be a feasible task. 

We start with the kaon-deuteron threshold amplitude Apcdi which is extracted 
from the analysis of KD data. This amplitude must be "unfolded" in terms of 
elementary KN scattering lengths, via multiple-scattering theory. Because the 
multiple-scattering series for the kaon-deuteron scattering length is likely not 
to converge at all, a re-summation is needed. A detailed investigation of the 
problem within an EFT framework has been carried out recently in Refs. [49, 279], 
see also Ref. [280] for a treatment of the same problem in potential scattering 
theory. The multiple-scattering series has been re-summed in Refs. [49, 279] to 
all orders, assuming that the nucleons are static (this is also referred to as Fixed 
Centers Approximation (FCA)), and neglecting derivative KN interactions. The 
resulting expression is given by 

(l + ^) Ar, = £ dr {u\r) + w\r)) Aj,,{r) , (12.12) 

where u{r) and w{r) denote the usual S- and D-wave components of the deuteron 
wave function, which are normalized via the condition J°° dr (u^(r) + u;^(r)) = 1. 
The NLO wave functions in the theory with non-perturbative pions [141,281] 
have been used in the calculations [49]. Furthermore, 

^ . A. + A.H2AA.-BD/r-2BlA.y 

l-A,An/r^ + BlAn/r^ 



with Bl = Al/il + Ao/r), and 

, Mk 



rup 



•^c,n,xfl •^c,n,x,0 1 \^^-^^) 



where Ac,n,x,o denote the threshold scattering amplitudes for K~p —+ K~p, 
K~n — > K~n, K~p -^ K^n and K^n -^ K^n, respectively, see appendix A. 
Retaining only the leading isospin breaking effects which are due to the unitary 
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cusps, these amplitudes can be expressed through the scattering lengths oq and 
ai (cf. with Eq. (11.2)), 

_ |(ao + Oi) + go«o«i 
•^^ - l + f(ao + ai) 

2 (ao — ai) 



1 - "T (oo + ai) 



A = ai + 


■ ■ ■ ; 


i (ao + ai) - iqcaoai 


1 «gc 

^ 2 


(ao + ai) 



;i2.i5) 



where 



gc = y/2ficA , go = \/2/ioA , A = m„ + M^o -nip- Mk , 

^^^= Tlf^' ^0 = ^f^~- (12.16) 

The isospin limit corresponds to gc = ^Zo = ^^. Finally, the quantity 5AkcL 
contains the 6-particle non-derivative LEG, which describes interaction of the 
kaon with two nucleons. A dimensional estimate yields a few percent systematic 
uncertainty due to the presence of this LEG [49]. In the numerical results, which 
are displayed below, this term is neglected completely. 

Remark: We comment on an important difference between the kaon-deuteron 
and pion-deuteron scattering in this respect. In the case of pion-deuteron scat- 
tering, a formula similar to Eq. (12.13) can be derived, see, e.g., Ref. [126]. The 
TiN threshold amplitudes that enter the expression are real, up to small isospin 
breaking effects. On the other hand, due e.g. to the process vrrf —> nn, the imag- 
inary part of the corresponding correction 5A.,^d does not vanish in the isospin 
symmetry limit. It is, therefore, expected to dominate the contributions from the 
first part, as a result of which the imaginary part of the pion-deuteron scattering 
amplitude at threshold cannot provide information on the vrA^ scattering lengths. 
In contrast to this, the imaginary parts of the KN amplitudes in Eq. (12.13) are 
of order one. As just mentioned, it is expected that SAk^ is small as compared 
to the first term, and is neglected here. This explains why kaonic deuterium 
does provide information on the imaginary part of the scattering lengths, while 
pionic deuterium does not. For an investigation of this issue in the framework of 
a potential model, see Ref. [269].) End of remark 



^^Here, we identify Ac,n,xfi in Eqs. (12.13) and (12.14) with the KN amphtudes evaluated 
at the pertinent physical thresholds. In the language of non-relativistic EFT, this amounts 
to neglecting diagrams which describe rescattcring of K on the same nucleon. Including the 
rescattcring contributions leads to the replacement of Axfi by the amplitudes evaluated at 
s ~ {Mk + ?Tip)^ instead of s = {Mf^o + "^n)^- 1^ addition, in the static approximation 
the mass difference A is neglected in the neutral kaon propagator (for comparison, see also 
Refs. [269,280], where the problem has been addressed within the framework of a potential 
model). Numerically, the effect may not be negligible. However, following Ref. [49], we do not 
take it into account at this point. 
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We would like to emphasize that, in order to reduce the systematic error on the 
calculated kaon-deuteron scattering length, it would be very important to extend 
the calculations beyond the lowest-order formula Eq. (12.12) by using the non- 
relativistic Lagrangian method of Ref. [49]. In particular, it would be interesting 
to estimate the effect of lifting FCA, as well as to evaluate the contribution of the 
Z^-waves (derivative interactions). The comparison with the Faddeev approach 
indicates that these corrections may not be negligible (see, e.g., Ref. [264,269]). 

The main goal of the DEAR/SIDDHARTA experiment is to determine indi- 
vidual KN scattering lengths oq and ai from the analysis of combined KH and 
KD data. Since the experimental results on the KD are still absent, in Ref. [49] 
such an analysis has been carried out, using synthetic input KD data. The 
procedure can be schematically described as follows. From Eq. (11.7) one de- 
termines one of the scattering lengths, say Oi. Substituting this expression into 
Eqs. (12.12), (12.13) and (12.15), one arrives at a non-linear equation for deter- 
mining oo with a given input value of Aj^^. 

The numerical solution of the above equation, carried out in Ref. [49], leads 
to an interesting conclusion. It turns out that unitarity and input DEAR data 
for KH impose severe constraints on the possible input values of A^fi, for which 
the solutions for the KN scattering lengths with Ima/ > do exist. The allowed 
region in the (Re^^^, Im^^rf)-plane is shown Fig. 17. It also turns out that the 
allowed values for A^d qualitatively agree with the value of the K^d scattering 
length extracted from the pp — > dK^K^ reaction [282] . Finally, we note that the 
region where solutions do exist is much larger in the case of the KpX input [153], 
than for the DEAR input of kaonic hydrogen [49] . 

The message of the investigation carried out in Ref. [49] is very clear: the com- 
bined analysis of DEAR/SIDDHARTA data on KH and KD is more restrictive 
than one would a priori expect. Moreover, if the corrections to the lowest-order 
approximate result (going beyond FCA and including derivative couplings) are 
moderate, they will not change the qualitative picture shown in Fig. 17. On the 
other hand, they constitute the largest potential source of theoretical uncertainty 
at present. 

We conclude with the expectation that the combined analysis of the forth- 
coming high-precision data from DEAR/SIDDHARTA collaboration on KH and 
deuterium will enable one to perform a stringent test of the framework used to 
describe low-energy kaon-deuteron scattering, as well as to extract the values of 
oq and Oi with reasonable accuracy. On the other hand, a considerable amount of 
theoretical work - related to a systematic calculation of higher-order corrections 
- is still to be carried out before this goal is reached. 
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Figure 17: The region in the (Re^x^^, Im^^^) -plane where solutions for qq 
and ai do exist. NLO wave functions [141,281] with the cutoff parameter A = 
600 MeV have been used in the calculations (for A = 450 MeV, the result changes 
insignificantly). For comparison, we also show the results of various calculations 
of Axd [266, 267, 269, 283] (squares). As we see, none of the calculated scattering 
lengths is located in the shaded area. The figure is the same as in Ref. [49]. 

13 Potential scattering theory 

In this section, we provide a condensed version of our view of the relation between 
EFT and potential model calculations. 

Many calculations of various hadronic atom properties were carried out within 
the framework of potential scattering theory since the seminal work of Deser 
et al. [8]. Potential model calculations (see, e.g., Refs. [22,66-68,122-124,126, 
127,191,229-231,265-269,284-300]) were basically the only theoretical setting 
to address the issue until rather recently, when methods of (effective) QFT came 
in use to tackle the problem. Already the very first example investigated in 
this new framework - the lifetime of the ground state of pionium - revealed a 
qualitative difference of the two approaches: A major contribution of the so-called 
isospin breaking corrections to the energy-level shift turned out to be of opposite 
sign in the two approaches (cf. Ref. [290] and Refs. [25,30]). As mentioned in 
section 10, similar discrepancies were later found in the calculation of the energy 
shift in ttH (Ref. [230] and Refs. [41,43]), see also the recent work Ref. [231]. 
These isospin breaking corrections are of the order of the envisaged experimental 
accuracy, and it is therefore important to reveal the reason for the difference. 
This section is devoted to a clarification of this point. 

We find it most instructive to illustrate the two approaches in the case of the 
strong energy shift of ttH . We display three expressions for the energy shift in 
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the ground state, in chronological order: 



A^f 



NReT 



NRer<l + (3i 



ReT 



0{r 



B 



ib) 



NReTrll 



4^F$ii;(i-«-i)R-^^} 



0(6' 



N 



aVc 



AiimpMTr 



The amplitudes T, T and T^ stand for the elastic 7r~p -^ ir'p amplitude at 
threshold, in specific settings. The formulae (a) , (6) and (c) were derived by 
Deser et al. [8], by Trueman [191] (properly adapted here), and by Lyubovitskij 
and Rusetsky [41], respectively. The first two are potential model ones, whereas 
the last one uses the EFT framework advocated in this report^'^. We now discuss 
these expressions in turn, and start with (a). Here, T is the elastic scattering 
amplitude in the absence of Coulomb interaction. Additional terms on the right 
hand side of this relation were not considered in Ref. [8]. The second relation 
corresponds to an expansion in powers of T/tb, where T is evaluated in the 
presence of the Coulomb interaction. The coefficient /3i depends on the dimen- 
sionless ratio tq/tb-, where ro(rB) stands for the effective range (Bohr radius). In 
order to compare these two formulae, one needs to know the relation between the 
amplitudes T and T. In Ref. [191], the relation [301] 



ReT = ReT <^ 1 



yU^a In a 
AnmpMT^ 



ReT + 0(a) 



(13.1) 



was invoked, which indicates that (a) is a first order approximation of (b) , in the 
sense of the expansion performed in (b). 

We now come to the relation displayed in (c), worked out in Ref. [41], see also 
section 10. Optically, the relation is similar to (a), (b). However, as is worked out 
at length in this report, (c) is based on a very general framework: quantum field 
theory. This framework allows one to include in a systematic manner strong and 
electromagnetic interactions. The amplitude 7^ is of the form 

ReTe = 87r(mp + M^) {a^^ + a^^) + aFi + (m^ - m„)F2 + o{5) , (13.2) 

where the first term on the right hand side is a particular combination of vrA^ 
scattering lengths, and where the coefficients Fi^2 depend on the underlying the- 
ory. In this report, this is taken to be QCD+QED, and the corrections encoded 



^•^As mentioned, there were very many potential model calculations over the last four decades, 
improving (6) in several respects. However, the basic difference between these approaches and 
EFT remained untouched. For this reason, we stick to the evaluation by Trueman for simplicity. 
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by -Fi,2 stand for the so-called short-distance electromagnetic contributions to 
the scattering amplitude. The quantities Fi 2 can be calculated in ChPT, in an 
expansion in the quark masses niu = rrid, in terms of a well-defined set of LECs. 

It is at this stage that the EFT framework and the potential model calculations 
differ: EFT allow one to provide a direct contact between the measured energy 
shifts and the calculated scattering lengths in QCD. For more details on the 
relation between the EFT and potential model framework, we refer the interested 
reader to Ref. [302]. In this reference, it is proved that the relation between the 
strong energy shift in the form (c) also holds in potential scattering, a property 
that goes under the name universality. Universality clearly shows that it is not the 
use of the potential framework that causes discrepancies with the EFT approach, 
but the choice of the short-range potential, which must include isospin breaking 
effects properly. 

In our opinion, for the reasons just outlined, potential model calculations are 
superseded by EFT methods in the case of simple systems like A2n, A-kK, t^H and 
KH . On the other hand, in more complicated systems like ttD and KD , the 
results of calculations performed within potential models may still be useful as a 
hint about the expected magnitude of different contributions. However, in order 
to avoid an uncontrollable systematic error, one would have to address the same 
calculations within the framework of low-energy effective theories of QCD as well. 

14 Summary and outlook 

In recent years, a general theory of hadronic atoms has been developed, as is 
outlined in this review. It relies in an essential manner on effective field theories: 
ChPT, which describes the interactions of hadrons and photons at low energy, 
and non-relativistic effective QFT, that are used to describe hadronic bound 
states. It turns out that the treatment of any hadronic bound system is pretty 
universal within this theory (i.e., independent of any details characterizing the 
constituents) and can be carried out with a surprising ease. The approach is sys- 
tematic: the terms that are neglected count at higher orders in a power-counting 
scheme as compared to the ones that are retained. Another important feature 
of the present approach is the fact that the needed effective field theories are 
built on top of QCD+QED by successively integrating out various high-energy 
scales. Consequently, calculating observables of a hadronic bound state within 
this approach, one ends up with a quantity that is evaluated in the fundamental 
theory of strong and electromagnetic interactions. 

The theory of hadronic atoms presented here merges several fields of theo- 
retical physics. As already pointed out, it heavily relies on ChPT and on non- 
relativistic effective field theories. Presently, the attention is shifted to hadronic 
atoms containing the deuteron or other (light as well as heavy) nuclei. As a result 
of this, methods of few-body physics and in-medium properties of hadrons are 
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needed in addition. 

We believe that the conceptual problems of a theory of hadronic atoms have 
now been clarified to a large extent. In this review, we have described several 
specific applications. There is, however, still room left for future investigations 
on the subject. We have collected in Table 9 our view of the status of the 
theoretical description of the various hadronic compounds. The first column in 
the Table displays the system under consideration, and the second one contains 
the status of the relation between the energy spectra and the pertinent scattering 
amplitudes. The third column displays information about the relation between 
the hadronic amplitudes and the various scattering lengths, whereas the last one 
refers to keywords related to the underlying physics. 

The main points displayed in the Table are the following. 

i) A precise calculation of isospin breaking corrections in pionic hydrogen 
energy shift and width should be performed. This implies calculations done 
solely in ChPT, in analogy with Ref. [43]. In particular, it would be useful 
to evaluate the charge-exchange amplitude at 0{p^) (for the width) and, 
possibly, at O(p^) for both 7r~p elastic and charge-exchange amplitudes 
(energy shift and width). 

ii) Using the results of these calculations, the analysis of the irH data should 
be done anew. In particular, it would be very interesting to update the 
value of aQ_^_ and of the pion-nucleon a-term. 

iii) In irK scattering, the structure of higher order ChPT contributions [105] 
remains to be understood in view of the low-energy theorem Eq. (2.7). An 
issue still to be investigated is the relation between the large/small 5^7(3) 
quark condensate scenario, and the pertinent scattering lengths. 

iv) A substantial progress in the precise quantitative description of pion-deute- 
ron scattering in the chiral EFT would be extremely desirable. This implies 
systematic calculations of various higher-order correction terms, including 
both, isospin conserving [Eq. (12.4)] and isospin breaking [Eq. (12.6)] con- 
tributions, see e.g. Ref. [133]. 

v) We believe that KH and KD , which will be investigated by the SID- 
DHARTA experiment at DA$NE [154, 155], represent the most challenging 
theoretical task at present. Many issues should be addressed in this con- 
text. For example, one has to gain a deeper insight into the incompatibility 
of the DEAR and scattering data. Is it possible to fit all data by using 
unitarized ChPT? What kind of additional experimental input could help 
to critically constrain the parameters of a fit? 

vi) Still in connection with the SIDDHARTA experiment, it would be a major 
breakthrough to present a systematic EFT calculation of the kaon-deuteron 
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Table 9: Status of the theoretical description of various hadronic compounds, 
ordered according to increasing complexity. The first column displays the system, 
the second the relation between the spectrum and the scattering amplitude T, the 
third column concerns the relation between the scattering amplitude T and the 
scattering lengths, and the last column displays issues of the underlying physics. 
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scattering length in terms of the threshold parameters of the KN interaction 
beyond the approximations used in Refs. [49,279]. 

On the experimental side, very precise data on nH are underway [118]. The 
DIRAC collaboration has improved on the uncertainty of the lifetime measure- 
ment of pionium, see Ref. [81], and an update of the central value of a^ — a2 is 
anxiously awaited. A successful measurement of individual scattering lengths via 
excited states, and ttK scattering lengths, would be extremely welcome. We are 
looking forward to precise data on KH and KD from SIDDHARTA. In conclu- 
sion, there are exciting times ahead on the experimental side. 
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A Notation 

In this appendix, we collect some of the notation used. 

A.l General 

The following Coulombic bound states are considered: 

Scalar QED / tt+tt- / n^K^ / n^p / K^p / n^d / R-d . 

In all these cases, at least one of the particles has spin zero. We always attach 
the label 1 to the particle with non-zero spin (if present). The masses of the 
particles are denoted by rrii and 1712- We further define 

mim2 mi 2 

S+ = mi + m2, ^ic = ^ — , ?7i,2 = ^^, (A.l) 

^+ ^+ 

where S_|_ and Hc denote the full mass and the reduced mass of a system of two 
particles, respectively. 

A. 2 Coulombic bound states 

For Coulomb bound states, we use 

En = ^+ - ^T^ , 7n = - , -f = afic, rB = 'j , 

in'^ n 

a~^ = 137.036, (A.2) 

where En stand for the unperturbed Coulomb energies and r^ is the Bohr radius 
for a given state. Further, in the case of the two spin-0 particles, the bound-state 
obeys the Schrodinger equation 

(Ho + Hc)|^„,„(P)) = hn + ^^ |*nim,(P)), (A-S) 

with the wave function given by 

f d^a 

|^„im(P)) = / 7^ ^nzm(q) |P, q) , 

|P, q) = al(r/iP + q)4(r/2P - q)|0) . (A.4) 

The same formulae in case of a particle with spin is 

(Ho + Hc)|^„,,,„(P)) = Un + ^) |^„Z,>n(P)), (A.5) 
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and 

— — ^ (jm|/(m - a)sa) ^„i(m_a)(q) |P, q, o") , 

|P, q, a) = bliviP + q, a)4(r/2P - q)|0) . (A.6) 

In the above formulae, a\ 2 and b\ denote creation operators of particles without 
and with spin. 

The wave functions '^nimi'i) are the Fourier transform of standard Coulomb 
wave functions \&„im(x) in coordinate space, 



^nlm{ci) = I d'^e-'^^'-^nlmi^) • (A.7) 

The normalization is 

j d'^\^>rarn{^)\^ = 1, / 7^ I ^n^m,(q) l' = 1- (A.8) 

Exphcitly [303], 

where Cl^{x) are Gegenbauer polynomials, with generating function 

oo 

{l-2xs + s')-"^ = Y,C':{x)s\ (A.IO) 

n=0 

and Nnim are normalization constants, chosen in conformity with Eq. (A.8). From 
this representation, 

[ H\ |q| ^0 
*nUq)~<^ • (A.ii) 



|q| ) |q| ~^ oo 



The ground-state wave function is 



(q' + 7') 



^ioo(q)= r„. , ^2 - (A-12) 



Due to rotational symmetry, the corrections to the energy levels, which are 
calculated with the use of the above wave functions, do not depend on the index 
m. Therefore, we use everywhere m = and introduce the shorthand notation 

*nz(q) = ^„w(q) , ^n«(x) = ^„io(x) . (A. 13) 
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We often need the wave function at the origin, 

|^no(0)P = ^. (A.14) 

In the bound-state calculations we use the notation 

s„(a) = 2(7^(n)-^(l)-- + lna-lnn), i){x) =T\x)/V{x) . (A.15) 

n 

A. 3 Master equation 

When only the Coulomb interactions are present, the resolvent [z — H)^^ is 
meromorphic in the complex 2;-plane, cut along the positive real axis above the 
elastic threshold. The poles are located aX z = E^. If other interactions are 
turned on, the bound-state poles are shifted from their pure Coulomb values into 
the complex plane. The degeneracy of the energy eigenvalues is lifted as well. 
The (complex) energy shift is given by 

/^E^ij = z^ij - En = (*„,z,|r"'^(^n)|*n;,) + o{5^) , (A.16) 

where the matrix element is calculated between unperturbed Coulomb wave func- 
tions and does not depend on the quantum number m. Further, z^ij denotes the 
new position of the pole, which was shifted from its original position at E^ and 
the pole-subtracted amplitude t"-'-^{z) is defined by Eq. (10.17). In order to sim- 
plify notations, we often do not attach explicit indexes nlj to the pole position 
in the text. In the case of two spin-0 particles, the index j must be suppressed. 
In the main text, we refer to equation (A.16) as master equation. 

A. 4 Energy shift of the atom 

The energy shift is further split into the electromagnetic and strong pieces, as 
well as the term corresponding to the vacuum polarization (whenever electrons 
are present) 



AE^i, = AE:T, + AEIT + 5i, (^AEf - ^ r.^ + o{5') . (A. 17) 



(If both particles have spin 0, the above formulae are modified by merely discard- 
ing the index j). The first term in this expression is given by Eq. (8.2), Eq. (9.6) 
or Eq. (10.19), for the case of 7i~^tt~ , tt^K^ and 7r~p atoms, respectively. The sec- 
ond term, which corresponds to the vacuum polarization contribution, in the case 
of the spin-0 particles is given by Eq. (5.41) and its generalization to the case of 
particles with spin is straightforward. Namely, this contribution does not depend 
on the total momentum j. Further, at the order in isospin breaking parameter 6 
we are working, only the S-wscve strong shift should be taken into account. The 
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strong shift in higher partial waves is suppressed by additional powers of S and 
does not contribute at this order. 

Further, note that if the splitting between the energy levels with the same 
orbital momentum / is tiny, it is convenient to define the averaged level energies 

l+s l+s 

AEra = AT-i Y, (2J + 1) ^^n^, , ^e= Y. ^^J + 1) • (A.18) 

j=\l^s\ j=\l^s\ 

When one speaks, e.g., about the energy of Is or 3p levels in ttH , this average 
is meant. 

A. 5 Threshold ampUtude 

Let us now consider the definition of the threshold amplitude. The scattering 
amplitude for the process 1 + 2^3 + 4 is given by (the trivial term without 
interaction is omitted) 

{picr',p2;out\qia,q2;in) = i(27r)V(pi + P2 - qi - q2)Ta'a{pi,P2;qi,q2) ■ (A. 19) 

The above amplitude describes the scattering in the channel with oppositely 
charged particles, as well as in other channels (e.g. in the "neutral" channel). 
The Condon-Shortley-de Swart phase convention is used. If there are no particles 
with spin, the indexes a', a should be dropped. 

At the order we are working, it suffices to evaluate the scattering amplitude 
Eq. (A. 19) at order a. Starting from Eq. (A. 19), one arrives at the threshold 
amplitude as a result of the following procedure. In a first step, the one-photon 
exchange piece is subtracted in the elastic scattering amplitude of two oppositely 
charged particles, 

Ta'a{Pl,P2] qi, q2) = T„'a{Pl,P2] ?! , ^2) + K'aiPl^ ?l) ^^ ^''{P2, ^2) , (A.20) 

where t = {pi - gi)^ = (p2 - q'2)^, and r^,^(pi,gi) and T''{p2,q2) denote the 
electromagnetic current matrix elements of particles 1 and 2, respectively (if the 
particle 1 has no spin, the T'^i^{pi,qi) is replaced by T^{pi,qi)). This procedure 
is demonstrated in Fig. A.l. In the inelastic channels, there is no one-photon 
exchange diagram and Tf^'cr = Tcr'a- Moreover, it is not needed to remove the 
one-photon piece in the elastic scattering amplitude with at least one neutral 
particle as well, because the pertinent contribution to the spin-nonfiip amplitude 
is not singular at threshold. 

We now list the decomposition of the matrix elements r^,^(j»i, gi). 

i) Spin (pion, kaon) 

T>^{p,q) = {p + qrF{Q^), Q = p-q; F(0) = 1 . (A.21) 
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Figure A.l: Subtraction of the one-photon exchange diagram, see Eq. (A. 20). 
ii) Spin I (proton) 

K'Ap, Q) = ^iP, <y') [I'FriQ') + la^^Q^^^l^ u{q, a) , (A.22) 

Fi(0) = l, F2(0) = «:p, (A.23) 

where Fi and F2 denote the Dirac and Pauh form factors, respectively, and 
Kp is the anomalous magnetic moment 

iii) Spin— 1 (deuteron) 

The decomposition is given, e.g., in Refs. [304,305]. 
In the next step, we define the spin-no nflip amplitude by averaging over spins 

T{pi,P2]quq2) = ^—Y^T^a{pi,P2]quq2)- (A.24) 

(7 

At the final step, we go to the two particle CM frame, defined by pi = — p2 = p 
and qi = — q2 = q, and remove the Coulomb phase from the (dimensionally 
regularized) spin-nonflip amplitude T 

e-^'^'^'^fip, q) = ;^ + ^2 In M + T + 0(|p|) , (A.25) 

IpI f^c 

where n = 2 for the scattering of two oppositely charged particles into two charged 
particles (example: 7r"'"7r~ —>■ tt^tt"), n = 1 for the transition of two charged 
particles into a pair of neutral particles or vice versa (example: tt^tt" —>■ tt^tt^) 
and n = for the scattering of two neutral particles into two neutral particles 
(example: tt^tt'^ -^ ti^ti^). The (infrared-divergent) Coulomb phase is given by 

ec = rif^'~'(j^ - l(r(l) + ln4vr) + ln^') , (A.26) 

IpI \a — 3 2 n J 

where fi stands for the scale of dimensional regularization. 

The equation (A.25) defines the relativistic threshold amplitude T in all cases 
of interest. 
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A. 6 Real part of the threshold amplitudes 

In particular cases, it turns out convenient to introduce the threshold amplitude 
A, which coincides with T (or the real part thereof) up to a normalization con- 
stant. Below we give the normalization for all amplitudes needed (for a general 
scattering process 1 + 2-^1 + 2 the notation A/""^ = 87r(?Tii + 777,2) is used, where 
mi and m2 stand for the mass of the particles 1 and 2, respectively). The ellipses 
stand for the isospin breaking corrections. 

— Pionium 

Scattering channels: 

TT+TT^ — > TT+TT^ (c) / TT+TT^ -^ TT^IT^ (x) 

-^c = 777^ Re% = - (2ao + 03) H , 

A = -i^TT- Re% = (ao - as) + ■ ■ ■ . (A.27) 

The scattering lengths Qq, 02 are denoted by aj], a^ in Ref. [84]. 

— ttK atom 

Scattering channels: 

K+n' -^ A'+TT- (c) / K+TT- -^ J?°7r° (x) 

Ac = AfReZ = -{2a^/^ + a'/^) + --- , 

A = -AfRe%/V2 = ]-{a^/^-a^/^) + --- . (A.28) 

The scattering lengths are normalized as in Ref. [228]. 

— Pionic hydrogen 

Scattering channels: 

pTC^ -^ piT" (c) / piT^ -^ niT^ (x) / nvr^ — > mr^ (0) / nir^ — > n7r~ (n) 

vAc,o,n = A/'Re7;,o,n, ^x = A/'Re7;/V2, 

yi„ ^ ^Q^ ag_|_ + ■ ■ ■ , yix- ^ a,g_|_ + ■ ■ ■ . (^A.zyj 

The scattering lengths are normalized as in Refs. [2,238]. 
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— Kaonic hydrogen 

Scattering channels: 

pK~ -^ pK~ (c) / pK- -^ nK^ (x) / nK^ -^ nK^ (0) / riK^ -^ riK^ (n) 

Ac,0,n = -f^ %,0,n 1 -Ax = N Tx , 

Ac= -{ao + ai)^ , Ao = -{aQ + ai) -] , 

An = ai^ , A = - (flo - fli) H • (A.30) 

The scattering lengths are normahzed as in Ref. [157]. 

— Pionic deuterium 

Scattering channels: 



-K d ^ -K d (c) 

— Kaonic deuterium 

Scattering channels: 
K-d -^ K-d (c) 



A = U% = a^d + --- ■ (A.31^ 



Ac = MTc = aKd + --- ■ (A.32) 



A. 7 Two— particle phase space 

In the text, we use the symbols 

^, ^ AV2(^^^n|^ ^ A(x,y,z)=x2 + y^ + .2-2x2/-2xz-2y., (A.33) 
with En given in Eq. (A. 2). 

A. 8 Panofsky ratio 

The Panofsky ratio is defined by 

P = "("'^ ^ "°^,^ , (A.34) 

evaluated at the ir'p threshold. Its value is P = 1.546 ± 0.009 [212] (P = 
1.546 ± 0.010 [213]). At the accuracy we are working, one does not need to 
consider electromagnetic corrections to the cross sections. P is a quantity of 
order S~^^'^. 
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B Generalized unitarity 

Let H be a (non-hermitian) Hamiltonian. We further write H = Hqr+H/, where 
Ho_R denotes one-particle Hamiltonian, which includes all relativistic corrections. 
Acting on the free two-particle states, it gives 

Hoij|qi,q2) = (g? + g2°)|qi,q2), (B.i) 



where q^ = yM^ + qf, i = 1,2. Next, we define the free resolvent 

Gor{z^) = {z-HoR±iO)-' (B.2) 

and the T-operator 

Tmr{z) = Hj + IljGoR{z+)T^n{z) = Tjvr{z)Gor{z+)IIj + H, , 
T]v^(^) = U\ + U\GoR{z-)T^j,j,iz) = T^j,j,iz)GoRiz-)U\ + U\. (B.3) 
The relation to the T-matrix elements, introduced in Eq. (3.10) is given by 
(pi,P2|TAr/j(2;)|qi,q2) = -(27r)'^5^(pi + p2 - qi - q2) TAr/j(pi, pa; qi, q2) , 

z = pl+pl = ql + ql. (B.4) 

Express now the Hamiltonian through the T-operator 

H, = TNR{z)(l + Gm{z+)T^NR{z] 

Hi = (l + Tj,^(z)Goif(z-)) Tj^^(^). (B.5) 

Subtracting these two equations and using the identity 

Gor{z^) - G^r{z~) = -2m5{z - Hoij) , (B.6) 

we finally arrive at the generalized unitarity relation for the T-operator. 
T^^(^) - T^^r{z) = -2m T\^{z) 5{z - YIor) Tnr{z) 

+ A + T\^{z)G,r{z-)\ (u: - Hi") (l + Gor{z'^)T^r{z)\ . (B.7) 

On the other hand, let us consider the non-relativistic multichannel scattering 
with hermitian Lagrangian. Unitarity condition in this case takes the form 

T,,{z) - Ti{z) = -2^ Y.T\^{z) b{z - Ho^) T^,{z) , (B.8) 

k 

where Tjj(2;) = Tat/? (2) describes the scattering in the elastic channel. From 
comparison of the above two equations it becomes clear that the term in Eq. (B.7), 
which contains H/ — H}, describes the "flow of the probability" into the shielded 
channels i ^ k. These channels are therefore encoded in the imaginary parts of 
the couplings in the Lagrangian Eq. (3.6). 
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C Matching and unitarity 

Here we express of the coupling constants gi and ei - introduced in subsection 10.3 
- through the threshold scattering amplitudes. The couplings gi,ei enter the 
expression for the energy shift and width of ttH at next-to-leading order, see 
Eq. (10.23). Taking the real and imaginary parts of this equation, we get 

r„ = —^lm{g, + A^le,-gl{g-c'iEn))) + o{6'/'). (C.l) 

For the real part, the procedure is straightforward since, as it can be checked, 
Reci = 0(1) that means that the contribution, proportional to Reei, does not 
emerge at the accuracy we are working. In is now clear that the answer - after 
adjusting the normalization accordingly - is given by Eq. (4.36), 

Matching of the imaginary parts of gi and ei is more complicated, because these 
break the naive counting rules in 6. A most straightforward way to address this 
problem is to use unitarity. The key observation is that, although the counting 
rules are modified, the modification at this order comes only from one term, 
namely the contribution of the n^n intermediate state, which can be easily singled 
out. 

An additional (technical) complication, which arises here, is related to the 
fact that we have to consider the unitarity condition at threshold in the pres- 
ence of photons and to subtract the infrared-singular pieces in accordance to the 
definition of the threshold amplitude. Most easily, this problem can be circum- 
vented by using the trick described in Ref. [45]: one analytically continues the 
unitarity condition for the negative values of the CM momentum squared p^ and 
approaches the threshold from below. We shall explain below, how this can be 
achieved. 

At the first step, it is necessary to establish the singularity structure of the 
physical amplitudes, when the threshold is approached from below (note that 
previously we have defined the threshold amplitudes through the limiting pro- 
cedure, when the threshold is approached from above). To this end it will be 
again useful to invoke the non-relativistic effective theory, which gives an ana- 
lytic representation of the amplitude in the vicinity of the threshold in terms of a 
finite number of simple loop integrals. In the following, we discard the transverse 
photons as well as non-minimal couplings of the Coulomb photons completely. 
The reason is that these do not affect the matching of luigi at 0(5^^"^) and Imei 
at 0(5^^/^), which suffices at the accuracy we are working. 
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Figure CI: Representative diagrams that describe the behavior of the non- 
relativistic spin-nonflip 7r~p -^ 7r~p amphtude near threshold, see Eq. (C.5). The 
vertex correction corresponds to Vc and the bubbles with and without Coulomb 
photons give Jc and Be, respectively. The ellipses stand for diagrams with two 
and more pion bubbles. These do not contribute to the matching condition at 
the accuracy we are working. The vertex corrections in the final state are not 
shown. 



We start with the non-relativistic amplitude for the 7i~p elastic scattering near 
threshold and remove the one-photon exchange contribution. In this manner, 
one obtains the counterpart of the amplitude Tjsfji, defined in the scalar case by 
Eq. (4.32). In our case, the scattering amplitude contains the spin-nonfiip as well 
as spin-flip parts. In the CM frame it is given by 



Ta'aiP, q) = (^a'aF{p, q) + Z (T„>^ " [q X p] G(p, q) 



(C.3) 



where p and q denote the CM 3-momenta for the outgoing/incoming pion- 
nucleon pair. For determining the coupling constants gi, ei it suffices to consider 
the spin-nonflip part of the amplitude and to project out the S'-wave component 



^(P,q) = ^I]T<,.(p,q), Fo(p') = ^y" rfcos^F(p,q) 



(C.4) 



where 6 is the angle between the vectors p and q. 

Discarding all but Coulomb photons, one is left with only a few diagrams that 
determine the behavior of -Fo(p^) at very small momenta, see Fig. C.l. At order 
a one obtains 



i^o(p') 



1 - 2V;) 91 - 4p2ei + {gi - 4p^e,)\J, + B,) 



(C.5) 



Note that, in order to keep the expressions as transparent as possible, we did not 
further expand Eq. (C.5) in 5 and p^, retaining some higher-order terms as well. 
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The loop integrals, entering Eq. (C.5), are given by 

Vc = -r-^^/-<A-i(r'(l) + ln4vr) + lln-^P' 



+ 0(|pU-3), 
Jc = -^{-pY' + 0i\p\^d-3), 

Be = 5,-^ln^ + 0(|pU-3), 



where p^ < 0. Further, the quantity Jc corresponds to the single loop with the 
proton and the charged pion, no photons. The quantities Vc and Be are similar 
to the quantities defined in section 4 and correspond to the vertex correction and 
to the internal exchange of the Coulomb photon, respectively. 

In Eq. (C.5) we now multiply both sides by (1 + 2Vc) (remember that the 
quantity Vc is real below threshold) and calculate the imaginary part. Using 
explicit expressions for the loop integrals, we see that the following is valid below 
threshold. 



2K) ImFo(p') = (ho + /I'oln -^ j + |p| (h + K In 



/^c 



.2 



+ |p|2(/,2 + /,'2ln-^)+o(p2), (C.7) 



where the coefficients /ij, h'^ can be expressed in terms of (71, ei. We shall further 
assume that Iragi = 0(5^^^), Imei = 0{S~^^'^) (these assumptions will be verified 
a posteriori). To determine the constants gi,ei from the matching, we need 
following relations, which can be established from Eqs. (C.5), (C.6) and (C.7), 

ho = lm{g,+glBc) + o{6^/^), h2 = -Umei + 0(5''^^) . (C.8) 

Let us now use the matching condition, which relates the imaginary part of the 
quantity Fo(p^) to the (one-photon exchange-removed) relativistic spin-nonfiip 
amplitude, defined in accordance with Eq. (A. 24). We denote the latter quantity 
as Tc, where the subscript "c" corresponds to the elastic channel 7r~p -^ ir^p. In 
the CM frame the matching condition reads 

1 /"i ('1 + 2^) 

(l + 2K)ImFo(p2) = - / dcose \ % ImT,(p,q), (C.9) 

2 J_i 2w^{p) 2wp{p) 
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where W7r(p) and Wp(p) are the relativistic energies of the charged pion and the 
proton in the CM frame, respectively. Further, one may replace T), by the full 
spin-nonfiip amplitude Tc, since the one-photon exchange piece is real. 

We further use unitarity in the relativistic theory to evaluate the imaginary 
part of Tc near threshold. Since the threshold is approached from below, the 
elastic contribution from the n~p intermediate state, as well as 7r~p plus any 
number of photons do not contribute. Then, at the order of accuracy we are 
working, only the contributions from nir^ and wy intermediate states should be 
retained. Substituting the result in Eq. (C.9), we finally get 

(1 + 2Vc) ImFo(p') = P*(p')^(p') + p;ip')H,{p') + ■ ■ ■ , (CIO) 

with 

^ X'/'{s{p'),mlM'^o) AV^(s(p^),m^,0) 

s(p2) = {wM+Mp)f, (c.ii) 

and H{p^), H^{p'^) denote the pertinent angular integrals containing the scatter- 
ing amplitudes into various intermediate states. We do not display the explicit 
expressions for H{p^) and H^{p'^) here. The threshold behavior of -Fo(p^) is 
determined by 

P*{p') = P*(o) + (/(o))'p' + o(p^) , p;ip') = p;{o) + oip') , 

H{p') = (^ko + K\n^^+0{\p\), H,{p') = H,{0) + O{\p\).{C.12) 

Note that at this order in 6 there are no logarithms in H^{p^). Various coefficients 
count as 

p^(o) = 0(6'/'), (p*(o))' = o(ri/2) , p;(o) = o(i), 

ko = 0(1), H^{0) = O{6). (C.13) 

The main property of the above representation is that the expansion coefficients 
in iJ(p^), iJ^(p^) are not enhanced in (5-counting to the order we are working. 
At this order, ~ 5~^^'^ behavior in ei comes solely from the derivative (p*(0))'. 
The matching therefore gives 

Im{g, + gfBc) = p%0)ko + p;{0)H,{0) + o{6'/^) , 

-4Imei = {p*{0))'ko + o{6'^/^). (C.14) 

From the above expression one may readily verify the power counting Im(7i = 
0{5^^^) and Imei = 0(S~^^'^). Further, for the particular combination of the 
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low-energy constants gi and ei, which enters the expression for the decay width 
Eq. (C.l), we get 

l + ^y{-ll)ko + o{6'/'). (C.15) 

In the above expression P denotes the Panofsky ratio which, at the order of 
accuracy we are working, is 



a{n p ^ 772) 



0(5^/2) _ ^Q_ig) 



^ ^ a[7r p^Tx^n) ^ p*iO)ko ^ ^^^y^ 



thr Pm^-yiO) 



Further, note that the net effect of the effective-range term at this order is to 
merely shift the argument of the phase space factor to the correct bound-state 
value p*^ = p*i--fl), since v^s^tI) = K + 0(5^). 

In order to evaluate the width at the required accuracy, it suffices to determine 
the coefficient feg. To this end, we recall that the quantity -ff(p^) is defined as 
an angular integral over IT^P = JReT^} + |lmTj,.} , where Tj;(p,q) denotes 
the scattering amplitude for the process 7i~p -^ ii^n. Thus, at lowest order, ko 
must be proportional to {Re 7^} . Further, the imaginary part of the amplitude 
gives the correction at 0{6) to this result. To find this correction, one may invoke 
unitarity once again and make sure that at this order only n^n intermediate state 
contributes. We skip all details here and display the final result. 



^0 = Tin 7TT i^eTA <l+ { r \, , Re To \+ oi5) . 



(C.17) 



From this equation we get 



Im {g, + glB^ + 47^6,) = ^^^^^ (l + i) ^^ |l + (p^(O)A)'} + o{5''') , 

(C.18) 

where we have used Eq. (A. 29) in order to rewrite the above expression in terms 
of the real amplitudes A^ Ax and ^o- Finally, using Eqs. (C.18), (C.l) and 
(10.28), we arrive at the complete expression for the decay width, displayed in 
Eq. (10.26). 

Remark: Although the modification of the counting rules, which was consid- 
ered in this appendix, seems a bit complicated at a first glance, it is in fact just 
a remnant of the unitary cusp effect. This effect has nothing to do with photons 
and can be established in the purely strong non-relativistic theory. Consider, for 
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example, the two-channel model like the one used for pionium (without trans- 
verse photons). The coupling constants Ci, 02,03 in this model obey the usual 
counting rules. Construct now an effective one-channel Hamiltonian, integrat- 
ing out the neutral channel and expressing new effective couplings in terms of 
ci, C2, C3 and the neutral bubble (see, e.g., Refs. [36,42]). The reader is invited to 
check himself that exactly the above modified counting rules for the new effective 
couplings emerge in the one-channel theory. End of remark. 
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